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Abstract. The observation of quark and gluon jets has played a crucial 
role in establishing Quantum Chromodynamics [QCD] as the theory of 
the strong interactions within the Standard Model of particle physics. 
The jets, narrowly collimated bundles of hadrons, reflect configurations 
of quarks and gluons at short distances. Thus, by analysing energy and 
angular distributions of the jets experimentally, the properties of the 
basic constituents of matter and the strong forces acting between them 
can be explored. In this review, which is primarily a description of the 
discovery of the quark and gluon jets and the impact of their obser- 
vation on Quantum Chromodynamics, we elaborate, in particular, the 
role of the gluons as the carriers of the strong force. Focusing on these 
basic points, jets in e + e~ collisions will be in the foreground of the 
discussion and we will concentrate on the theory that was contempo- 
rary with the relevant experiments at the electron-positron colliders. 
In addition we will delineate the role of jets as tools for exploring other 
particle aspects in ep and pp/pp collisions - quark and gluon densi- 
ties in protons, measurements of the QCD coupling, fundamental 2-2 
quark/gluon scattering processes, but also the impact of jet decays of 
top quarks, and W ± ,Z bosons on the electroweak sector. The pre- 
sentation to a large extent is formulated in a non-technical language 
with the intent to recall the significant steps historically and convey 
the significance of this field also to communities beyond high energy 
physics. 



1 Introduction 

Quantum Chromodynamics [QCD], the theory of the strong interactions within the 
Standard Model of particle physics [1121314] , describes the building blocks of strongly 
interacting particles, like proton, neutron and many others, and the forces acting 
between them. The fundamental building blocks of these particles are the spin-1/2 
quarks q |5l6j . which come in three families. Their masses cover a large range [7]. 
The three lightest quarks u, d, s weigh only a small fraction of the proton mass, the 
charm quark c just about the proton mass while the two heavy quarks b, t weigh more 



* In Memoriam: Hans Joos 



2 



than 5 and 180 times the proton mass, respectively. Baryons, like proton and neutron, 
are composed of three quarks qqq, while mesons, like pions, are bound states qq of 
quark-antiquark pairs. 

Since the spin and the spatial S'-wave functions of the lightest baryons are sym- 
metric under the exchange of quarks, the Pauli principle demands the quarks to be 
labelled by new charges, called colours, which discriminate between the three compo- 
nents of the baryons [819110) . Rephrased within the SU(3)c symmetry group for three 
colour degrees of freedom, the colour wave function is antisymmetric. This three- 
fold antisymmetric combination of colours renders baryons non-coloured, i.e. they 
are white SU(3)c singlets; summing up the quark colours symmetrically in mesons, 
these hadrons are white too. By reducing the lifetime of the neutral pion by a factor 
3 2 = 9, the three-colour extension reconciles the prediction of the quark model with 
the experimental measurement, a crucial point in establishing the colour charges. 

Equivalently to the electric charges in electrodynamics, the colour charges of 
the quarks can serve as sources for force fields, which bind the quarks within the 
mesons and baryons [TT]. Eight such gluon fields g are predicted by the non-abelian 
gauge group SU(3)c- Like the photon field they are vector fields with spin = 1, but 
in contrast to the photon they carry colour charges themselves, mediating colour 
flips of quarks by absorption or emission. This theory, Quantum Chromodynam- 
ics, is theoretically described by a non-abelian Yang-Mills gauge theory [T2]. Gluons 
couple to each other, giving rise to three- and four-gluon interactions. These self- 
interactions of the gluons have profound consequences for the QCD coupling. While 
virtual fermionic quarks render the vacuum colour-diamagnetic, the large contribu- 
tion of virtual bosonic gluons renders the vacuum finally colour-paramagnetic. Thus, 
in contrast to the electric coupling, the QCD coupling decreases with decreasing dis- 
tance and the quarks and gluons become asymptotically free [314] . Quarks and gluons 
therefore interact weakly at short distances while the strength of their interactions 
grows with increasing distance, suggesting the permanent confinement of particles 
carrying non-zero colour charges |13) . 

Quarks can be seen in the scattering of electrons or neutrinos off nucleons. The 
final-state pattern of these processes reveals that the leptons scatter off point-like, 
nearly massless spin-1/2 constituents of the nucleons which carry the electric and 
weak charges of the quarks. Gluons inside nucleons, which do not carry electric nor 
weak charges, manifest themselves only indirectly. Half of the momentum of fast 
moving nucleons cannot be accounted for by the spectrum of the quarks alone, and 
it must be attributed to gluons as flavour- neutral constituents [13]. In addition, the 
quark spectrum is modified by gluon bremsstrahlung if the energy of the impinging 
leptons is raised from low to high values |15) . 

However, QCD suggests another method to unravel its basic constituents. As a 
result of asymptotic freedom, quarks and gluons move as quasi-free particles, called 
partons [TB], at short distances. When these coloured objects are accelerated in scat- 
tering processes, they develop bremsstrahlung cascades of narrowly collimated gluons 
and quark-antiquark pairs, which finally transform to equally well collimated hadrons 
at distances at the colour confinement radius of about 1 fm [10 -13 cm]. Thus, the 
quarks and gluons at small distances map themselves into jets of hadrons at large 
distances. Since the quanta associated with the confinement forces are soft, their im- 
pact on the energies and momenta of the jets is small so that the configurations of 
high-energy quarks and gluons at short distances are truly reflected in the energy and 
angular distributions of the jets. Since these jets can be observed experimentally, the 
properties of quarks and gluons can be determined experimentally by jet analyses, 
such as their spins, flavour and colour charges, and their basic interactions. 

It should be stressed here that the field of jet physics and QCD owes a great deal 
of gratitude to the development and successful operations of high energy colliders, in 
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particular, electron-positron colliders. Starting from SPEAR at SLAC, which started 
the physics runs in 1972 and had a maximum beam energy of 4 GeV, the subsequently 
built e + e~ colliders DORIS (physics start 1973; maximum beam energy 5.6 GeV) 
and PETRA (physics start 1978; maximum beam energy 23.4 GeV) at DESY, PEP 
(physics start 1980; maximum beam energy 15 GeV) and SLC (physics start 1989; 
maximum beam energy 50 GeV) at SLAC, TRISTAN at KEK (physics start 1987; 
maximum beam energy 32 GeV), and LEP (physics start 1989; maximum beam energy 
104.6 GeV) at CERN, saw the main jet activity and detailed tests of QCD. The results 
from these machines are the primary focus of this review and are discussed in the first 
six chapters. However, in a long epilogue, described in chapter 7, entitled jets as tools, 
we have discussed some selected results related to QCD and jets, which have come 
out from the electron-proton collider HERA (physics start 1992; maximum e + /e~- 
beam energy 27.6 GeV and maximum proton energy 920 GeV) at DESY, and the 
hadron colliders Tevatron (physics start 1987; maximum p and p energy 980 GeV) at 
Fermilab and finally the LHC (physics start 2010; maximum proton beam energy so 
far 3.5 TeV). It is not our mandate to discuss the technical aspects of these machines, 
which will take us too far afield from the main focus, namely historical development of 
jets and QCD from the theoretical and experimental points of view. For the interested 
readership of this review, the high energy machine related aspects are summarised 
concisely in Reviews of Particle Physics by the Particle Data Group (PDG). Many of 
these colliders, in fact all the e + e~ colliders, are no longer working in particle physics, 
and for these we refer to the 1996 PDG review [T7], while for the others to the 2010 
PDG review [7]. 

Quite early, the final states in e + e~ annihilation to hadrons had been predicted 
to consist [primarily] of two jets evolving from a quark-antiquark pair produced in 
the collision process [18119120] : 

e + e~ — !> qq^ 2 jets . (1) 

Experimental evidence for these quark jets was first provided at the e + e _ collider 
SPEAR [21 22 by demonstrating that the hadrons in the final states were not isotrop- 
ically distributed but accumulated near the central event axis specified by the mo- 
menta of the quarks [33]. At PETRA energies (12 < y/s < 46.6 GeV) the two jets 
could be recognised without any sophisticated analysis, cf. FigfT] (left-side frame). 
Angular distributions and charge analyses finally proved the jets to be associated 
with spin- 1/2 quarks indeed. 

First indirect evidence of gluons was provided by the PLUTO collaboration at the 
e + e~ collider DORIS 21] from the decay T(1S) — > ggg. However, as T(1S) has a mass 
of 9.46 GeV, significant non-perturbative contributions had to be taken into account. 
PLUTO used their 2-jet data below the resonance to extract the q* —> hadrons frag- 
mentation and used this to estimate also the fragmentation g* ~ > hadrons. With this, 
their analysis was in agreement with the expectations from the underlying process 
T(1S) — ¥ ggg. Gluon jets were later discovered unambiguously at the e + e~ collider 
PETRA [25126127128] running at higher energy (typically 30 GeV). A 3-jet event from 
the very early PETRA data [33] is shown in Fig. [1] (right-side frame) . Such events had 
been predicted theoretically |30j for configurations in which the quark pair produced 
in e + e~ annihilation radiates a hard non-collinear gluon: 

e + e~ — > qqg —> 3 jets . (2) 

This bremsstrahlung mechanism is characteristic for gauge theories and it is familiar 
from electrodynamics where charges accelerated in collision processes emit photons, 
as in electron-positron scattering e + e~ — > e + e~j, for example. Bremsstrahlung glu- 
ons in QCD which transform to hadron jets generate characteristic patterns in the 
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Fig. 1. Observation of (a) 2-jet final states in electron-positron annihilation to hadrons: 
e + e~ — > qq — > 2 jets (TASSO [25]); and (b) 3-jet final states in gluon bremsstrahlung off 
quarks in e + e~ annihilation: e + e~ — > qqg — > 3 jets in the TASSO detector [29| . 



final states which allow one to prove the existence of gluons: With increasing energy 
the primary quark jets become narrower; the events become flat and "Mercedes-Star- 
like" (V-shaped); and finally three well separated jets emerge. Detailed comparison 
of the event structure with the underlying theory (QCD) required apart from the 
perturbative (hard) processes also modeling of the non-perturbative (soft) features of 
the quark and gluon fragmentation. Hence, event generators, incorporating the per- 
turbative and non-perturbative aspects of QCD were necessary to relate the emerging 
jet distributions to the predictions derived from gluon bremsstrahlung in QCD. Ref- 
erences [31132] illustrate the early use of such event generators, taking the form of 
Monte Carlo simulations to match the theretical calculations with the experimental 
measurements, which were state-of-the-art tools at that time, and which helped in es- 
tablishing the properties of the quark and gluon jets. To avoid any confusion, Monte 
Carlo in the present context is a numerical computational technique to calculate 
multi-dimensional distributions of a process in which events are generated randomly 
but weighted to reflect the underlying dynamics. Dedicated experiments at PETRA 
and PEP and theoretical progress in the 80's greatly consolidated jet physics and 
led to quantitative tests of QCD. A more modern view of the use of Monte Carlo 
programs, in particular, their role as tools in hard hadronic collisions can be found in 
recent reviews, for example [33] . 

The program to establish QCD in studying quark and gluon jets was naturally 
continued at the e + e~ collider LEP, see, for example, [33], where the increased energy 
could be exploited to measure the gluon self-interactions in multijet events, 

e+e~ -> qqq'q', qqgg -> 4 jets , (3) 

with the production amplitudes dominated by the qqgg states, which included the 
virtual gluon splitting, e.g. g* — > gg. By measuring energy and angular distributions 
of these 4jet-events the colour charge of gluons could be determined, the crucial 
element for generating asymptotic freedom in QCD. Correspondingly, the variation 
of the quark/gluon coupling could be examined for a large range of energies, though 
experiments at PEP, PETRA and TRISTAN had already confirmed the running of 
a s (Q 2 ) in agreement with the renormalisation group (RG) equation. 
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The quark/gluon jet phenomena were also indicated at the pp collider ISR [35j . 
before high-energy jets were unambiguously isolated at the SppS [36] ■ Since then, 
jets in hadronic collisions have become precision tools in not only testing QCD and 
the electroweak physics at the highest available energies (such as at the Tevatron and 
the LHC), but also in searching for phenomena beyond-the-Standard-Modcl (BSM), 
such as dijet resonances and quark substructure. By the same token, jet phenomena 
observed at hadron-hadron and lepton-hadron colliders have provided fundamental 
information on the quark and gluon densities (parton distribution functions) of the 
proton. We shall review this towards the end of this paper, but for now concentrate 
on the general development of jet physics and QCD which took place in the context 
of e + e~ colliders. 

Quark and gluon processes at short distances can be treated, due to asymptotic 
freedom of QCD, in perturbative expansions for the weakly interacting fields. There- 
fore the basic short-distance processes as well as the evolution of the quark/gluon 
cascades are well controlled theoretically. However, the final transition from the 
quark/gluon configurations to hadrons is a long-distance process governed by large 
values of the QCD coupling which cannot be treated in perturbation theory and 
which, so far, cannot be analysed rigorously. Instead, QCD-inspired models have been 
developed which parametrise the transition phenomenologically. Two alternative ap- 
proaches have been worked out in detail. In the first picture a quark moving out of 
the short-distance regime generates a string-like gluonic flux tube which breaks up re- 
peatedly by spontaneous quark- antiquark creation when its length approaches about 
1 fm. This mechanism generates a jet of collimated hadrons with energy and direction 
corresponding to the initial high-energy quark |37] . Gluons had been treated analo- 
gously |31l32j . or they were assumed to generate kinks, local depositions of energy 
and momentum in the strings stretched between quarks and antiquarks [38 . Alterna- 
tively in cluster fragmentation, after splitting all final gluons in a quark/gluon cascade 
to quarks and antiquarks, qq pairs with low invariant masses transform to hadronic 
resonances which eventually decay to the standard low-mass mesons and baryons [39^ . 

After the important work of Ref. [3D] (see, also [41142143] ). numerous methods 
have been proposed, with steadily increasing refinement, to reconstruct the jets ex- 
perimentally. One class consists of algorithms based on sequential jet recombination. 
Particles are sequentially combined if their distance in momentum space falls below 
a pre-set minimum. Typical examples are the JADE algorithm |44j . where the dis- 
tance is defined by the invariant mass of pairs, developed later to algorithms based on 
transverse momenta k t . A second class is built by cone algorithms in which particles 
belonging to pre-defined cones are grouped into jets. Originally introduced to reg- 
ulate singularities in infrared and collinear quark-gluon configurations [45146] . they 
have been developed to a standard method in hadron collider analyses. 

The original jet analyses at PETRA were based on independent-jet fragmenta- 
tion [31132] . providing a valid tool for reconstructing the quark/gluon configurations 
at small distances in e + e~ annihilation. Subtle effects observed later in the hadron 
distributions between jets, were interpreted as string effect [38147] . or explained al- 
ternatively by additional soft gluon radiation with angular ordering [48]. PYTHIA 
HMD], HERWIG [ST] and SHERPA [S3] are modern versions of Monte Carlo pro- 
grams which are used in present jet analyses. 

The connection of jets with QCD has been extensively treated in the literature 
under theoretical and experimental aspects, see e.g. [53]- [55]. This review will sum- 
marise the basic concepts of jet physics, intended to describe how jet physics has 
been exploited to establish QCD as the non-abelian quark/gluon gauge field theory 
of the strong interactions. Addressing also communities outside high energy physics, 
the review is presented mostly in a non-technical language, giving a qualitative ac- 
count of theoretical and experimental developments which have dramatically changed 
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the earlier picture of the strong forces in particle physics. In doing this, we have in- 
cluded some landmark measurements in a chronological order as they were reported. 
The same remark applies to the discussion of the theoretical aspects, and we have 
emphasized only works which were contemporary with the discoveries. The picture 
now is based on a few fundamental principles summarised succinctly in Quantum 
Chromodynamics . 

The topics on which we concentrate are the non-perturbative and perturbative 
elements of quark/gluon jets, including experimental and phenomenological methods 
to define the jets. Early evidence and indirect indications of quark and gluon jets in 
e + e~ annihilation to quarks at SPEAR and T decays to gluons at DORIS will be 
reviewed. In the central core of this paper, we will describe the theoretical basis of the 
discovery of gluons in the three-jet events at PETRA and the measurement of their 
properties. The picture will be completed with LEP. Finally we will demonstrate in a 
few examples how jets can be used as tools for measuring other parameters and funda- 
mental processes of QCD, the gluon content of nucleons, QCD Rutherford scattering, 
etc., but also how to exploit jets for identifying electroweak W, Z and Higgs bosons, 
top-quark physics, and search for new phenomena, in particular possible substructure 
of partons. Such problems have been addressed at HERA and the Tevatron, and they 
will play an important role at the LHC. However, despite discussing some of the most 
recent measurements in jets and QCD, this is not a review of the up-to-date theoret- 
ical advances. We have included some of these topics to introduce the readers to the 
vast areas of particle physics research in which jets and QCD have branched out, but 
emphasize that this article aims primarily at providing a historical perspective. 

This paper is organised in 8 sections and the main topics discussed are as follows: 
Fragmentation properties of quarks and gluons (section 2), discovery of quark jets 
at SPEAR and the first application of perturbative QCD to derive the 2-jet cross 
section in e + e~ annihilation (section 3), gluon jets in T decays and the basic partonic 
process T — > ggg (section 4), jets in QCD and at PETRA and PEP (section 5), jets 
and QCD studies at LEP (section 6), jets as tools, with applications in Deep Inelastic 
Scattering Processes, 77 collisions, and hard hadronic collisions at the Tevatron and 
the LHC (section 7). A brief summary (section 8) will conclude this review. 



2 The fragmentation of quarks and gluons 

Quarks and gluons move, due to asymptotic freedom of QCD, as quasi-free particles 
at short distances of the order of 10 -15 cm (10 -2 fm) in the femto-universe. When 
these coloured objects separate to more than of the order of 1 fm, confinement forces 
become effective which bind the quarks and gluons in hadrons. The hadronisation 
proceeds through the formation of jets in high energy processes which is driven by 
two dynamical mechanisms. These mechanisms can be explicated most easily in e + e~ 
annihilation to hadrons, e + e~ — > qq,qqg, cf. FigJ5J (i) Quarks which are suddenly 
accelerated in the production process at short distance and time of the order of 1 /E <C 
1 fm, will radiate gluons preferentially into a cone of small aperture, dN/dO 2 ~ 
I/O 2 . Subsequently the gluons may split into two gluons or quark-antiquark pairs, 
and, repeatedly, quarks and gluons again into quark and gluon pairs, so that the 
original quark fragments finally into a quark/gluon cascade within a narrow cone, (ii) 
When the coloured quarks on the way out of the femto-universe to large distances 
separate to more than 1 fm, a gluonic flux tube of narrow transverse dimensions 
builds up which fragments into ordinary hadrons. Similar mechanisms lead to the 
hadronisation of gluons. In total, the perturbative quark/gluon cascade as well as 
the partons fragmenting non-perturbatively into hadrons generate jets of particles 



Fig. 2. (a) Quark fragmentation to hadrons induced by confinement forces; (b) Quark/gluon 
cascades at high energies in QCD. 



preserving, in momentum and energy, the original kinematic characteristics of their 
parent partons. 

2.1 Quark fragmentation 

When quarks and antiquarks in high energetic processes, like e + e~ — » qq, separate 
from each other, a linear gluonic flux tube is expected to build up, with energy density 
of about 1 GeV/fm and small transverse size, which will confine the two coloured ob- 
jects. For sufficiently large separations R ~ 1 fm, enough energy will be accumulated 
in the flux tube so that new quark- antiquark pairs can be created spontaneously and 
the flux tube breaks up. This expectation is borne out by lattice analyses [SS] which, 
in static approximation, support this picture qualitatively. Beyond the short-distance 
regime, the potential between heavy quarks rises linearly with distance, V(R) = gR 
with a ~ 1 GeV/fm. However, when the distance between the heavy quarks exceeds 
a value of about 1.2 fm, light quark pairs qq are created and the heavy-quark [QQ] 
system breaks up into two separate mesons [Qq] and [Qq] . 

Repeating this break-up process, adjacent quarks and antiquarks coalesce to hadrons 
with small transverse momenta of the order of 350 MeV so that narrow jets of colli- 
mated hadrons are generated |37l38j . If the partition of energies in q — > h[ q gn +q' scales 
with the energy of the primary quark, the number density of hadrons D{z) observed 
with energy z = E h /E q obeys, for a single species, the recursion formula |37j 



with /(C) denoting the break-up probability for fractional energy £. This equation 
states that the primary meson might be the first in rank primary meson, with prob- 
ability /(l — z)dz, or if not, then the first-rank primary meson has left a momentum 
fraction rj with probability f(rj)dr), and in this remaining cascade the probability to 
find z in dz is F(z /rf)dz jr\. Dividing out by dz gives the above equation. The prob- 
ability /(C), which must be determined experimentally, is generally parametrised as 
a polynomial ~ (1 — C) or a s an exponential ~ C _1 (l ~ C) ex P[ — a /C] m string 
fragmentation. 

From this picture two important consequences can be derived. 
(i) Solving the evolution equation generates a pole in z — > 0, most easily seen for 
polynomial probabilities, 





const 



for z — » . 



(5) 



z 



Thus the fragmentation mechanism predicts a large number of soft (low energy) 
hadrons in the jets, i.e. a long constant plateau in the rapidity y — \ogz~ 1 . 
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(ii) Summing up the hadron charges in the jets reflects the charge of the parent 
quark. If u,d,s quark pairs were created in the flux tube spontaneously with equal 
probabilities, the sum would measure the charge exactly. However, since s-quarks 
weigh a little more than u- and d-quarks, the probabilities for spontaneous quark- 
pair creation deviate from 1/3 by a small amount and a small fraction of the charge 
leaks into the plateau: 

{Q q )=Y t e h = e q - 1 . (6) 

h£jet 

In the parton model language [37] one finds 7 rj 0.067, i.e. (Q u ) — 0.60, (Qd) = 
(Qs) = —0.40. The close relation to the ideal values +2/3 and -1/3 therefore allows 
the efficient tagging of the parent quark charges in the jets. In practise, things are more 
involved in perturbative QCD. Jet-charge studies have been undertaken extensively 
at LEP [57], where a lot of data are available at the Z peak and flavour-tagged 
results |58j . distinguishing between the light-quark, charm and bottom contributions, 
have been obtained. 

The light quark fragmentation to mesons can effectively be described by the frag- 
mentation function 

D{z) = {1 + /3)^-(l - zf with /3~ 0.2 , (7) 

for small jet energies ~ 7 GeV. For higher energies QCD predicts a stronger fall-off 
of the spectrum. [For a detailed overview of quark fragmentation to various types of 
mesons and baryons see [59160] .] 

The fragmentation function of the heavy c, &-quarks behaves rather differently. It 
was recognized very early 61 62 that a heavy flavoured meson (containing a charm 
or bottom quark and a light antiquark) retains a good fraction of the momentum 
of the primordial heavy quark. Thus, due to the inertia of the heavy quarks, the 
fragmentation function of a heavy quark should be much harder than that of a light 
hadron. Estimating the size of the transition amplitude by the energy transfer in the 
break-up process Q — > h^Q^ + q, the fragmentation function behaves, for example, as 

m 

D Q (z) ~ — ^ with 6 Q ~ A 2 1 M% , (8) 

l_I_«a_ 

z 1—z 

with A ~ 200 MeV denoting the strong interaction scale. The spectrum develops a 
narrow maximum near z$ ~ 1 — y/tQ. This form describes the essential characteristics 
of the hard spectra of Q-flavoured mesons in the heavy c, b jets with M c ~ 1.5 
GeV and M& ~ 4.5 GeV. For more recent works on the heavy quark fragmentation, 
see [64I65I66I67I68I69] . 

It should be pointed out that at higher energies, where the heavy quarks are 
produced with momenta much larger than the heavy quark mass, one expects im- 
portant perturbative QCD corrections, enhanced by the powers of the logarithms of 
the transverse momenta to the heavy quark mass, which modify the shape of the 
fragmentation functions. These effects can be implemented using the framework of 
an evolution equation which is discussed later in this review. They have to be incor- 
porated in the analysis of data. 

After the discovery of quark jets in 1975 in e + e~ — > qq at SLAC, detailed studies 
in understanding the hadronisation process, and hence the energy-momentum profiles 
of the quark jets, were initiated in 1977 by Feynman and Field [37 . In their approach, 
the initial quarks and antiquarks produced in e + e~ annihilation fragmented indepen- 
dently in a cascade process, q — » q + (q'q 1 ) — > ht q q>) + q', conserving the charge and 







other flavour quantum numbers at each step of this cascade. To determine the energy- 
momentum profile, light-cone variables p = (jj+,p_,pt) were used with p_ <C p+, 
where p± = E±p«. The fragmentation q — > h+q' is then affected through a primordial 
fragmentation function 

/ g "(z) = l-a + 3a(l-z) 2 , z= {E + Pll)h , (9) 

0+P||)q 

with a an adjustable energy-independent parameter, fixed by data. As already dis- 
cussed, this gives rise to an scale-invariant longitudinal energy distribution of hadrons 
in a jet. Heavy quark fragmentation (for example of a charm quark into a D meson) 
is encoded by a different primordial c — >• D fragmentation function, as already dis- 
cussed in this section above. The px-distribution (px is the transverse momentum 
measured with respect to the jet-axis, which can be identified with the direction of 
the fragmenting quark, for the time being) was implemented in terms of a Gaus- 
sian function, characterised by a q ~ 0.35 GeV, also determined phenomenologically: 
9{Pt) — (2cr 2 )~ 1 e~ PT / 2CT 9 . Like the flavour quantum numbers, pr is locally com- 
pensated, implying that an r th -rank primary meson has a momentum ky(r), with 
kx(r) = qiv — QT(r-i)- The striking feature of the Feynman-Field jet is its simple 
algorithm with the phenomenological profile determined in terms of a few parame- 
ters, which provided an adequate description of the non-perturbative aspects of jets 
initiated by quarks. 

2.2 Gluon fragmentation 

The fragmentation of gluon jets follows rules similar to quarks. Two paths had been 
chosen in the analysis of PETRA jets. The properties of g-jets may be described 
either as a nearly flavour-neutral average |31] of u, d and, with reduced probability, 
s-quark jets, or, alternatively, gluon jets may be reinterpreted as a superposition 
of quark and antiquark jet pairs [33] with the spectra derived from the g — »• qq 
splitting function. In any case, the transition from gluons to quarks g — > qq in creating 
the leading particle will soften the gluon fragmentation compared with the quark 
fragmentation, accounted for effectively by raising the power fall-off towards z — > I 
of the fragmentation function of the order of ~ 1.5. 

2.3 Hadronisation Models 

Quark and gluon configurations created at small distances must transform to bundles 
of hadrons due to confinement at large distances. The transformation requires non- 
perturbative mechanisms and, therefore, cannot be treated rigorously at the present 
time. Instead, models have been developed which describe this step at various levels 
of sophistication. 

2.3.1 Independent jet fragmentation 

Gluonic flux tubes, built up when coloured objects separate, may hadronize into a 
jet of collimated hadrons as argued earlier. While the basic picture had first been 
described for quark jets [37], gluon jets can be analysed similarly when the gluons 
are either treated globally as partons [3T] or split into quark-antiquark pairs, frag- 
mented incoherently again [32] . Implementing energy- momentum conservation in the 
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overall event was an unsatisfactory element of the model. Nevertheless, independent 
jet fragmentation has a simple and transparent structure including a small number of 
parameters. The picture could account quite successfully for the essential properties 
of two- and three-jet events in e + e _ annihilation at PETRA and PEP. Thus, it had 
initially been the right theoretical tool for proving experimentally the gluonic nature 
of the third jet in three-jet events. 



2.3.2 String model 

Apart from the different choice of the primordial splitting function /(C), motivated by 
covariance and side-independence of the beginning of the break-up sequence, quark 
jets in the string model [38] are not very different from independent fragmentation 
schemes. However, gluons are incorporated quite differently. They generate kinks 
which locally transfer energy and momentum to the strings stretched between quarks 
and antiquarks. A small number of hadrons is boosted from the segment between 
quark and antiquark jets to the segments between quark or antiquark and gluon jets. 
This string effect has been observed experimentally as reshuffling of hadrons between 
jets, discussed later. 



2.3.3 Cluster hadronisation 

Quite a different hadronisation mechanism is based on colour pre-confmcment |70] . 
Neighbouring coloured partons in cascades arrange themselves in colour-neutral is- 
lands with preferentially small invariant masses. In practise, the quark/gluon par- 
tons in cascades are evolved down to low invariant masses of the order of several 
^Iqcd = 0(200 MeV), where Aqcd is the scale parameter specific to QCD and ap- 
pears in the argument of a s (Q 2 ). Splitting the gluons in the final step into qq pairs, 
neighbouring quarks and antiquarks form the colourless clusters which may finally 
decay to standard hadrons according to phase space rules |39j . The reduced number 
of radiated gluons off heavy quarks and the small number of large invariant masses 
in the tail of the distribution of the colour-neutral clusters can be covered by non- 
perturbative string-type compliments to the cluster hadronisation scheme. 

Based on these schemes QCD event generators have been constructed which de- 
scribe hadron spectra at a high level of accuracy. While the prototypes had origi- 
nally been developed for hadron production in e + e~ annihilation, the event gener- 
ators have been expanded to proton-(anti)proton collisions and complimented by 
programs for lepton-nucleon collisions. The modern versions of PYTHIA |49l50j . 
HERWIG [51], SHERPA [52] and others involve the cascading of quarks/gluons in 
the final and the p/p initial states, and string or cluster hadronisation in the fi- 
nal states. For multijet final states, frequently produced at high energies in collid- 
ers, elaborate techniques have been developed, based on the relation [71] PS(Q 2 ) = 
ME(Q 2 ) x Sudakov factor \Q 2 max — s- Q 2 ], to accomplish smooth transitions between 
quark/gluon parton showers (PS) and well separated multijet final states described 
by fixed-order perturbation theory matrix elements (ME) squared. 



2.4 Inclusive jet measures 

In this section we discuss some inclusive jet measures which have played an important 
role in the quantitative tests of QCD. The first of these is the observable called 
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sphericity which played a central role in the discovery of quarks jets at SPEAR. In 
its tensorial form it is defined as follows [23] : 

Sep = , (10) 



E l Pi 



which can be diagonalised obtaining the principal axes ni , ri2 and ri3 with corre- 
sponding eigenvalues Q\,Qi and Q3. The Qi can be ordered Qi < Q2 < Q3 and 
normalised so that Q\ + Qi + Q3 = 1. Then the squares of the transverse momenta 
are minimal with respect to the axis 113 and the sphericity S is given by 

S = \{i-Q 3 ) = \(Qi + Q2) , (11) 

with the sphericity axis equal to 113. For events with two particles with equal and 
opposite momenta (ideal two-jet event) we have 5 = and S — > 1 for completely 
isotropic events. Because of the normalisation Qi + Qi + Q3 = 1 only two of the 
eigenvalues are needed to characterise an event. For example one can take in addition 
to S the so-called aplanarity Ap, which is 

A ^ ^ . J2i |PiT,out| , 10 h 

Ap = -Qi = -mm = r — . (12) 

z z Z^i Pi 

The aplanarity Ap minimises the transverse momenta with respect to a plane. Events 
with small Ap are almost planar. The jet variables of an event, Qi, Qi and Q3 can be 
plotted inside a triangle as shown in Fig. [21 in which events obtained by the TASSO 
Collaboration at PETRA are shown. In this plot planar events are found in the strip 
with small Ap, 2-jet events have in addition also small S. Alternatively, sphericity 
can be defined as [23] 

2 E^F ' (13) 

Here, pix are the transverse momenta of all produced hadrons in the final state event 
relative to an axis chosen such that the numerator is minimised. 

The method based on the sphericity tensor, first applied to the analysis of the 
3-gluon decay of the Y resonance [75] and to the analysis of qqg events [30]; has the 
advantage that the eigenvalues Qi and the principal axes rii and from this S and 
Ap can be calculated quite easily. Since in these jet measures the momenta enter 
quadratically, the high momentum particles are weighted stronger in the calculation 
of S and Ap. Also, these variables are not invariant against clustering of particles 
and depend strongly on details of the fragmentation of quarks (and gluons) into 
hadrons. This has, for example, the effect, that the sphericity changes if a particle 
momentum splits up by decay, as for example, p° — > ir + ir~ or by fragmentation, for 
instance q — > q' + meson into two or more momenta. Therefore these variables are 
also sensitive to the emission of soft or collinear gluons. 

There exist other variables which are infrared safe (this term is used for observables 
which are free of divergences when calculated in perturbation theory in the limiting 
case of low energy radiation) and which depend on linear sums of momenta. Known 
examples are thrust T and acoplanarity A which are defined by 

T = max^- — - , 14 
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For thrust T, which was introduced in [73174) . the thrust axis n is obtained by max- 
imising the longitudinal momenta piL along it. T varies in the range 0.5 < T < 1.0, 
where the lower limit corresponds to isotropic events and the upper limit to com- 
pletely collinear configurations. In a similar way for spherocity S', defined as [75] 

S = ( « )mm {zM\) ' (16) 

the |piT | is minimised with respect to a unit vector n. It lies between and 1 for 
configurations from collinear to fully isotropic events. Similar to Ap the acoplanarity 
is obtained in such a way that the PiT.out is minimal with respect to a plane. Planar 
event must have small A values. For massless particles A varies between and 2/3. 

Various other jet measures have been proposed: For example a generalisation of 
thrust to three clusters instead of two, called triplicity [75] or jettiness [3D]. A variable 
introduced for the analysis to verify the existence of four-jet events is the variable 
tripodity D 3 [75] , These and other jet variables will be defined explicitly when they 
are used to interpret specific final state data in e + e~ annihilation in later sections. 



2.5 Jet algorithms 

Classifying multi-particle final states qualitatively in jets with high energies is straight- 
forward for a coarse picture. However, when the picture is refined to a high level of 
precision, algorithms must be employed which define the jets in a rigorous way. In ad- 
dition, when experimental measurements are compared with theoretical predictions, 
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the algorithms used in the experimental analyses must conform with the algorithms 
adopted in the theoretical analyses. 

A multitude of algorithms [77] has been developed to describe jets at high en- 
ergies. A few representative examples should characterise the two classes, sequential 
recombination and cone algorithms. Recombination algorithms have been introduced 
originally in e + e~ annihilation, while cone algorithms have been used frequently at 
hadron colliders so far. We shall discuss some of these algorithms later while discussing 
jets in hadronic collisions. 

2.5.1 Sequential recombination algorithms 

The JADE algorithm [44] is a prominent representative for recombination algorithms 
applied in e + e~ annihilation. Particles are clustered in a jet iteratively as long as their 
distance remains less than a prescribed cut-off value. The distance of two particles is 
defined by the invariant mass of the pair: 

y ij =2E i E j {l- C os6 ij )/El m . (17) 

If the criterion yij < y cut is fulfilled, the particles i and j are combined to a compound 
by adding up energy and momentum, for instance, and the recombination continues 
by pairing the compound with the next particle k. The procedure stops after all 
particles are associated with jets. The cut-off value y cu t is generally chosen in the 
range from 1CP 1 down to 10~ 3 . 

The presence of EiEj in the numerator of yij means that two soft particles moving 
in opposite directions can get combined into a single particle in the early stages of 
clustering, which is counter-intuitive to the idea of a jet as consisting of particles 
restricted in the angular dimension. Apart from this, JADE algorithm leads to a 
structure in higher orders of pertrurbation theory which does not allow itself to be 
expressed in a compact resummed form. Technically, this means that the double 
logarithms of the type a" In 2 ™ y cut (n = 1,2,...), which arise in higher orders of 
perturbation theory and which should be resummed to have the correct perturbative 
form in a limited kinematic region, named after Sudakov, are either not present or 
not discernible easily. To rectify both of these shortcomings, the JADE algorithm has 
been improved by substituting EiEj — > mm[Ef, Ej] in the DURHAM algorithm [75]. 
This amounts to defining the distance by the minimal transverse momentum kt of 
the particles in nearly collinear pairs. 

The concept has been transcribed to hadron colliders 79 80 , where the total sub- 
energy is experimentally not well defined, by switching to un-normalised measures 
and replacing the angles between particles by the differences of the rapidities y = 
1/2 \og(E +p z )/ (E — p z ) along the beam axis and the azimuthal angles tfi in the plane 
transverse to the beam axis, 

da = min[p 2 f ,p 2 J] [( W - y,f + - fafj/R 2 , (18) 

with pu(j) denoting the transverse momenta of the particles with regard to the beam 
axis. The jet parameter R is chosen of the order of one. Since the individual quantities 
(Hi {4>i~4>j)i Pu and p t j are all invariant under longitudinal boosts, the distance 

measure dij is also longitudinally invariant. Recombination with the beam jets is 
controlled by the observable diB = p^f, included parallel to the measure dij when 
recombining all the particles to jets with non-zero transverse momenta and beam jets. 
Originally, the power parameter p had been chosen 1 in the k t algorithm [78] and in 
the Cambridge/ Aachen algorithm [5T]. However, clustering involving hard particles 
are favoured by choosing p = — 1 in the anti-kt algorithm [82] . This algorithm, which 
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makes jets grow outwards from hard seeds as intuitively expected, is the preferred 
tool for LHC analyses. 



2.5.2 Cone algorithms 

Cone algorithms had been introduced in QED to tame infrared and collinear singu- 
larities encountered in photon radiation off charged particles. The concept has been 
translated to QCD in formulating the Sterman- Weinberg jets [IS]. Defining 2-jet 
events as hadron configurations in which all but a fraction e of the total energy is 
contained in cones of half-angle 8 around the central event axis, the ratio 

CT 2 32 Ct s 1 1 

— = 1 - — — log - log - . (19) 
a 6 z7r o e 

describes the 2-jet fraction of hadronic events in e + e~ annihilation in the leading 
logarithmic approximation. 

The transition to hadron collisions has been formulated again by adopting the 
definition of distances based on rapidities and azimuthal angles. The clustering re- 
quires a seed location; the 4-momentum of the cluster is determined by summing 
the 4-momenta of all the objects within a distance R = y/(y — y c ) 2 + ((/) — cj) c ) 2 
from the seed (y c ,<j) c ). In one variant, used in the analysis of the Run II Tevatron 
data, the 4-momenta are summed using the so-called E-scheme [53] (this should not 
be confused with the E-scheme in the analysis of the jets in e + e~ annihilation), 
(E,p x ,p y ,p z ) — *^2(E,p x ,p y ,p z )i, and the various variables are defined as 



1 , (E 



PT = \jpl +P 2 y , V = 2 l n {^ZjJ > * = tarr^/fg . (20) 

This differs from the Snowmass clustering algorithm [53] , used in the analysis of the 
Tevatron I data, in which the clustering centroid was defined as the i/^-weighted 
averages of the rapidity y and the azimuthal angle <fr. The cones are either centred 
around seed particles (defined as those particles setting the initial direction and one 
sums the momenta of all the other particles around these seed particles within a 
specified jet measure), an approach which is not infrared safe, or they are defined 
by grouping all particles into subsets such that the subsets correspond exactly to 
pre-defined cones. For further discussion of these and related issues, we refer to the 
works of Seymour [85) and the comprehensive review of jet measures by Salam [77) . 



3 Discovering quark jets 
3.1 Quark jets at SPEAR 

The notion of jets in e + e~ annihilation is closely connected with the discovery of 
Bjorken scaling in deep-inelastic electron-nucleon scattering in 1968 at SLAC. As 
mentioned in the introduction, inelastic electron scattering on protons and neutrons 
at large spacelike (q 2 < 0) momentum transfer and large inelasticity v can very 
well be described in terms of an interaction of the spacelike virtual photon with the 
pointlike constituents of the nucleon, the partons, identified as the u and d quarks 
inside the proton and neutron |86| . The analogous process with a timelikc (q 2 > 
0) virtual photon is e + e~ annihilation into a quark-antiquark pair, e + e~ — > qq, as 
shown in Fig. [4] where q stands for the quarks u, d, s, c, b. As explained already in the 
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introduction, in this simple model the virtual photon from the annihilating electron 
and positron creates a quasi-free quark and antiquark. The occurrence of real quark 
and antiquark particles in the final state is prevented by the the fact that they carry 
non-trivial colour quantum numbers. The quarks and antiquarks transform themselves 
into hadrons with unit probability under the confinement forces, which act at much 
later times t ~ 1 GeV -1 . These hadrons should appear in the final state aligned 
roughly along the momentum direction of the original q and q, so that two distinct 
hadron jets with oppositely directed momenta appear in experiments. This discussion 
mirrors the early "outside- in" picture of jet formation, which was used in the formative 
stages of jet physics. This approach was later replaced by the so-called "inside-out" 
description where quark-antiquark pairs were created out of the vacuum before the 
step of hadronisation. We will discuss the salient feature of this development later. 

This simple quark model [18 1 9120] was supported by the fact that the total an- 
nihilation cross section for hadron production er(e + e~ — > hadrons) is given by the 
square of the quark charges Q / multiplied with the number of colours Nc of each 
quark q 

a(e + e~ —> hadrons) = do = Nc Q 2 t , (21) 

where the sum over '/' is over all active flavours which can be produced at a given 
center-of-mass energy y/s; a is the fine structure constant a ~ e 2 /137. Dividing by 
the cross section for the production of a pair, a(e + e~ — > one obtains 

the famous Drell-ratio R, defined as 

Z ^SJ*!"^ =NcJ2Qh (22) 



cr(e+e~ — > n+u~) 



f 



which has the numerical value 2 (for f — u, d, s), and assumes the values 10/3 (for / = 
u,d,s,c) and ll/3_(for ,/ = u,d, s,c,b), as the threshold for the processes e + e~ — > 
cc and e + e~ — > bb are crossed. A recent compilation of the hadronic cross section 
a[e + e~ — > hadrons) and the corresponding ratio R is shown in Fig. [5] (taken from the 
Particle Data Group [7]) where the various resonances (ft, u>, 0, J/ip, ...) encountered 
in e + e~ annihilation and the transition regions are clearly visible. Away from the 
resonances, the ratio R is almost flat, increasing as a new quark-antiquark threshold 
is crossed in agreement with the values quoted above. Note that the it threshold (at 
around 350 GeV) lies beyond the energies of the e + e~ collider rings operated so far. 

The production of hadron jets in e + e~ annihilation as a signature of the process 
e + e~ — > qq was suggested by Bjorken and Brodsky already in 1970 [53]. However, it 
was not until 1975 that they were discovered experimentally at SLAC's e + e~ storage 
ring SPEAR by the SLAC-LBL Collaboration using the MARK I detector [21] when 
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Vs [GcV] 

Fig. 5. Measurements of the ratio R as a function of the e + e~ centre-of-mass energy y/s 
[From PDG [7]]. 




Fig. 6. The process e + e — > qq — > jets, with the jets defined using the jet-cone angle 5, 
arising from limited pr of the hadrons. 

high enough centre-of-mass energies y/s became available. At low energies, for example 
at the ADONE ring at Frascati or the original DORIS ring at DESY, it was not 
possible to see jets because the jet cones were too broad. This is easy to understand 
if we assume that the transverse momentum px with respect to the jet direction 
(which, theoretically is the momentum direction of the original quark q or antiquark 
q ), which are emitted back-to back in the c. m. system, is limited and that the hadron 
multiplicity < n > increases only modestly with y/s. The jet cone becomes narrower 
with increasing y/s, characterised by the mean half angle < 6 > of the jet cone (see 
Fig. [6]). At y/s = 4 GeV the particle multiplicity is about 6, so that with < px >— 
0.35 GeV the half-angle of the jet-cone < S >~< p T >< n > /y/s ~ 0.53 ~ 30°. This 
shows that at this energy each of the two jets is broader than 60° on average. 

For establishing the jets it is necessary to determine the jet axis along which 
the transverse momenta of the produced hadrons are minimised, In the early work 
of the SLAC-LBL collaboration, the jet axis was defined in terms of the sphericity 
variable defined earlier. In the SLAC-LBL experiment the mean sphericity was found 
to be approximately constant as a function of the total e + e~-energy E c . m . = y/s up to 
4 GeV and then it decreases with increasing E c _ m _. A detailed comparison is shown in 
Fig. El in which the measured sphericity distributions da/dS at E c , m , = 3.0, 6.2 and 
7.4 GeV are compared with the calculated distributions based on a two-jet model 
and the phase-space. At 3.0 GeV there is no distinction between the two models 
and the data agree with both. At 6.2 and 7.4 GeV the S distributions are peaked 
towards low S favouring the jet model. But the S distributions are still very broad. 
This comparison shows quite clearly that (i) the E c . m . must be high enough to see 
the production of jets in e + e~ annihilation, and (ii) that even at the higher E c . m . 
energy range of the SPEAR storage ring, the jet structure is visible only through 
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Fig. 7. Observed sphericity distributions for data from MARK I detector, jet model (solid 
curves), and phase-space model (dashed curves) for (a) E c . m , — 3.0 GeV, (b) E c . m . = 6.2 
GeV, and (c) E c . m . = 7.4 GeV. (From Ref. [21].) 



a detailed comparison with the prediction of an appropriate jet model. Observing 
the jet structure was easier at PETRA energies, where most of the events have a 
two-jet topology, which, because of the higher energy had much narrower angular 
jet-cones. An example of such an event measured by the TASSO collaboration at 
E c .m. = 31.6 GeV, is shown in Fig. Q] (left-hand frame). 

Further tests of the underlying quark structure of the jets in e + e~ annihilation 
were undertaken at SPEAR. One such test is the measurement of the angular distri- 
bution da jd cos 9 of the jet axis with respect to the beam direction. This distribution 
for the production of massless spin 1/2 particles is [57] 

-^-~l + cos 2 tf. (23) 
dcosO 

The first data came from the SLAC-LBL Collaboration at SPEAR. They did the 
measurement with transversely polarised e + and e~ beams available at the lower 
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Fig. 8. Observed distributions of jet-axis azimuthal angles from the plane of the storage 
rings for jet axis with | cos 8\ < 0.6 for (a) E c . m . = 6.2 GeV and (b) E c . m . — 7.4 GeV. (From 
Ref. HI]-) 



cm. energies of the SPEAR ring. With transversely polarised beams the angular 
distribution has the following form 



da 
dJ2 



1 + a cos 2 9 + oP+P_ sin 2 9 cos 2</> , 



(24) 



where 4> is the azimuthal angle of the jet axis with respect to the storage ring plane 
and P+ and P_ are the polarisations of the e + and e~ beams, respectively. The 
measured cf> distributions (averaged over 8) for 6.2 and 7.4 GeV are seen in Fig. [5] At 
6.2 GeV the beam polarisations are P + = P_ = and therefore the <fi distribution 
is isotropic. At 7.4 GeV, where P+P- = 0.5 the characteristic cos2</> behaviour is 
observed. From this measurement at SPEAR, the value a = 0.97 ± 0.14 [21122] is in 
agreement with the expectation for spin 1/2 quarks, a = 1. 

Similar, but less accurate, results were obtained by the PLUTO Collaboration at 
DORIS for P c . m . = 7.7 and 9.4 GeV [55]. Measurements of the angular distribution 
of jets at higher energies were also performed at the e + e~ storage rings PEP and 
PETRA. Although the beam energies were much higher, yielding a much better de- 
fined jet axis, the result a = 1.04 ± 0.16 [SjJ] does not have a better accuracy than 
the SPEAR measurement, which had the benefit of polarised beams. This test of the 
spin 1/2 nature of the quarks produced in e + e~ annihilation is very much the same 
as the verification of the Callan-Gross relation [SU] i n deep inelastic lepton-nucleon 
scattering: F2{x) — 2xFi{x), which is also very well satisfied experimentally. 
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3.2 Sterman-Weinberg Jets 

The existence proof of jets in QCD was provided by Sterman and Weinberg [35]. As 
already noted, they calculated the cross section eT 2 -jet(e, 8) for the process e + e~ — > 
2 — jets, where the jets are defined by two cones in opposite hemispheres with half- 
angular size 8, having all but a fraction e of the total cm. energy. In the general field 
theory context, jets were anticipated due to the Lee-Nauenberg theorem [5T], which 
states that the transition probability in a theory involving massless particles is finite 
provided summation over degenerate states is performed. 

The Feynman diagrams which contribute in order a 8 (Q 2 ) are shown in Fig. [5] For 
small e and 8, and to leading order in a s (Q 2 ) one has 



02 -jet (e, 8) = (t 



l + C F as( - Q ^ ( -41n2eln<5-31n<5- ^- + ~ + 0{e) + 0(8) 

7T V 3 2 



(25) 

where ero is the lowest order cross section given in Eq. (|2~T1) . Cf =4/3 and a s (Q 2 ) is 
the QCD coupling constant defined in the lowest order 

as[Q ) = m 9 M Q> ' (26) 
(33 - 2?i/) In 



where n/ is the number of quark flavours. The terms O(e), 0(8) neglected by Sterman 
and Weinberg are all finite, essentially proportional to the phase space and have 
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been subsequently worked out [92.. Here A is a scale parameter, to be determined 
experimentally, typically of 0(200) MeV. As Q 2 — > A 2 , a s {Q 2 ) -> oo, signaling the 
breakdown of perturbation theory. The above expression for a s (Q 2 ) also states that 
a s {Q 2 ) — > as Q 2 —> oo, implying that QCD is an asymptotically free field theory. 
In the range of Q 2 where a s (Q 2 )/n <C 1, one has a controlled perturbative region. 

Since, in leading order in a s (Q 2 ), the inclusive hadronic cross section for e + e~ — > 
7 — > hadrons is 

a(e+e~ -> 7 -> hadrons) = a Q (l + ^9A) , (27) 
the complement of <72-jet(e,^) is the 3-jet cross section 

a 3 _ jct (e, 6) = g " a (Q2) C f ( 41n2e lng + 3 \n8 + ^- - ~ + 0(e) + 0(g)) . (28) 
7r \ 34 y 

This implies that for typical jet resolutions, a small fraction of hadronic events should 
consist of three-jets. They were found subsequently in e + e~ annihilation at PETRA 
and we shall discuss them later quantitatively. 

Another example of a jet measure which can be used with ease to characterise 
jets is in terms of the invariant mass of the partons j/y emerging from a hard process, 
defined in Eq. ([TT]). Requiring yij > y m ; n > 0, one avoids both infrared and collinear 
singularities in a perturbative calculation. The first of such ?/ m i n -dependent 2-jet cross- 
section was calculated in [S3], yielding (y m m — v) 



&2-jet — Co 



1 + Cf 2^ ( -21n y-31ny + 4t/lnt/-l + y + 0(y) 



(29) 



The 0(y) terms have been derived in [93] . 

The two prescriptions just discussed have been used in the experimental analysis 
of data concerning jets. Thus, for example, the JADE algorithm [53], used mostly in 
the analysis of the e + e~ data at PETRA and PEP, is based on the y m i n -prescription, 
which can be used to classify also muli-jet events. This was subsequently replaced 
by the fc t -jet algorithm, as discussed in the preceding section. The modified form of 
the the cone-prescription is widely used in the analysis of jets in hadroproduction 
processes. 



4 Gluon jets in Upsilon decays 

The T meson first produced in proton- nucleus collisions at FERMILAB 95 96] and 
identified by the T — > /i + /i^decay was later observed as a narrow resonance with mass 
mr — 9.46 GeV and width JY = (40^g 3 ) keV in the process e + e~ — s- T — > hadrons 
[97198199110011011102] . This resonance is a bb bound state and has the quantum num- 
bers J PC = 1 . As the BB threshold lies above mr, the T(9.46) state is predicted 
to decay mainly into 3 gluons (<?) in QCD, the massless colour-octet vector parti- 
cles [103110411051106110711081109] in complete analogy with the decay of orthopositro- 
nium into 3 photons [llOj . While T(9.46) — > ggg is the dominant decay mode, with 3% 
probability it can decay also into a photon and 2 gluons, T(9.46) — > jgg. Average en- 
ergies of the three partons were measured as (E\) ~ 4.1 GeV for the most energetic of 
the three gluons, with the other two having energies (E2) — 3.4 GeV and (£3) ~ 2.0 
GeV, respectively [111] , in approximate accord with the lowest order QCD matrix 
elements. However, only the fastest of the three partons yielded a collimated jet of 
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Fig. 10. Processes contributing to hadronic final states in the T(9.46) resonance region: 
(a) from direct decays of the X, (b) from the T vacuum polarization and (c) from the non- 
resonant continuum. The fj, + fi~ final state can be produced from the X vacuum polarization 
(d), and from the continuum (e). (from Ref. Ill ). 

hadrons and its detailed phenomenological profile was studied by PLUTO [1111112] , 
Phenomenological models, which included the lowest order QCD matrix elements and 
the fragmentation of the partons (quarks and gluons) , were found to be in conformity 
with a number of inclusive measurements undertaken by PLUTO. 

The contributions to the multi-hadron events from the T mass region originates 



from three sources [111] , as shown in Fig. 1101 The first row in this figure shows 
the direct decay of the T (a), decay through the vacuum polarisation (one-photon 
decay) (b), and the non- resonating continuum (c). Denoting the corresponding cross 
sections as a on (cross-section in the T(9.46) energy range), <r vp (cross-section for the 
T vacuum polarisation), and a oS (cross-section for the non- resonant continuum), the 
cross-section for the T(9.46)-production with direct decay is <7 dlr = cr on — cr oS — a vp . 
Since for final states, a 'direct' production does not exist, the term er vp can 

be obtained by scaling the /^-pair production on and off-resonance to the hadron 
production level. This yields: 

a on _ a oS 

er dir = cr on - er off - er vp = er on - cr off - a oS MM MM . (30) 

Using (cr™ - a°^)/a°^ = 0.24 ± 0.22 [M] and the number of events in the two energy 
regions, the Y direct decay cross section is obtained. This is evaluated as a function of 
sphericity S. The results are shown in Fig. ITTIa. b, c, separately for the off-resonance 
data, the data at the T resonance and the subtracted distribution for the T direct 
decay. These experimental results are compared to the two-jet model based on the 
Feynman- Field model, already discussed (dash-dotted line in Fig. [11] a) and the 
predictions of the phase-space model in Fig. [TT] c (dashed line) together with the 
three-gluon decay model (solid line). The off-resonance data are well described by the 
two-jet model in agreement with the earlier findings at SPEAR. The distribution for 
the direct decay is shifted to larger sphericity values and is best reproduced by the 
three-gluon decay model. 

The T meson is produced at rest. Therefore, the scaled momenta Xi = 2ki/mr 
obey the relations (xi = |xi|): xi + x 2 + x 3 = , X\ + x 2 + x 3 = 2. 

Another possibility to describe the configuration of the final state uses the angles 
between the gluons. For massless gluons the relation between the Xi and the angles 
(I, ■:>: •<•, 2sin8 > 
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Fig. 11. Differential sphericity distributions and the sphericity angular distributions. The 
dash-dotted line in (a) represents the two-jet model. The dashed and solid lines in (c) 
represent respectively phase-space and the three-gluon decay models. The dash-dotted line 
in (d) is proportional to 1 + cos 2 9 and the solid line in (f) to 1 + 0.39 cos 2 8 (from Ref. I24J ). 



This relation allows one to characterise the final gluon configuration in the corners 
of the Dalitz plot (x% = X2 = X3 = 2/3 is the "Mercedes-Star" -like configuration, 
x\ = X2 = 1, X3 = is a two-jet configuration with the third gluon perpendicular to 
the direction of the first two and x\ = 1, x% = X3 = 0.5 is the configuration, where the 
fastest jet recoils against the two others with x-i — X3). The momentum distribution 
of the gluon as calculated in leading-order (LO) QCD is |105I106I107I108I109| 

;a£s - W~- Kfl - «■)' + »fr - + *8fl - *»)') ■ (3D 

The above cross section formula is the basis for Monte-Carlo model calculations men- 
tioned above. In these models the hard cross section for Y — > 3g must be supplemented 
with the additional fragmentation of the 3 gluons into hadrons. 

To compare the decay T — > 3<? with vector gluons as follows from QCD, also models 
with scalar gluons have been considered. The momentum distribution corresponding 
to scalar gluons was derived in |113) leading to the result that they peak at the corners 
of the Dalitz plot and have zeros in the middle of each boundary. In contrast, vector 
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Fig. 12. Corrected experimental distribution in | cos#| from the decays T — > hadrons by the 
PLUTO collaboration, where 9 is the angle between the thrust axis and the beam axis. The 
curves are theoretical distributions for T-decay into vector (solid curve) and scalar (dashed 
curve) gluons, respectively (from Ref. [111] ). 



gluons populate nearly uniformly the Dalitz plot. As the majority of the events have 
one gluon with very low energy, a 2-jet structure is expected for scalar gluon theories 

urn 

The distribution of the gluon direction in space is essentially determined by the 
QCD theory |i05ll06ll07ll08ll09] . For vector gluons QCD predicts 

W(cos 6) ~ 1 + 0.39 cos 2 9 , (32) 

where 9 is the angle between the most energetic gluon and the momentum of the 
incoming initial electron (see Fig. [TT]) . Scalar gluons would give rise to the angular 
distribution |113) 

W(cos8) ~ 1- cos 2 6> . (33) 

The PLUTO collaboration has measured a number of observables to strengthen 
the 3-gluon interpretation of the hadronic T decay. The test of vector gluons versus 
scalar gluon has been done using the angular distribution in cos 9, where 9 is the 
angle between the thrust axis and the beam momentum. Theoretical distributions 
predicted for T decay into vector and scalar gluons, respectively, are shown in Fig. [12] 
compared with the PLUTO measurements |lllj . The data clearly prefer the vector 
gluon decay. Similar conclusions have been reached by the LENA collaboration on 
the basis of their measurements [115] . 

Further comparisons with the 3-gluon model presented in [111) consist of topo- 
logical tests with the jet variable thrust, defined earlier, and the variable triplicity T3 
[75] . an extension of thrust to 3-jet configurations, where the final state particles are 
grouped into 3 classes with total momentum P(C/), 1 — 1,2, 3. The values of triplicity 
vary between T 3 = 1 for a perfect 3-jet event and T 3 = 3y/3/8 for completely spherical 
events. PLUTO data have been analysed also in terms of the fractional energies Xi 
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Fig. 13. Experimental distribution of the reconstructed reduced energy x\ of the fastest 
triplicity jet for the T-direct and off-resonance data taken by the PLUTO collaboration 
compared to Monte Carlo calculations for various models (from Ref. [111] ). 



of the jets, where the jet axes are obtained from the triplicity analysis. If the three 
jets would be completely separated in space, the fractional energies would be inde- 
pendent of the fragmentation of the gluons and would depend, in the lowest order 
perturbation theory, only on the QCD matrix element W(x\, x%, £3) given above (see 
Eg. !?!]) . In Fig.[T3]thc projections of the two-dimensional histograms, spanned by the 
axes X3 and (xi — £2)/ v3 5 on the x\ axis is shown (this is also the distribution of 
the most energetic triplicity jet). The prediction of the 3-gluon Monte Carlo model is 
compared to the data and impressive agreement is obtained, whereas two versions of 
the phase-space Monte Carlo model fail to do so. (See [112] for a recent reappraisal of 
the PLUTO experimental analysis). All this information taken together demonstrates 
that the T-direct decay data are very well reproduced by the 3-gluon model while all 
the other models disagree with T-direct data. These findings were further backed up 
later by the CLEO collaboration at the CESR storage ring at Cornell [116] and the 
ARGUS collaboration at the DORIS ring at DESY [TT7]. 
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5 Jets in QCD and at PETRA and PEP 

To put the contents of this chapter in historical perspective, we would like to intro- 
duce the main detectors which played an important role in the development of jets 
and detailed tests of QCD at PETRA and PEP. In doing this, however, we will be 
very brief and refer the interested readers to the review by Gidal et al. |118j . which 
is a compendium of the properties and performance characteristics of the major high 
energy physics detectors in that epoch, and the review by Lynch |119j . In alphabeti- 
cal orders, these detectors were CELLO, JADE, MARK-J, PLUTO and TASSO (all 
located at the PETRA e+e" rings at DESY, Hamburg), DELCO, the High Reso- 
lution Spectrometer HRS, the Magnetic Calorimeter MAC, MARK II, MARK III, 
and the Time Projection Chamber TPC (all located at the PEP e + e~ ring at SLAC, 
Stanford). As already mentioned in the introduction, PETRA started data runs in 
1978 with the maximum beam energy of 23.6 GeV and PEP came a little later in 
1980 having the maximum beam energy of 15 GeV. These detectors were involved in 
measurements for almost a decade, ending their runs as LEP started taking data at 
higher energies. 

5.1 Jet-like distributions from the weak decays of heavy quarks 

The process e + e~ — > qqg leads to p^-broadening of the quark jets, leading even- 
tually to three-jet topologies as the centre-of-mass energy increases. There is an- 
other source of pT-broadening in e + e~ annihilation due to the production of a heavy 
quark-antiquark pair e + e~ — > QQ, and the subsequent weak decays of the heavy 
quarks/hadrons. For the centre-of-mass energies available at the e + e~ colliders PEP 
and PETRA, the heavy quarks whose production and decays had to be correctly 
taken into account were the charm- ant icharm (cc) and bottom-antibottom (bb) pairs. 
Sampling the theoretical predictions of the top quark mass in the PEP and PETRA 
era, most guesses put it around 10 - 15 GeV [1201121] . Hence, the production of a top- 
antitop pair was widely anticipated at these colliders [122 123 , and their characteris- 
tic jet topologies were worked out in the context of the Cabibbo-Kobayashi-Maskawa 
(CKM) 6-quark model 124 125 . However, as subsequent developments showed, there 
were no top quarks to be seen at PETRA and PEP (or at LEP). Thanks to the 
Fermilab-Tevatron [126] . the top quark has a measured mass of about 173 GeV. 

The event topology in e + e~ annihilation is sensitive to the onset of QQ threshold. 
The data in the center-of-mass energy in the range 9.4 GeV < ^fs < 17.0 GeV 
was analysed [127] in terms of the measures of jettiness, (S) and (1 — T), which 
showed a clear step as the BB threshold is crossed. The data taken by the PLUTO 
collaboration [H] at 9.4 GeV at DORIS, and at 13.0 and 17.0 GeV at PETRA by the 
PLUTO [T2"g] and TASSO [T2"9l collaborations were well described by a theoretical 
Monte Carlo [130. taking into account the production processes e + e~ — > cc and 
e + e~ — ¥ bb, with the subsequent non-leptonic decays c — > sud and b — > cud, following 
the CKM theory of weak decays. 

The effects of heavy quark production and decays above their respective thresholds 
on the jet distributions are taken into account by a three-step modifications of the 
light quark pair production and subsequent fragmentation |130j . The heavy quark 
mass enters the Lorentz-invariant density matrix for e + e~ — » QQ (here Q 2 = s): 

\M\ 2 = ^[(i +Pl )(£_ P2 ) + {l +P2 )(i_ Pl ) + m 2 Q Q 2 /2] , (34) 

where £-(£+) is the electron (positron) momentum and Pi(p2) is the momentum of 
Q(Q), and the quark mass is denoted by tjiq. In the second step, the heavy quark 
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(antiquark) fragments into a heavy hadron and a number of light hadrons, determined 
by a function fg(z), which peaks increasingly near z — > 1, as toq increases. In the 
third step, the heavy hadrons decay, dominantly non-leptonically, modelled on the 
quark transitions Q(p) — >• 91(91) + 92(52) + 93(93) [ZHZ]< These effects are important 
quantitatively for jet physics for the lower PETRA energies (typically < 30 GeV). 

5.2 3-jet events and cross sections at PETRA 

As discussed earlier, the characteristic feature of the process e + e~ — > qq with the 
subsequent fragmentation of the quarks and the antiquarks into a jet of hadrons is 
that it leads to a two-jet configuration. In QCD, the diagrams shown in Fig. mod- 
ify this picture. These corrections being proportional to a s (Q 2 ), the QCD coupling 
constant at the scale Q 2 , are small. However, the process e+e~ — > qqg may reflect 
itself in a structure of the final states that topologically is different from the dominant 
process e + e~ —¥ qq. The radiated gluon provides a new (non-local) mechanism for 
producing large-pT hadrons, which, unlike the px of the hadrons generated in the pro- 
cess e + e~ — > qq, is expected to increase with the e + e~ centre-of-mass energy. Thus, 
broadening of the transverse momentum of the hadrons with increasing centre-of-mass 
energy is a consequence of gluon bremsstrahlung. It was argued in |30] , that a corol- 
lary of this phenomenon is that a third jet should exist in the direction of the large px 
particle. In particular, if there is enough phase space available, i.e. for large enough 
Q, a three-jet topology in the shape of "Y" (Mercedz-Benz symbol) should emerge, 
clearly distinguishable from the (dominant) oblate cigar topology corresponding to 
two-jet events. 

The calculation for the process e + e~ — > q(pi) + 9(^2) +5(^3), shown in the upper 
two Feynman diagrams in Fig. |H] leads to the following (Dalitz) distribution: 

1 d 2 v = a s (Q 2 ) c x\+xl 
ao dx\dx2 27r (1 — x\){l — #2) 

where Q 2 = 4E 2 , Xi = Ei/E = 2Ei/Q, and Ei are the energies of the quark, anti- 
quark, and gluon, with x\ + X2 + £3 = 2, and <To is the lowest order e + e~ — > hadron 
cross section given in Eq. (|2Tj) . The differential cross section in Eq. ([33)) diverges near 
the end-points x\^ — > 1, and indeed has infra-red and collinear divergences. We shall 
discuss finite 2-jet and 3-jet cross sections in the next subsection, but for the present 
discussion these divergences can be removed by a reasonable cut-off procedure, such 
as a cut-off Qq on the invariant masses S13 = Q 2 (l — X2) and S23 = Q 2 (l — £1), 
yielding a finite lowest order three-jet fraction. 

5.3 Experimental evidence of three-jet events at PETRA 

While valuable tests of QCD were performed in studies of the T decays, based on the 
underlying mechanism T — > 3<? and the subsequent fragmentation of the gluons, three- 
jet events were first observed in e + e~ annihilation at PETRA in 1979 by the four 
experimental collaborations: TASSO [23], MARK- J [25], PLUTO [27] and JADE [25]. 
The process e + e~ — > qqg leads to planar events, the search of three-jet events in these 
experiments was concentrated mainly in demonstrating the excess of planar events 
compared to the estimates based on the 2-jet final states around ^/s = 27 GeV, where 
most of the early experiments at PETRA were carried out. Such quantitative analyses 
were backed up by topologically well separated 3-jet events. Fig.[T3]shows momentum- 
space representation of a representative two-jet and three-jet event measured by the 
TASSO collaboration, analysed on the basis of sphericity tensor and jettiness |40j . 
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Fig. 14. Momentum space representation of a two-jet event (a) - (c) and a three-jet event 
(d) - (f) in each of the three projections, (a),(d) fi2 - plane; (b),(e) h\ - fi2 plane; (c), (f) 
hi - fi3 plane. Here hi are the three axes of the sphericity tensor. (From TASSO [25]). 



The MARK-J measurement of the distributions in oblateness (denned below) at 
y/s = 17 GeV and at higher energies (27.4 + 30 + 31) GeV are shown in Fig. ITSTa) 
and Fig. [15] (b), respectively. For this measurement, the coordinate system is defined 
by the thrust axis ei, the major axis e^, which is in the plane perpendicular to ei, 
and is in the direction along which the projected energy in that plane is maximised, 
and the minor axis, e3, which is orthogonal to both ei and e2. Oblateness is then 
defined as 

C -^major -^minor ; (*^) 

where F major = £\ Pi- e 2/ J2i \Pi\ and ^minor = Z)iPi-ea/Z)i The two frames 
on the r.h.s. of this figures show the energy flow in the event plane defined by the 
thrust and major axes (upper frame) and by the thrust and the minor axes (lower 
frame). These measurements were compared with the qq (two-jet) and qqg (three-jet) 
Monte Carlo models [31 32 , and clearly favoured the qqg description, in a statistically 
significant way. 

PLUTO studied the averages of the momenta of the charged particles (p\\), where 
Py is the longitudinal momentum, (p±) and (|>j_)j measured relative to the thrust 
axis of the event as a function of the cm. energy. Their analysis showed that the 
quantities (py) and (pj_) are not very discriminative between the qq and qqg, but the 
energy dependence of (pj_) is better described if gluon bremsstrahlung is included. To 
study this effect in more detail, they distinguished for every event the two jets which 
are separated by a plane perpendicular to the thrust axis. The jet with the lower 
(higher) average (p±)is called the slim (fat) jet. Fig. [TBTa) shows (p\) of the charged 
particles as a function of the cm. energy, where the average is taken over the charged 
hadrons in all slim (fat) jets. For the slim jet the qq and qqg predictions from the 
Monte Carlo [3T] are very similar and the data are in agreement with both. For the 
fat jet, however, the data clearly favour qqg, and qq is ruled out. Fig. ITBT b) and [TBI c) 
show the so-called "sea-gull plot", obtained by plotting the variable x p = p/pbe&m 
and (pj_), at lower cm. energies 13 and 17 GeV and at higher energies 27.6, 30 and 
31.6 GeV, respectively. At the lower energy (Fig. ITBT b)). there is very little difference 
between qq and qqg predictions. For the higher energies(Fig. [TBT c)). qqg predicts a 
genuine one-sided jet broadening caused by the gluon jet; the effect is quite dramatic, 
especially at high x p . TASSO collaboration [35] has done a very similar analysis. 

The JADE analysis is based on the normalised sphericity tensor S a p (defined 
in Eq. ([TO]) ) and the resulting eigenvalues Qi,Q2,Q3 obtained by diagonalising this 



28 




90° Major 



( b) 



27.4 + 30 +31.6 GeV 
446 EVENTS 





Fig. 15. Left-hand frames: Normalised Oblateness distribution at y/s = 17 GeV (a), and at 
yfs = 27.4 - 31.6 GeV (b). The solid curves are the predictions based on a QQg model and 
the dashed curves are based on the QQ model with (pr) = 325 MeV (denoted as qqg and 
qq, respectively, in this review). The dashed-dotted curve in (b) is the QQ model prediction 
with (pt) ~ 425 MeV (Q = it, d, s, c, b). Right-hand frames: Energy flow in the event plane 
defined by (a) the thrust and the major axes, and (b) by the thrust and the minor axes 
with the events satisfying the cuts thrust < 0.8 and oblateness > 0.1 at yfs = 27.4 - 31.6 
GeV. The energy value is proportional to the radial distances.; dots are the experimental 
measurements (From MARK- J [26]). 



tensor on an event by event basis. The variables which play a central role in this 
analysis are the sphericity = 3/2(Qi + Q2) and planarity = {Q% — Qi). Fig. [T71 shows 
the planarity distribution dN/d(Q2 — Qi) measured by JADE at s/s = 27.7 and 30 
GeV. Their data are compared with a qq model, with cr q — 250 MeV and 350 MeV, 
both of which fail to describe the data. The qqg model describes the data well. 

The results reviewed in this section were the first measurements through which the 
effect of a third (gluon) jet was convincingly established in e + e~ annihilation. This 
is an important milestone in the confirmation of QCD in which jet physics played 
a central role. From a theoretical point of view, observation of the gluon jet was 
inevitable. Like many other discoveries in particle physics, this discovery needed high 
energy e + e~ beams, particle detectors well equipped to measure the characteristics 
of the hadrons, and data analysis techniques. This was the work of dedicated teams 
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Fig. 16. Average observed (p?) of charged particles in the slim and fat jets as a function of 

the cm. energy (a). Sea-gull plots {(p%) as a function of x p = p/pbcam, where p = yjp\ +Pr) 

for slim and fat jets in two separate energy range (b), (c). The solid and dashed curves are 
qqg and qq predictions, respectively. In (c), the dashed curve corresponds to a q — 0.247 GeV 
(default value) and the dash-dotted curve to a q = 0.35 GeV. (From PLUTO [27]). 



of machine builders and experimental physicists who should be credited with the 
discovery. For the interested readers we refer to individual accounts leading to the 
discovery of the gluon jets [112113211331134113511361137] . but stress that this list of 
references is by no means exhaustive. 



5.4 Quantitative studies of QCD at PETRA and PEP 

Subsequent to the discovery of the gluon jet, the four PETRA collaborations, JADE, 
Mark-J, PLUTO (later replaced by CELLO) and TASSO collaborations made many 
more measurements in e + e~ annihilation to hadrons, in which further evidence for 
the gluon jet was presented, These prompted quantitative studies of QCD for inclu- 
sive jet-observables, like thrust and the Fox- Wolfram shape variable |138j etc., and 
for jet topology, like the 2-jet and 3-jet rates etc. Also the gluon spin was deter- 
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Fig. 17. The planarity distribution compared with the model predictions for e + e" 
and e+e" -> qqg at yfs = 27.7, 30 GeV. (From JADE [2H]V 
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mined, following a suggestion in |139] , One important issue was the universality of 
the quark-gluon coupling a s (Q 2 ), i.e. to check whether the same value is obtained 
independent of the observables and the measurements of a s (Q 2 ) at various values of 
s = Q 2 were consistent with the evolution anticipated by the renormalisation group. 
These attempts to obtain a s (Q 2 ) required the calculation of next-to- leading order cor- 
rections to the topological jet-rates and inclusive jet-observables, and also required 
a better understanding of the non-perturbative models used to interpret the experi- 
mental data. These theoretical and phenomenological studies often took the form of 
detailed Monte Carlo programs without which no realistic comparison of theory and 
experiment was possible. In fact, since the days of experimentation at PETRA and 
PEP, Monte Carlo based theoretical frameworks have become indispensable for the 
quantitative analysis of data, as witnessed later at LEP, HERA and the Tevatron, 
and now at the LHC. 

In 0(a s (Q 2 ), 2-jet cross sections defined by a jet-resolution criterion, such as the 
Sterman- Weinberg jet-cones or the jet invariant mass, receive contributions from the 
virtual corrections to the process e + e~ — > qq, and soft or collinear configurations from 
the processes e + e~ —> qqg. In 0(a 2 (Q 2 )), the 3-jet cross sections receive contribution 
from the virtual corrections to e + e~ — > qqg and soft and collinear configurations from 
the 4-parton processes e + e _ — > qqgg and e + e~ — > qqqq. The hard and non-collinear 
configurations in the 4-parton processes give rise to 4-jet cross sections, with the 
leading contribution arising in 0(a 2 (Q 2 )), whose rates were calculated in [1401141) 
including the quark mass effects. They were important to check the non-abelian char- 
acter of QCD, as discussed later. The first complete next-to-leading order correction 
to event shapes up to order a 2 were undertaken by Ellis et al. [1421143] . They pre- 
sented their results in terms of the tensor 

^=E(pipi)/b-i)(Eb-ir 1 . ( 37 ) 

a a 
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where p l a are the components of the centre-of-mass three-momentum of hadrons a, 
and the sum runs over all hadrons. The eigenvalues of 9 are determined by the char- 
acteristic equation 

A 3 -A 2 + ic*A-L>/27 = 0, 0<C,D<1. (38) 

The quantities C (also called the Fox- Wolfram shape variable 138 ) and D are sym- 
metric functions of the eigenvalues, defined as 

C = 3(AiA 2 + A 2 A 3 + A 3 Ai), D = 27AiA 2 A 3 , (39) 

Integrating ~ in the range ^ < C < 1 yields 

i/'V^C^IW*!). (40) 
CT o Jo. 5 dC TT TT 

Numerically, C x = 2.8 and C 2 = 18.2 ± 0.7 for five quark flavours [1421143] . Thus, 
large corrections are obtained for the Fox- Wolfram shape variable, C . 

Another, and experimentally widely studied, example of an inclusive distribution 
is thrust. In 0(a 2 ), this was first calculated by Vermaseren et al. [144) . and verified 
subsequently by Ellis and Ross [145] and by others |146I147I148| . using the earlier 
work reported in [1421143] . In next-to- leading order NLO in a s , thrust-distribution 
in e + e~ — » hadrons is given by the following expression 

d ° - M T)^l + M T) { ^lf, (41) 



dT 



TT 



where the functions Ag(T) and A\(T) are shown in Fig. [18] (note that the variable t 
used in these plots is the same as T used in the text). The shapes of Aq(T) and Ai(T) 
are rather similar, but the 0(a 2 s (Q 2 ) corrections to the thrust-distributions are also 
numerically large. Integrating the distribution in Eq. (j4Tj) up to T = 0.85 yields 



if^^,*),,^,*,, (42) 

°0 Jo. 5 dl TT TT 

Numerically, K\ = 1.156, K 2 = 17.6 ± 0.3 for five quark flavours, which for a s (Q 2 ) = 
0.13 at -y/s = 35 GeV yields a correction of about 70%. Theoretical calculations 
from [142] were later implemented in the independent jet Monte Carlo 32 and used 
to determine a s (Q 2 ) from the inclusive measurements. The first such determination 
using the thrust and oblateness distributions, measured by the TASSO |149j and 
MARK-J [T5D] collaborations, respectively, yielded [HH] a s (Q = 35GeV) = 0.128 ± 
0.013 from the TASSO data and a s (Q = 35GW) = 0.120 ± 0.010 from the MARK-J 
data. 

Subsequently, an enormous effort has gone into estimating the effects from the 
jet resolutions, choice of jet variables, and fragmentation models. Also, the statistical 
significance of the data from the experiments at PETRA and PEP increased enor- 
mously over the time. An observable studied intensively in theory and experiments at 
PETRA is the energy-energy correlation (EEC) and its asymmetry (AEEC). EEC is 
a measure of the energy flow involving two calorimeters subtending solid angles J? and 
fl' with respect to the incoming e + e~ axis. Keeping the orientation between the two 
calorimeter cells fixed (= x), the differential distribution in cosx can be expressed as 

ldS EEC l^f da 



- f - — 7-^-, XiX-jdxidxj , (43) 



a dcosx & J dxidxjd cosx 
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Fig. 18. The scalar functions Ao(t) and A\ (t) for the Thrust distribution as defined in Eg. 1411 
Solid histograms (based on Ellis, Ross and Terrano |142| 1. dashed histograms (Vermaseren 
et al. [Hi]). (From [135] V 



where Xi are the scaled energies in terms of the cm. energy y/s. The experimental 
configurations with fixed angle between the calorimeters \ and the polar angle of 
one of the calorimeters 9 are calculable in perturbative QCD [1511152] , However, 
most experimental measurements were carried out for the averaged EEC, obtained 
by integrating over cos 6, for which perturbative QCD yields the following expression 
(for m q = 0) 

1 dS EEC a (O 2 ) 

- -^-F(0 , (44) 



where £ 



1— cos x 



(To G?COSX 7T 

and F(0 is given by [1511152] 



F{0 = 6^(l-g) [2(3 " 6C + 2?2) + ln(1 " + 3C(2 ~ 3?)] 



(45) 



The averaged (obtained by integrating over cos 9) AEEC cross section has an obvious 
definition 

1 dS AEEC _ 1 dS EEC (7T~ X ) 1 dS EEC ( X ) 



Co d cos x 



<jq d cos x 
a s (Q 2 ) 



[(F(l-0-F($] 



dcos x 
_ a s (Q 2 ) 



M0 



(46) 



Effects of quark masses in the EEC and AEEC cross sections were calculated in [T53 154 155 1156] , 
The 0{a 2 s (Q 2 ) corrections to these distributions were calculated numerically {153115411571158] . 
Restricting the angular range to —0.95 < cosx < 0.95, where the non-perturbative 
effects are relatively small, the NLO corrections to the EEC cross-section were found 
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Fig. 19. The asymmetric part of the energy energy correlation cross section measured 
by the JADE [159] and [TASSO] [161] collaborations at PETRA and comparison with the 
perturbative QCD calculations including 0{a 1 s {Q 2 )) corrections from [153] . The distribution 
in the upper left-hand frame from the JADE collaboration shows the corrected asymmetry 
A(9) vs. 9, measured at v / s = 14, 22 and 34 GeV. The upper right-hand frame from TASSO 
shows the measurements at ^/s — 43.5 GeV and comparison with the perturbative QCD 
(solid curve). 



to be moderate, typically 0(35%), but the corresponding corrections to the AEEC 
were small, typically 0(10%), giving reasons to be optimistic about the convergence 
of perturbative QCD in these variables, particularly the AEEC. 

Measurements of the EEC and AEEC were undertaken by all four experiments 
at PETRA: JADE, MARK-J, TASSO, and PLUTO. The AEEC measurements have 
been used to determine a s (Q 2 ) by comparing them with the NLO expression [153]. 
The extracted values of a s {Q 2 ) are found to be: a s (Q 2 = (34 GeV) 2 ) = 0.115 ± 
0.005 [JADE] [E3], a s (Q 2 = (34 GeV) 2 = 0.13 [MARK-J] [HO], a s {Q 2 = (34.8 GeV) 2 ) 
= 0.125±0.005 [TASSO] [MI], and a s {Q 2 = (34.6 GeV) 2 ) = 0.125±0.005 [PLUTO] [1521 . 
Within errors, these values of a s (Q 2 ) are consistent with each other, and with the 
ones from oblateness and thrust distributions, given earlier. Representative distribu- 
tions from the JADE [T^H] and [TASSO] [TEE] collaborations are shown in Fig. [U 
in which A{9) vs. 9 and 1 / adS A /d cos x vs. 1 — cosx are plotted, respectively. These 
measurements are compared with the perturbative QCD expression, calculated to 
0{a 2 {Q 2 )) and the agreement is impressive for (9,x) > 30°. For (9, x) < 30°, one 
needs to implement non-perturbative effects as well as the resummation of the large 
logs to all orders in perturbation theory. 

A lot of experimental effort went also in studying the topological cross sections (jet 
multiplicity) in e + e~ annihilation experiments at PETRA, PEP, TRISTAN and later 
at LEP. Theoretical distributions to these topologies were calculated in a series of pa- 
pers [93I163I164I165I166I167] . Making use of this theoretical work, the JADE collab- 
oration measured a s (Q 2 ) in a limited range of \fs using the three-jet rate and estab- 
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lished the running of a s (Q 2 ). Defining the fractional three-jet rate R3 = era-jet/ &tot 
as a function of y m i n , which is a cut-off parameter such that yij > y m in for any pair of 
partons i and j and y^ — M^/s, the measured jet-rate was fitted to the expression 

i? 3 (2/min) = C ias (Q 2 ) + C 2 a 2 s (Q 2 ) , (47) 

where C\ and C2 are y m i n -dependent constants calculated by Kramer and Lampe 
(called KL below) in |166| . The JADE measurements for Rs(y m in) as a function of 
y/s in the range 20 < y/s < 44 GcV are shown in Fig. (left-hand frame) [168] . 
They follow nicely the RG-prescribed running of a s (Q 2 ) with Aj^-g — 205 MeV for 
0.04 < j/ m ; n < 0.12 using KL, with almost the same value Ajj^ = 210 MeV using 
a calculation by Gottschalk and Shatz (called GS) [169 . An even more convincing 
measurement of the running of a s (Q 2 ) was presented by the MARK II collabora- 
tion [170] at PEP and SLC. They determined a s (Q 2 ) from the differential three-jet 
rate ga{yz)\ V3 = Vcnt = s^-Ms/cut), where / 2 (j/ C ut) is the fraction of two-jet events 
defined by the jet resolution y cu t- Their result for 33(2/3) is shown as a function of ys 
in Fig. HOI (right-hand frame) for two values yfs~ = 91 GeV (SLC) and y/s = 29 GeV 
(PEP). The three curves shown are the predictions of KL |166j for three different 
values Ajjg = 0.1, 0.3 and 0.5 GeV. Here Ajjg refers to the QCD scale parameter in 
a specific renormalisation scheme, the so-called modified minimal subtraction scheme 
M S |171j . For further reading of the technical issues of renormalisation and schem- 
dependencies at a nonspecialist level, we refer to a review on perturbative QCD [172] . 
These measurements yielded a 8 (Q 2 ) = 0.123 ± 0.009 ± 0.005 at Q = yfs = 91 GeV 
and a s {Q 2 ) = 0.149 ± 0.002 ± 0.007 at Q = y/s = 29 GeV, The running of a s {Q 2 ) 
is clearly established. A comparison with the values of a s (Q 2 ) determined from the 
measurements of the AEEC cross section at PETRA energies, discussed earlier, also 
shows that non-perturbative effects at these energies are observable dependent and 
not negligible. 

These investigations were extended to jet rates of higher multiplicity, i.e. four-jet 
and five-jet. An earlier paper along these lines is due to the JADE collaboration, in 
which n-jet rates (n — 2,3,4,5) were presented [94 . At this time, NLO corrections 
to the 4-jet rates and even LO predictions for the 5- jet final states did not exist. 
The data were compared with the leading-logarithmic- approximation (LLA) model 
. Similar studies based on the MARK II data at PEP (y/s = 29 GeV) are found 
in [173] and |174) using the so-called "optimised perturbation theory", i.e., by fitting 
the scale. 

An earlier attempt to establish the non-abelian nature of QCD from a study of 
multijet events was made by the AMY collaboration at the TRISTAN e + e~ storage 
ring at the KEK laboratory [1 75] . Their data showed a clear preference for QCD in 
contrast to an abelian model. In addition, they showed the running of a s (Q 2 ) by 
measuring the 3-jet rate R3 at y/s — 50 to 57 GeV by comparing their measurements 
with those of the JADE collaboration [168 and the TASSO collaboration [176] at 
PETRA taken at lower cm. energies. Other publications towards a determination of 
a s from PEP and PETRA are for example by MARK II [T77] and CELLO [178] . 

5.5 String- and String-like efFects in Jets 

The data taken by the experimental collaborations at PEP and PETRA have been 
used also to investigate non-perturbative effects in the jet profiles with the view of 
testing various phenomenological models available in the 1980's. This was important, 
since depending on the observables considered, non-perturbative effects influenced 
also the measurement of a s . Several groups [179118011811182] have used three-jet 
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Fig. 20. Left-hand frame: Three-jet event rates measured by the JADE collaboration 
as a function of the cm. energy E cln [GeV] for the indicated values of the jet resolution 
parameter y C ut, together with the predictions of the order a 2 s perturbatve calculations by 
Gottschalk and Shatz (GS) and Kramer and Lampe (KL). Right-hand frame: Experimental 
distribution (73(2/3) as a function of 2/3 at (a) y^s = 91 GeV and ^/s = 29 GeV measured by 
the MARKII collaboration [17U] . The 2/3 range used in the fit for the determination of a s is 
defined by the two dashed lines. The curves are second order perturbative calculations with 
the indicated values of A-rvw. 



(qqg) events to study the impact of hard gluon bremsstrahlung on the hadronisation 
process. In these studies they observed the so-called string effect 07], predicting a 
depletion of particles in the angular region between the quark and antiquark jet 
relative to the particle flow in the regions between the quark and gluon jets and the 
antiquark and gluon jets. In Fig. 1211 (left-hand frames), we show the measurements of 
the normalised energy flow (l/E)dE/d8 in planar three-jet events and the normalised 
charged particle flow in these events undertaken by the JADE collaboration [180] 
between yfs = 30 GeV and 36 GeV at PETRA. These distributions allowed one to 
distinguish between a hadronisation model |31j in which the fragmentation proceeds 
along the parton momenta (the independent jet IJ model) and the model in which 
the fragmentation takes place along the colour-anticolour axes (the LUND string 
model |38j). discussed earlier. Only the leading order (0(a s )) matrix elements were 
taken into account for the gluon bremsstrahlung process (e + e~ — > qqg), which were 
encoded in these fragmentation models. As seen in this figure, JADE data on the 
energy and charged particle flow are better reproduced by fragmentation along the 
colour axes [35]. 

A similar analysis was undertaken somewhat later in 1985 by the TASSO collab- 
oration |182j . In this case, the three-jet events produced in e + e~ annihilation into 
hadrons at 34.4 GeV were compared with the 0{a 2 s ) perturbative QCD calculations 
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Fig. 21. Left-hand frames: (a) The normalised energy flow l/EdE/d9 in the three-jet events 
compared with two model predictions, (b) The normalised charged particle flow 1/ndn/dO. 
(c) XjndnjdQ with py ut > 0.3 GeV. Here n is the total number of particles used in each plot 
(JADE collaboration [180] '). Right-hand frames: Ratios of particle densities in the angular 
gaps between the jet axes, defined by 0.25 < ip'j < 0.75 as a function of X[ n . The calculation 
of the IJ (g — q) and of the LUND models are shown as shaded bands, a) N(2)/N(3) and 
b)N(l)/N(3) for the three-jet event sample (TASSO collaboration [182] ). 



convoluted with two different models of fragmentation (IJ and Lund). The analysis 
was undertaken in terms of the "reduced azimuthal angles" tp'j and x- m = p- ul / -Ebcam, 
where p; n is the particle momentum projected into the event plane. The ip'j are defined 
as 

ib — 

ib'i = , i, j, k = 1,2,3 and cyclic , (48) 

J <P k - <Pi 

where the particle under consideration is located between jets i and k (<Pi < ib < ^&). 
The reduced angles ib'j run from and 1. The subscript j denotes the angular region 
opposite to the jet j. The analysis was restricted to X[ n < 0.1 and the data were 
divided in two samples x [n < 0.04 and 0.04 < x- m < 0.1. The result of the TASSO 
analysis is displayed in Fig. [5T] (right-hand frames) showing that the distribution of 
low energy (soft) hadrons in the 3-jet plane is better described by the LUND colour 
fragmentation model than by the independent jet model. The opposite is true for 
more energetic particles flowing between the 3 jets. 

The "string effect" was subsequently attributed to the coherence of soft gluon 
emission from the qqg system - a characteristic feature of the non-abelian nature of 
QCD |4B. • This is illustrated by contrasting the case of a soft gluon emission (assumed 
here as g{p2)) in e + e~ — > q{p+) + q(j>-) + g(pi) + g{p2) from the process in which the 
gluon g{p\) is replaced by a photon, i.e., e + e~ — > q{p+) + q{p-) +7(pi) + <?(P2)- The 
angular distribution of the soft gluon (antenna pattern)in the case of e + e~ — > qqj is 
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given by 

W + ^ 2 ) = 2C F a + -V(a,p) = - - l=f) , (49) 

cos a — cos p \ sin a sin p J 

where a = 02 and /3 = 0_| — 02 (see the kinematics shown in the upper left-hand 

frame in Fig. 1221) ; a%k = 1 — (n^n/j), with rij being the unit vector in the direction 

of and 6*_| is the angle between the q and q directions. Replacing ~f(pi) with a 

gluon g(pi) changes the angular distribution essentially due to the antenna element 
g(pi) participating in the emission as well. One now obtains (7 = 8+i + fo)'- 

W ±1 (<t> 2 ) = N c [a + iV(a,j)+ai-V(a, 7 )} + (2C F -N c )a + -V(a,p) . (50) 

The (soft) particle flow according to these two configurations is illustrated in Fig. [22] 
(upper right-side frame) showing that the flow opposite to the direction of ni is 
appreciably lower for the case of a gluon than for a photon due to the destructive 
interference in the case of QCD (qqgg). 

This phenomenon can be qualitatively understood . Omitting the small contri- 
bution from the second terms in Eq. (|50|) . one reduces this equation to the sum of 
two independent quark antennas (+, 1) and (— , 1). Therefore, in this approximation, 
the total particle flow can be obtained by the simple incoherent composition of two 
"annihilations" e + e~ — ¥ qq, boosted from their respective rest frames to the over- 
all qqg cm. frame. It is clear that the angular region between the q and q will be 
depopulated as it is opposite to the boost direction of both two-jet configurations. 
This perturbation theory based scenario (3 = 2 + 2 + Lorentz boost) then coincides 
with the fragmentation of the gluon in the process e + e~ — > qqg events in the LUND 
fragmentation model. The independent jet model misses this, as the gluon fragments 
independently on its own. Consequently, the Lorentz boost effect is absent. 

The colour coherence study of e + e~ jets by Azimov et al. [35] suggested an inter- 
esting experimental test in the form of particle flows in three-jet (qqg) and radiative 
two-jet (qqj) events by observing the negative contribution of the third antenna. This 
test was carried out by the MARK II collaboration at PEP at y/s = 29 GeV [183 
with the result that in the angular region between the quark and antiquark jets fewer 
charged tracks were observed in the two-jet events than in the radiative three-jet 
events. Their result is shown in Fig. I2"2l flower two frames). 

To end this review of the studies of jets at PETRA and PEP, we briefly discuss 
the angle ordered perturbation theory, as this approach has been used to develop 
a parton shower Monte Carlo [39]. In this approach, the phase space of soft gluon 
emission is restricted using an angle ordering criterion, which allows one to take 
into account the interference (colour coherence) approximately, and hence it repro- 
duces the string and string-like effects discussed above. Both the LUND fragmenta- 
tion model (PYTHIA in its modern incarnation) and the parton shower Monte Carlo 
models (such as HERWIG) describe the e + e~ data adequately. However, the main 
drawback of these models is that they do not (easily) match with the fixed order 
perturbation theory in next-to-leading and higher orders. The main obstacle is that 
fixed order perturbation theory has soft and collinear singularities that give rise to 
logarithmic enhancement of higher order contributions. These enhanced terms should 
be summed to all orders. However, there is no unique way of doing this. For example, 
the pT-ordered and the angular ordered showers can both be arranged to resum these 
logarithms. Matching with a fixed order perturbation theory is more easily achieved 
in pt ordered showers which, however, do not have the colour coherence needed by 
the low-energy e + e~ data. It is the other way around with the angle-ordered showers. 
We will discuss these aspects further in the next section. 
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Fig. 22. Upper left-hand frame: Kinematics of non-jet radiation in three-jet events; Upper 
right-hand frame: Directivity diagram of the soft particle flows, projected on to the qqy 
(dashed lines) and qqg (solid line) event planes. Dotted circles show the constant levels of 
density flow [W(<^>) = 1,2,4] (from Azimov et al. [48] ). Lower frame: The charged-track 
density as a function of the event-plane angle <f>. The angular region between <f> = 0° and 
4> — 150° separates the q and q for the qq-y events and for 65% of the qqg events (from the 
MARK II collaboration [183]). 



6 Jets in QCD and at LEP 



In this section we review the salient features of jets at LEP which were helpful in 
testing some of the basic elements of QCD more precisely. Just as in the preceding 
section, we recall the main detectors at LEP, ALEPH, DELPHI, OPAL, and L3, 
which collected large data samples (typically, 4 million hadronic events around the Z 
resonance for each of the four LEP experiments. In the second stage, LEP2, the beam 
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Fig. 23. The three basic splitting processes of quarks and gluons into pairs of quarks and 
gluons. 



energy was increased to about 103 GeV. These detectors proved to be very powerful 
tools in carrying out precision electroweak and QCD physics. 



6.1 Quark/Gluon cascades 

Electric charges which are accelerated, reduce their energy by radiating photons pref- 
erentially collinear with the flight direction of the charge. This is a general feature 
of gauge theories and, specifically, collinear radiation is predicted in QCD processes 
in which quarks and gluons are produced with high energies. If the observed partons 
carry away a fraction z of the parent partons, the splitting functions [184] . cf. Fig. [231 
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universally predict collinear splittings, with Q 2 ~ z(\ — z)E 2 2 denoting the invariant 
mass of the final parton pair. The notation defines a distribution such that 

for any sufficiently regular function f(z), 



dzf(z)[F(z)l 



dz(f(z)-f(l))F(z) 



(52) 



The bremsstrahlung splittings q qg and g — > gg preferentially generate soft radia- 
tion spectra in the limit z — > 0. The group characteristics are Cf — 4/3, Ca — 3 and 
T R = 1/2 for SU(3) C of QCD. 

Repeated splittings generate cascades of collimated quarks and gluons. Since the 
lifetime of the final-state pair in the splitting processes is long, t* ~ E/Q 2 , the cascade 
is expected to be described by a sequence of probabilities and not by interfering 
quantum-mechanical amplitudes. If the branching occurs at a value Q 2 without any 
radiation between the initial maximum value Ql^ax and Q 2 , the probability is given 

by 
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where the exponential Sudakov factor |71] accounts for the non-radiation probability. 
Here P a -^Vcf (z') are the same functions as P a ^b'c'{z'), defined in Eq. (|51"j) except for 
the regularization terms at z' — > 1. 

However, the branching probability Eq. (|53[) is refined by an important coherence 
effect. If in electrodynamics a photon splits into an electron— positron pair, the pair 
can emit photons only at angles less than the angle between the charged pair as 
photons propagating at larger angles would see coherent electron+positron states 
which, being neutral, cannot radiate. As a result, the emission angles are ordered 
in the sequence 0% > 02 > ■■■ This effect is also predicted in QCD [33]. The only 
difference arises from the fact that the coherent superposition adds up the colour 
charges of the daughter partons to the non-zero value of the parent colour charge so 
that wide-angle splitting is generated at a non-zero rate. In addition to the angular 
ordering, non-resolved infrared radiation restricts the energy fractions of the partons 
in the cascades. These restrictions on energies and angles can be mapped into the 
boundary values of the Sudakov integral after re-expressing the invariant mass by the 
angle between the momenta of the daughter partons. 

The cascading of the primordial quarks and gluons affects the observed hadron 
distributions within the jets. In particular, energy spectra are softened through the 
cascading mechanism, multiplicities increase strongly with energy, and quark and 
gluon jets will develop different profiles. Formulated for simplicity by neglecting the 
change in kr, and restricting to one parton species, the energy dependence of the 
fragmentation function is described by the DGLAP equation [184 185 1 861187] : 

dD(z,Q 2 ) _a s (Q>) r l % Pm (> l Q A. (54) 



d\ogQ 2 /A 2 2tt J z C V ' VC 

The splitting function P(Q consists of two parts (see, Eq. (j5l"j) . The first part P 
describes the standard component and accounts for the accumulation of particles at 
z generated by the splitting of partons at £ > z, the second part accounts for the 
loss of particles at z due to splitting to smaller energy values. In the parton model, 
the function D depends only on the variable | but not on Q 2 . This would yield an 
scale-invariant fragmentation function D(^). In QCD, this is obviously modified. The 
solution of the above equation leads to striking effects which modify the predictions 
of the scale-invariant parton model. In particular: 

(i) For large z values beyond 0.2 the spectrum decreases with increasing energy 
while the particles accumulate at small z. The loss of particles by splitting at large z 
is bigger than the gain by splitting from yet higher £ values. This is naturally opposite 
at small z values. 

(ii) The constant plateau (characterized by D(z) as in Eq. (5)) in the parton model 
generates a multiplicity of particles which increases logarithmically with the length 
of the plateau ~ \og(^/s/mi l ). Multiple splittings raise the multiplicity much more 
strongly. Solving Eq. (|54[) for the multiplicity, given by the integrated fragmentation 
function, predicts a rise with energy stronger than exponential: 

n(s) ~ exp(a; 1 / 2 (s)) ~ exp(log 1/2 s/yl 2 ) . (55) 

In addition, the flat plateau in logz -1 is deformed to a humpback of Gaussian char- 
acter (188j . with centre [logz _1 ] rooa: ~ logs/ A 2 and width a ~ log 3 ^ 4 s/yl 2 . Experi- 
mental proof for the energy dependence of the fragmentation function D(z, Q 2 ) and 
the formation of the humpback at small z is presented in FigJM] 

(Hi) While the first splitting of a quark jet q — > qg is determined by Cf = 4/3, 
the first splitting of a gluon jet g — > gg is determined by the bigger Casimir invariant 
Ca = 3. Thus, the large average colour of gluons compared with quarks should 
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Fig. 24. (a) Scaling violations in quark and gluon jets (using the DURHAM algorithm) . The 
solid curves show the DGLAP evolution and the dotted parts of the curves are extrapolation 
outside the fit range [DELPHI [189] ]. (b) Measurements of the ln(l/a;p) distributions for the 
center-of-mass energies between 14 and 91 GeV and comparison with the MLLA 188 and a 
Gaussian 190 . This figure is referred as the humpback plateau of the gluon fragmentation 
function for small x values in the text [OPAL [191] ]. 

increase the multiplicity of gluon jets with regard to quark jets asymptotically in 
the ratio Ca/Cf — 9/4. Similarly, since dN g / q ~ Ca/f dlogO, the angular widths 

of quark and gluon jets are different, O g — Oq F , viz. gluon jets are wider than 
quark jets. Even though the asymptotic limit has not been reached yet, the particle 
multiplicities in gluon jets have been shown significantly larger than in quark jets 
|192j : n g /n q > 3/2. For a recent discussion of particle multiplicities, we refer to [193] . 

(iv) Small-angle gluon radiation off heavy quarks Q — c, b is suppressed com- 
pared to light quarks 194 , and the logarithmic enhancement of the particle yield is 
restricted to infrared gluon configurations. 



6.2 Multijets at LEP 

Increasing the energy from the PETRA regime of about = 46 GcV to the LEP 
regime by factors of two and five in the two phases of LEP, Z-boson runs with -y/s = 
91 GeV and beyond with y/s up to 206 GeV, provided two opportunities: the exper- 
imental analysis of multijet final states [M] and the study of the jet properties over 
a large range in energy |195) . This allows a more precise measurement of two funda- 
mental characteristics of QCD [196], the running of the QCD coupling with energy 
as predicted by asymptotic freedom, and the three-gluon coupling, a fundamental 
ingredient of asymptotic freedom. We discuss these measurements below in turn. 



6.3 Inclusive jet observables and determination of a s (Mz) at LEP 

All four experiments at LEP, DELPHI Q£7], OPAL [US], L3 [US], and ALEPH [2TJD] , 
undertook measurements of the inclusive jet (or event shape) variables and their 
moments. In these analyses, the next-to-leading order (O(a^)) perturbative QCD 
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calculations for the inclusive observables and event shape distributions, discussed in 
the context of three-jet events at PETRA, were augmented by theoretical estimates 
in the next-to-leading-log approximation (NLLA) [2011202] and others in which the 
0(a 2 s ) and NLLA schemes were combined [203] . It is helpful to explain the NLLA in 
more detail. For a generic shape variable, y, well-defined in pertrubation theory (i.e., 
infrared and collinear safe), the typical leading behaviour at small y is 

I^L^ a »I m 2«-iI . (56) 
a dy y y 

The normalized event shape cross section R(y) defined as 

f v 1 da 

R{y)= / dy-^, (57) 

then has the expansion 

i2 Tl (y)~a?ln a "(l/y) = a,L 2n . (58) 

Whenever L is large, one can improve the range and accuracy of perturbative pre- 
dictions by identifying these logarithmically-enhanced terms and resum them to all 
orders. For a class of variables for which the leading logarithmic contributions expo- 
nentiate, i.e., variables which in the small-?/ range yield the logarithm of the shape 
cross-section in the form \xiR(y) ~ Lg\(a s L), where g%(a s L) has a power series ex- 
pansion in a s L, one has [203] 

R(y) = C(a s )E(y, a.) + D(y, a s ) , (59) 

with C{a s ) = l+Y,n=i C„(a s /27r)™ andlnZ'(y, a s ) = Lgi(a s L)+g 2 (a s L)+a s g 3 (a s L)- 
.... The function gi(a s L) resums all the leading contributions of the form a s L n+1 
(defining the LAA approximation) and g2(a s L) contains the next-to-leading loga- 
rithmic corrections a"L™ (defining the NLLA) and so on. 

For the hadronisation effects, either the Monte Carlo based hadronisation models, 
PYTHIA [35150] , HERWIG [5TJ and ARIADNE [201] were used, o r, alternatively, 
non-perturbative power-correction formulae derived in p05 206 207] were employed. 
This latter ansatz provides an additive term to the perturbative QCD estimate in 
mean event shape variables. Studying the energy-dependence of these mean variables 
(/), defined as 

(/) = —/#= (/pert) + (/pow(ao)) , (60) 

0tot J da 

yielded a measurement of a s (ji) and ao, the non-perturbative parameter characteris- 
ing the power corrections. We have discussed the 0{a 2 s ) calculations of / pcrt for several 
observables (thrust, the Fox- Wolfram shape variable C, etc.) in section 5.4. Explicit 
formulae for / p0 w are given by Dokshitzer et al. [205 206 207 , and can also be seen, 
for example, in the DELPHI analysis [197] of the LEP2 data. For related discussion 
of event shapes in e + e~ annihilation and deep inelastic scattering, we refer to |208j . A 
typical measurement along these lines is shown in Fig. 1251 in which the measured mean 
values of (1 — T), and the scaled heavy jet mass (M?/E% is ) are shown as a function of 
the center-of-mass energy together with the results of the fits. The dotted curves in 
these figures show the perturbative QCD part only. It is obvious from this figure that 
even at the highest LEP2 energy, non-perturbative power corrections are not small. 
The fits yield a s {M z ) = 0.1191 ± 0.0015 ± 0.0051 from (1 - T) and a very consis- 
tent value from the other observable (M%/E 2 is ). However, the value of ao(2 GeV), 
the measure of power corrections, differs by about 20% from the two measurements, 
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Fig. 25. Measured mean values of the observables (1 — T) and (M^/E^ ls ) as a function 
of the center-of-mass energy. The solid lines present the results of the fits including power 
corrections and the dotted lines show the perturbative part only. (From DELPHI |197j ). 



showing considerable non-universality in the parametrisation of / pow by Dokshitzer 
et al. [2U 5 206 207 . Along the same lines, the L3 collaboration measured the mean 
values of several global shape parameters. For all these variables, the same theoretical 
frameworks l 201 202 203 as discussed above in the context of the DELPHI measure- 
ments were used. To compare these calculations at parton level with the experimental 
measurements, the effects of hadronisation and decays were corrected for using the 
JETSET PS (parton shower) Monte Carlo program. We display in Fig. [26] the mea- 
sured distributions in thrust and the variable C at yfs = 206.2 GeV and comparison 
with the QCD fits, showing excellent agreement. To determine a s at each energy, the 
formalism in gH] is used, which yielded a s (M z ) = 0.1227 ± 0.0012 ± 0.0058 [155] . 



6.3.1 Jet rates 

Due to the high energy at LEP1, up to four jets could be resolved experimentally. 
The number of resolved jets depends strongly on the criterion by which the jets are 
defined. Early definitions had used the JADE recombination scheme which combined 
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Fig. 26. Measured distributions of thrust, T, (left-hand frame) and the C-parameter in 
comparison with QCD predictions at \fs =206.2 GeV [From L3 |199| ]. 

particle pairs on the experimental side, and equally quark/gluon parton pairs on the 
theoretical side [166] . for scaled invariant masses Mfj/s = y below a cut-off value 
< y C ut- In the DURHAM scheme 78l the invariant mass was replaced by Mfj — 
2 min(£f, Ej) [1 — cos(9y], essentially the transverse momentum between the particles 
or partons for small angles. The cut-off value y cu t was chosen typically from 10 _1 down 
to 1CT 3 . Small values of y cut naturally lead to large numbers of jets while the number 
of jets is reduced if y cu t is increased. The cross section for 3-jet events, 



has been calculated up to third order in the QCD coupling 209]. NLO corrections 
to the four-jet rates in the process e + e~ — » 7*, Z — > 4 jets were done around 1996 
by Dixon and Signer [210121 lj and subsequently by Nagy and Trocsanyi |212j . The 
experimental number of jets in Z decays is displayed in Fig. 1271 and compared with a 
parton shower Monte Carlo prediction (Jetset partons), and including hadronisation 
effects (Jetset hadrons). Evidently, for y cut below 10 -2 up to four jets can clearly 
be identified at LEP1 |213j . A dedicated effort to test QCD and determine a s was 
undertaken by the combined JADE (at PETRA) and OPAL (at LEP) collabora- 
tions [198] . giving a considerably larger lever arm in energy from 35 GeV to 189 
GeV. The observables used in this (JADE + OPAL) analysis are exclusively based 
on the multiplicities of hadronic jets. The n-jet fractions, R n were defined using the 
JADE [44], Durham [78] and Aachen/Cambridge[8"Tl algorithms. We show in Fig. [28] 
the three-jet fraction R3 obtained with the JADE and Durham jet algorithms versus 
the center-of-mass energy. (The result from the Aachen/Cambridge algorithm can be 
seen in [198] .) Here, the data from PETRA and LEP are compared with the 0(a 2 s ) 
prediction. The renormalisation scale dependence is shown by the scale parameter 
Xfj, = fji/s/s, with the solid lines corresponding to a fixed value = 1, and the 
dashed lines are the results obtained with a fitted scale, indicated on the figure. This 



and related analyses reported in |198] yield a rather precise value for the QCD cou- 
pling constant a s {M z ) = 0.1187° S°o 3 4 19 . At LEP2 (up to ^/s = 206 GeV), the highest 
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Fig. 27. Relative production rates of n-jet events defined in the Durham jet algorithm 
scheme [78] as a function of the jet resolution parameter y cu t- The data are compared 
to model calculations before and after the hadronisation process as indicated on the fig- 
ure [OPAL [215]]. 
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Fig. 28. Measured 3-jet fraction as obtained with the JADE (left-hand frame) and Durham 
schemes (right-hand frame) at parton level versus the c.m.s. energy *Js. The data shown are 
from the JADE and OPAL collaboration, and the curves are the 0{a 2 s ) predictions at a fixed 
scale (solid lines) and for the fitted values of the scale (dashed lines). [From OPAL [198] .] 



jet multiplicity measured is five, obtained using the variable j/ C utj and inclusive mea- 
surements are available for up to six jets. To match this data, NLO QCD corrections 
to five-jet production at LEP have been carried out by Frederix el al. [214] . and 
the fixed-order perturbative results have been compared with the LEP1 data from 
ALEPH [200 . Two observables have been used for this comparison: 
(i) Differential distribution with respect to the five-jet resolution parameter 2/45, the 
maximum value of y cu t such that a given event is classified as a five-jet event by the 
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DURHAM jet algorithm [75]: 

J dy45 J^ = o£d rt (»cut) , (62) 

where crf n ~J et is the inclusive five-jet production cross section in e + e~ annihilation, 
(ii) Five-jet rate i?5(y C ut), defined as follows: 

R 5 (y cut ) = ^xcT^cut) (63) 
Ctot 

where cr^j Ct (?/ cu t) is the exclusive five-jet production cross section. This is also cal- 
culated using the Durham jet algorithm by requiring that exactly five jets are recon- 
structed. Both observables, a^ \ da jd In and i?5(y C ut), can be written as a series in 
a s (/i), with the leading contributions starting in O(a^). A comparison of the leading 
order and next-to leading order predictions for (1 / 'a)do / 'din 2/45 vs. 0(3/45) and the 
exclusive 5-jet fraction R$ vs. ln(j/ cut ) is shown in Fig. [25] with the ALEPH data in 
a limited range of these variables (perturbative regime). Hadronisation effects have 
been estimated using the SHERPA MC [52 . As is typical of NLO calculations, scale 
dependence is significantly reduced compared to the LO calculations. Agreement be- 
tween data and NLO theory is impressive and has been used to extract a value of a s , 
obtaining a s (Mz) = 0.1156to'oo34i which is in excellent agreement with the world- 
average discussed below. The limitation of using fixed-order perturbative QCD in 
describing the e + e~ data can be seen in the ALEPH data shown in Fig. [201 which 
show a characteristic turnover shape around — lnj/45 ~ 7.5. In this region, pertur- 
bation theory fails and a resummation (equivalently showers) have to be included to 
describe the data. This underscores the importance of having MC generators which 
include showers. 

In Fig. [301 the ALEPH LEPI data for c t ~ ] da jd In are compared with the 
hadron level predictions of three event generators (the numbers denote the various 
versions of these Monte Carlo programmes), PYTHIA6.1 [US], HERWIG6.I [2To] 
and ARIADNE4.1 |204j . Agreement between data [200J and these MCs is generally 
good. However, as shown in the upper frame of fig. 1301 hadronic corrections are large, 
varying from 0.5 to 1.5 in this range. In addition, differences between hadronisation 
corrections among the MCs are as large as 25%. This deficiency can be overcome to 
some extent by matching the parton shower and high multiplicity matrix elements, 
as, for example, proposed in [217] . This matching procedure has been implemented 
in the SHERPA event generator [52] and results in improved agreement between the 
MC and fixed-order perturbative description. 



6.3.2 The gluon self-coupling 

The study of the three-gluon coupling in gluon splitting to two gluons requires four 
(or more) jets in e + e~ annihilation. A variety of angular correlations and energy 
distributions, see [218] - (22i] . can be exploited to signal the three-gluon coupling of 
QCD. 

The sensitivity to angular distributions may be illustrated in a transparent ex- 
ample [218] , A virtual gluon, radiated off the quarks in the process e+e~ — > qqg* , is 
polarised preferentially in the production plane. The subsequent splitting of the vir- 
tual gluon into two real gluons or a quark-antiquark pair is sensitive to the azimuthal 
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Fig. 29. ALEPH LEP1 data 200 compared to leading and next-to-leading order predictions 
for 1/crdcr/dln J/45 plotted against In j/45 (left-hand frame) and R — 5 plotted against ln(y cut ) 
(right-hand frame) without hadronisation corrections. The uncertainty bands are obtained 
by varying the scale in the interval 0.15Mz < /i < 0.6Mz, and the solid lines refer to NLO 
QCD evaluated at a B (M z ) = 0.118 and fi = 0.3M Z . [From Ref. [213] ]. 
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Fig. 30. ALEPH data 200 for the j/45 distribution at LEP1 compared to 
PYTHIA,HERWIG and ARIADNE Monte Carlo results. The upper frames show detector 
and hadronisation corrections, respectively. The lowest frame shows the relative difference 
between data and event generator predictions. [Figure attributed to H. Stenzel in Ref. [214| ], 
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angle <fi between the g* polarisation vector and the decay planes: 

n gg = _ h z(l - z) cos 20 , 

z(l z) 

n qq = i [z 2 + (1 - z) 2 ] - z(l - z) cos 20 . (64) 

As a result, the polarisation vector and the decay plane tend to be aligned in gluon 
splitting to two gluons. In contrast, if the gluon splits to a quark-antiquark pair, the 
decay plane tends to orient itself perpendicular to the polarisation vector. 

The azimuthal distribution can be studied experimentally by measuring the angle 
between the planes spanned by the two hardest and the two softest jets. In an abelian 
theory the <fi asymmetry is large and the two planes would orient themselves perpen- 
dicular to each other. By contrast, since the ^-independent term in gluon splitting 
to two gluons in QCD is large, the azimuthal asymmetry in this process is predicted 
to be small, but the two planes should have a tendency to orient themselves parallel 
rather than perpendicular. This is borne out by experimental analyses |222) indeed, 
as demonstrated in Fig. [3T] 

Quite generally, four jets are produced in e + e~ annihilation by three mechanisms: 
double gluon bremsstrahlung, gluon splitting to two gluons, and gluon splitting to a 
quark-antiquark pair. The cross section [1401141] can be decomposed accordingly: 

(ot \ ^ 
—J C F [C F a bb + C A o gg + n f T R a qq ] . (65) 

The first term accounts for double gluon bremsstrahlung q — > qg and q — > qg, the 
second for gluon splitting to two gluons g — ¥ gg, the third for gluon splitting to 
rif quark pairs g — > qq. The Casimir group characteristics of the splitting vertices 
are [Cf,Ca,Tr] — [4/3,3,1/2] in QCD, while the corresponding characteristics are 
[1,0,3] in an abelian theory. Measurements of their ratios yield [223] 

Ca/Cf = 2.29 ± 0m[stat.} ± O.U[syst] , 

T r /Cf = 0.38 ± 0.03[stat.] ± Q.06[syst.] , (66) 

compared with the theoretical predictions Ca/Cf = 9/4 and Tr/Cf = 3/8 in QCD. 
Again we observe a strong signal of the three-gluon coupling in Ca, far away from 
zero in the abelian theory. 



6.3.3 QCD coupling and asymptotic freedom 

A large range of energies can be covered in the measurement of the QCD coupling 
a s (Q 2 ), extending from 29 GeV (at PEP) to 46 GeV (at PETRA) up to about 206 
GeV at LEP; the lever arm can be elongated down to 1.8 GeV by including r decays. 

The QCD prediction for the running coupling a s (Q 2 ) has been determined up 
to the 5th power 224 225 . Keeping terms up to 2nd order leads to the following 
expression 

(n2 s i _ Aiogiog(Q 2 M 2 ) 

(3 Q \og{Q 2 /A 2 ) /3^og 2 (Q 2 /A 2 ) + "" l&0 

with f3 a = (33 - 2n / )/12 7 r, ft = (153 - 19n / )/24 7 r 2 , ... and A « 200 MeV denoting 
the QCD scale at which the coupling grows indefinitely. An ensemble of observables 
has been calculated in perturbative expansions in next-to-leading order NLO up to 
N 3 LO. Most accurate are the totally inclusive observables, like total cross sections, 
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Fig. 31. The distribution of the azimuthal angle xbz (called <j> in Eq. (59)) between the 
planes formed by the two most energetic jets and the two least energetic jets 222 ; the 
predictions of QCD including the gluon self-coupling are compared with an abelian theory 
without self-coupling of the gauge fields. 



followed by jet cross sections and hadronic shape variables, like thrust. The estimates 
still depend significantly on the models used for calculating the shape variables in the 
non-perturbative region, see [2261227] . for instance. 

Combining the experimental measurements with the theoretical apparatus, the 
knowledge of the QCD coupling and its evolution with the energy is summarised in 
Fig. [35] The lever arm extends from the hadronic decays of the r lepton through- 
out the PETRA range up to the highest energy values in the second phase of LEP. 
Including deep-inelastic lepton-nucleon scattering and jet production in hadron colli- 
sions, all the analyses are in remarkable agreement with the theoretical expectation 
from asymptotic freedom [3J. It has become customary to quote the value of a s (n) 
measured in experiments by scaling the result to the scale fj, = Mz using the RG 
equation. This yields the current world average |228] 

a s (M z ) = 0.1184 ± 0.0007 . (68) 



This ends our discussion of jets in e + e~ annihilation experiments and in QCD. In 
summary, essential parts of QCD jets can now be controlled at the level of typically ten 
percent (a s (Mz) is known better than 1%). Notabene the basic interactions and the 
strength of the quark-gluon coupling are proven to be asymptotically free. The high 
level of accuracy achieved in measuring the gauge couplings - weak, electromagnetic 
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Fig. 32. Left-hand frame: Summary of measurements of a s (Mz), used as input for the 
world average value. Right-hand frame: Evolution of the QCD coupling with energy. Among 
other methods, analyses based on hadron production in e + e~ annihilation play a leading 
role up to the highest LEP energies. Both frames are from Bethke [228| . 

and QCD - is now a diagnostic tool to probe physics at scales as high as the the grand 
unification scale. 



7 Jets as tools 

In the preceding sections, we discussed the impact which jet physics in e + e~ an- 
nihilation experiments had in establishing QCD quantitatively. This progress owes 
itself to some extent to the fact that in e + e~ annihilation the initial state is pre- 
cisely known. This is not the case in other high energy collisions, such as the electro- 
and photoproduction processes ep and jp, as well as the gamma-gamma and hadron 
hadron collisions, 77, pp and pp. Here, jets could be used as powerful tools for study- 
ing other aspects of high energy collisions. Examples are the partonic composition of 
the proton, i.e., quark and gluon densities of the proton (and antiproton), the par- 
ton distribution functions (PDFs) of the photon and the QCD coupling at HERA, 
Tevatron and the LHC. Yet other applications of jet physics include analyses of the 
electroweak sector and searches for new heavy particles in many extensions of the 
Standard Model (SM) - the QCD and electroweak theory of particle physics based 
on the groups SU(3) C <g> SU (2)/ <g) U(l)y Thus the prominent role of jets in studying 
QCD phenomena extends to quite different areas in particle physics. 

Before we embark upon illustrating the use of jets as tools in ep, 77, pp and pp 
collisions, it is worth pointing out that in these processes, QCD is at work in both 
the initial and the final states as opposed to the e + e~ annihilation processes, where 
it influences only the final state distributions and rates. As seen in Fig. [33] for the DIS 
process, the cross section depends on three components: (i) the probability of finding 
a parton in the proton having a fractional longitudinal momentum x (or £Bj), (ii) the 
interaction between these partons and the virtual photon, and (iii) the transition of 
partons to jets, which theoretically involved the recombination of two partons into 
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one jet. While perturbative QCD provides a framework to evolve the PDFs and the 
fragmentation functions FFs from a low scale y 2 = Qq to a high scale /i 2 = Q 2 , 
non-perturbative inputs for the PDFs and FFs are required at the lower scale. This is 
obtained by parametrising an ansatz at lower scale. The theoretical tool which enables 
this is called factorization, a key concept in the application of QCD to high energy 
processes [229j . Simply stated, factorization of a process (such as inclusive- or jet- 
cross sections in deep inelastic scatterings) allows it to be expressed as the product of 
a short-distance part, calculable in perturbative QCD, and a long distance part, com- 
prising of non-perturbative matrix elements or PDFs. The universality of the PDFs 
and the evolution from a lower scale to a higher scale are process-independent. The 
division into the short- and long-distance parts is governed by the factorization scale 
Hf- We illustrate the applications of factorization on the example of deep inelastic 
scattering processes, discussed below. 



7.1 ep Collisions 

In DIS, described here by the process ep — > e + X, an electron e with momentum k 
emits an off-shell photon with momentum q which interacts with a proton of momen- 
tum P. For virtualities of the photon (Q 2 = —q 2 > 0) far above the squared proton 
mass (but far below the Z-boson mass), the differential cross section in terms of the 
kinematic variables Q 2 , x = Q 2 /(2P.q) and y = (q.P)/(k.P) is 

^-^[(i+d -»)■)«(«.«•) -^CfJ-)] ■ m 

where F%(x, Q 2 ) and Fl(x, Q 2 ) are proton structure functions, which encode the inter- 
action between the photon and the proton. The structure functions are not calculable 
in perturbative QCD. In the lowest order, i.e., keeping only the Born contribution as 
shown in Fig. [33] (a) where x = XBj , the structure functions are given by 

F 2 (x,Q 2 )=xJ2c 2 q f q/p (x) : F L (x,Q 2 )=0, (70) 
i 

where f q / p {x) is the PDF for quarks of type q inside the proton. This result in which 
fq/p (x) are independent of the scale Q is the "quark-parton model" picture. Hence, in 
this picture, also the structure functions F% and Fl are independent of Q 2 . Including 
higher order perturbative QCD corrections, the structure function F%(x, Q ) has the 
form [230] 

n=0 ^ 17 ' Jx Z Z 

Here we have a series in powers of a s (/i 2 ^ ) , each term involving a coefficient C^} , which 
can be calculated using Feynman graphs. The scale /ir is called the renormalization 
scale at which the QCD coupling constant a s (n 2 R ) is calculated. The leading order in 
a s QCD Compton scattering and the boson-gluon fusion contributions are shown in 
Fig. [231(b) and[351(c). respectively. An important point to note is that the momentum 
of the quark when it interacts with the photon and the momentum when the quark is 
extracted from the proton may differ. The ratio of these two momenta is z. The 
coefficients depend on the ratio z, and hence one has to integrate over z as indicated 
above. In lowest order, i.e. without including any perturbative QCD corrections, one 
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Fig. 33. Deep-inelastic lepton-proton scattering at different orders in a s : (a) Born contri- 
bution O(l), (b) example of the QCD Compton scattering 0(a 3 ) and (c) boson-gluon fusion 
O(o.). 



has C^l = SqS(l — z) and C^l — 0, and one recovers the "quark-parton model" 
result. 

The PDFs fi/ p {-,lip) depend on the factorisation scale (x 2 Fl and this dependence 
is governed by the DGLAP equation. In leading order in a s ([i F ), this reads as follows 



3 x 



dfi% 



where, for example, Pq}\ g (z) = Tr(z 2 + (I — z) 2 ), and the others can be extracted from 
Eq. (1511) . The coefficient functions are also //^-dependent, and in the leading order in 

a„ one has C^x, Q\ n% /4) = C 2ti (x, Q 2 , Q 2 ) - hi(§)E j St 
In the above expressions, the choice of the renormalization and factorization scales is 
arbitrary. Varying \i p and [Ir provides an estimate of the scale-dependent uncertain- 
ties. In inclusive DIS processes, the default choice is fijt = = Q- 

The extension of the factorization formalism to the processes with two initial-state 
hadrons follows very much along the same lines, and we shall discuss this somewhat 
later, as we discuss high energy hadronic collisions. 



7.1.1 Jets in DIS processes 

Jet production in neutral current (NC) deep inelastic scattering at the HERA collider 
at DESY provides an important further testing ground for QCD. While inclusive DIS 
gives indirect information on the strong coupling via scaling violations of the proton 
structure functions, the production of jets allows one a direct measurement of the 
strong coupling constant a s . The Born contribution to DIS (Fig. 155b ) generates no 
transverse momentum in the "f*p centre-of-mass frame, where the virtual boson and 
the proton collide head on. Significant transverse momentum in the 7*p frame is 
produced at leading order (LO) in the strong coupling a s by the QCD-Compton 
(Fig. [35b) and the photon-gluon fusion (Fig. 155b ) processes. In LO the momentum 
fraction of the proton carried by the parton is given by £ = xsj(l + Mf 2 /Q 2 ), where 
xsj is the Bjorken scaling variable xsj = Q 2 /(Q 2 + W 2 ). Here W is the total cm. 
energy W 2 = (q + P) 2 , q = is the momentum of the virtual photon, Myx is the 
invariant mass of two jets of highest pr and Q 2 is the negative four-momentum 
transfer squared of the ingoing and outgoing electron. In the kinematic region of low 
Q 2 and low £, the 7*-gluon fusion dominates the jet production imparting sensitivity 
to the gluon component of the parton density functions (PDFs) of the proton, whereas 
the contribution of the QCD-Compton process yields information on the various quark 
(antiquark) components of the proton. 
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Fig. 34. Proton structure function ^(Q 2 ,^) given at two Q 2 values (3.5 GeV 2 and 90 
GeV 2 ). The various data sets have been normalised by the factors shown in the brackets in 
the key to the plot, which were determined in the NNLO MSTW2008 global analysis 12321. 
(From Amsler et al. in Ref. [2^5] ). 



In order to make theoretical predictions on jet productions in neutral current DIS 
scattering one needs the PDFs of the proton, provided mostly by the global analysis 
collaborations (i.e., analysis in which all high energy physics measurements relevant 
for testing QCD and determination of various parameters and non-perturbative func- 
tions are undertaken), such as CTEQ [23 1) and MSTW [232 . In addition one must 
have infrared and collinear safe parton cross sections, which are known now up to 
NLO in a s 233 234] . An example of the inclusive DIS measurements at HERA to- 
gether with data from lower energy fixed target experiments is shown in Fig. 1341 A 
striking feature of the HERA data is the dramatic rise of the proton structure func- 
tion F 2 (x, Q 2 ) for increasing Q 2 and low values of x (typically x < 10~ 2 ). Almost all 
of this rise of F 2 (x, Q 2 )is due to the gluon density in the proton. This has profound 
consequences for high energy pp (at the Tevatron) and pp collisions (at the LHC). 

DIS jet production depends in general on two large scales Q — \[Q 2 and the px of 
the produced jets. The ep jet data were collected by two detectors at HERA: HI and 
ZEUS, resulting from the collision of electrons or positrons with energy E e — 27.6 
GeV with protons of energy E p — 920 GeV, providing a center-of-mass energy of 
yfs ~ 320 GeV. In the more recent analysis the inclusive kr jet algorithm [75] is used 
to combine the particles in the hadronic final state into jets. Theoretical predictions 
(at next-to-leading order) have been corrected for hadronisation effects, which are 
calculated via Monte Carlo Models with parton showers. The most recent publication 
on jet production in DIS comes from the HI collaboration at HERA, where data up to 
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Fig. 35. The normalised inclusive jet cross sections measured as a function of the jet trans- 
verse momentum in the Breit frame Pt in two regions of Q 2 indicated on the frames. The 
points are shown at the average value of Pt within each bin. (From HI Collaboration 236 ). 



2007 are included and Q 2 spans the range 150 < Q 2 < 15000 GeV 2 236 . Inclusive jet, 
2-jet and 3-jet cross sections, normalised to the neutral current (NC) deep inelastic 
scattering cross section, are measured as functions of Q 2 , jet transverse momentum 
and proton momentum fraction. We show in Fig. [35] the normalised inclusive jet 
cross section as a function of the jet transverse momentum px in the Breit frame 
(defined as the frame in which 2xp + q = 0, where p and q are the 3-momenta of 
the proton and virtual photon, respectively) for two ranges of Q 2 , 700 < Q 2 < 5000 
GeV 2 (shown in the left-hand frame) and 5000 < Q 2 < 15000 GeV 2 (shown in the 
right-hand frame). Agreement between data and theoretical predictions 233 23"4"] is 
excellent. The ratio R (of data over theory) lies near 1 (shown at the bottom of these 
frames). Similar plots for bins with smaller Q 2 can be seen in .236 . HERA data 
on inclusive jet production in DIS 237 238 239 have constrained the gluon density 
in the range 0.01 < x < 0.1 The strong coupling a s {Q 2 ) has been determined and 
translated into a s (M z ) = 0.1168 ± 0.0007(exp.)±°°°^(theor.) ± 0.0016(PDF) using 
the usual renormalisation group equation. This result is competitive with those from 
e + e~ data and is in good agreement with the world average [235] . 

A similar recent analysis of the ZEUS collaboration has less integrated lumi- 
nosity, as it is based only on the data taken from 1998- 2000. However, their data 
include also results for rather large Q 2 . Therefore, Z exchange is included in ad- 
dition to the 7-exchange. The analysis is done in a similar fashion as that of the 
HI collaboration described above. They also studied the inclusive one-jet cross sec- 
tion as a function of Q 2 and E T e B (the transverse energy of the jet in the Breit 
frame). In addition they also measured this cross section for three different radii, 
R = 0.5,0.7 and 1.0, used for combining hadrons into jets with the help of the in- 
clusive kx cluster algorithm |79|80j . The results are shown in Fig. [35] for da/dQ 2 
for E T % > 8 GeV as a function of Q 2 . Further kinematic constraints, namely 
|cos7/j| < 0.65, where 7/, corresponds to the angle of the scattered quark in the 
quark-parton model and is constructed using the hadronic final state, and the pseu- 
dorapidity range —2 < rf B < 1.5, are indicated on the figure. The pseudorapidity 
variable, defined as r\ — — lntan(0/2) where cos0 = p z /\p\, is approximately equal to 
the (longitudinal-boost-invariant) variable called rapidity y = \ ln( E ^ ) in the limit 
\p\ 3> m, and can be measured when the mass m and the momentum of the particle \p\ 
are unknown. The NLO QCD predictions with scales p,R = E T e B , pf = Q, corrected 
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Fig. 36. Measured differential cross-section da/dQ 2 for inclusive-jet production with 
-®ts 8 GeV and —2 < i/ B °' < 1.5 (dots) for different jet radii, in the kinematic range 
given by | cos7 h | < 0.65. (From ZEUS Collaboration [230] ). 



to include hadronisation and Z effects, are compared with the measurements 240 . 
The calculations reproduce the measured differential cross section quite well for all 
three jet radii considered. In this work also a s (Q 2 ) has been determined. The result is 
a s (M z ) = 0.1207 ± 0.0014(stat.)tg °^(exp.);g ™g(theor.), which is also consistent 
with the world average. 

We now turn to photoproduction. At HERA the largest cross section is due to 
photoproduction, where the beam (electron or positron) interacts with the proton via 
the exchange of a virtual photon with a small virtuality Q 2 ~ 0. The spectrum of 
the ingoing virtual photon can very well be described by the well-known Weizsacker- 
Williams formula |24ll242l243j . 

The photoproduction of single jets, dijets and triple jets with high transverse 
momenta can be calculated also within perturbative QCD if the transverse momentum 
of the jets is large enough to provide the hard scale. Besides the larger cross section, 
as compared to the DIS jet production, the photoproduction of jets does not depend 
on the additional scale Q. The contributions to the theoretical cross sections which 
have been calculated up to NLO come from two processes: (i) the direct process in 
which the photon enters the hard sub-processes directly by coupling to the quarks, in 
the same way as in deep- inelastic ep scattering (see Fig. 155b. c in LO), and (ii) the so- 
called resolved process in which the photon fluctuates into partons, quarks or gluon, 
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and one of them participates in the hard parton-parton scattering process [2441245] . 
This latter process is equivalent to jet production in hadron-hadron collisions, of 
which the LO hard scattering cross sections for qq' — > qq' , gg — > gg and gq — > gq are 
written below. The only difference is that the PDFs of one of the hadrons is replaced 
by the photon PDF. This process, therefore, is sensitive to the parton structure of 
the proton and the photon. It is one of the few processes which can give information 
on the gluon content of the photon. 

The basic 7-parton processes which enter the calculation of the direct process in 
LO are the following: QCD Compton process: 79 — > gq, and the photon-gluon fusion: 
7<7 —¥ qq. These 7-parton cross sections have the following simple forms 

da 2 aa s li / u s\ 

"/q^gq-^ — ~ < 7: ---- , (73) 

da 9 a a s 1 / u t\ 

™^ q « : d^¥~ e k^2-s[i + i) < (74) 

where e q is the charge of the quark with flavour q, s = 4, i = —2(1 — cos^*), u = 
—2(1 + cos#*) and 8* is the angle of the dijets in their centre-of mass system. | cos 8* \ 
is related to the pseudorapidities of the two jets, 771 and 772 by 

I cos 0*| = |tanh(»?i -»ft)/2)| • (75) 

There are many observables which have been measured and which can be used to 
test the basic parton-parton cross sections for the direct and resolved process up to 
NLO [24612471248 249 . We shall present only a few taken from the most recent HI 
250 and ZEUS |251) publications. The da/d\ cos0*| distribution has been studied as 
a function of |cos#*| by the HI collaboration |250] with and without an additional 
cut on the invariant mass of the two jets Mjj for the direct (resolved) enhanced 
contribution. This analysis is done in terms of a variable x 7 defined by, 

1 2 

7 1=1 

where £t,i and £?t,2 are the transverse energies of the two jets with the two largest 
Et's. In LO the direct contribution is at x 1 = 1 and the resolved contribution has 
x 7 < 1. Therefore, by selecting events with x 1 > 0.8 (a; 7 < 0.8) the direct (resolved) 
parts of the cross sections are dominant. The results of the HI analysis are shown 
in Fig. [37] for the two bins of x 7 as a function of cos 0*, with the upper two frames 
without a cut on Mjj and the lower two frames with Mjj > 65 GeV. 

The cross section da /d\ cos 8* \ with the Mjj cut is sensitive to the dynamics of 
the underlying 7-parton and parton-parton hard interactions. The cross section in 
the resolved sample ie 7 < 0.8 rises more rapidly with | cos#*| than that in the direct 
sample due to the dominance of the virtual gluon exchange in the resolved processes 
(see formulae for parton-parton cross sections below). The dependence on | cos#*| and 
also the absolute normalisation are well predicted by the NLO calculations. Similar 
results have been presented by the ZEUS collaboration [252 by varying the dijet 
mass Mjj and their results have been presented in [253 . Another example is the 
cross section da/dEx, where Et is the mean transverse energy of the two jets 

Et = ^(Ef 1 + Ef 2 ) . (77) 

An example of such a cross section for x 1 > 0.75 and x 1 < 0.75, respectively, as 
presented by the ZEUS collaboration |251] is shown in Fig. [38] The cross section is 
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Fig. 37. Bin averaged cross sections as a function of | cos#*| for data (points), NLO QCD 
calculations with (solid line) and without (dashed line) hadronisation corrections <5had and 
for the PYTHIA Monte Carlo predictions (dotted line) scaled by a factor of 1.2. The inner 
(hatched) band of the NLOx (l+5h a d) result is the scale uncertainty, the outer (shaded) band 
is the total uncertainty. The cross sections are shown for two regions in x ly with and without 
an additional cut applied on the invariant dijet mass (Mjj). (From HI Collaboration 250 ). 




Fig. 38. Measured cross-section da/dE T for (a) x° ha > 0.75 and (b) a;° bs < 0.75 compared 
with NLO QCD predictions using the AFG04 [2"M] (solid line) and CJK [2"55] (dashed line) 
photon PDFs. The predictions using AFG04 are also shown with their associated uncer- 
tainties (shaded histogram). The ratios to the prediction using the AFG04 photon PDF are 
shown at the bottom of the figure. (From ZEUS Collaboration [251p . 
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measured up to Et — 80 GeV, i.e. further out in Et than in DIS jet production. In this 
figure also results for two different photon PDFs (namely, the so-called AFG04 [254 
and the CJK [255 ) are shown. The most sensitive cross section concerning direct and 
resolved separation is the cross section da/dx° bs , where x° hs is the x 7 defined above 
with the sum over the two jets (therefore "obs" in the notation of x 7 since not all 
jets are included in the sum). The result from the ZEUS collaboration is shown in 
Fig. 1391 from which one can see how the data compare with different photon PDFs 
assumed in the NLO prediction. An appreciable dependence on these PDFs is seen in 
the small x 1 region as one would expect. In an earlier analysis the ZEUS collaboration 
determined also the strong coupling a s , just from jet production in jp interactions 



alone. The result is a s [M z ) = 0.1224±0.0001(stat.)^;^g(exp.)t°;^(theor.) [2~51t] . 



and the variation of a s with the scale Et has been found in good agreement with the 
running of a s as predicted by QCD. 

Summarising the DIS and photoproduction processes at HERA, we see that QCD 
and jets have made very significant impact on the profile of the proton and the pho- 
ton in terms of their respective PDFs, which determine the luminosity functions of 
the parton-parton scatterings at high energies and hence the hard scattering cross 
sections of interest. Detailed studies of F2(x, Q 2 ) at HERA have also rekindled the- 
oretical interest in the small- x region. The evolution in ln(l/a;) at fixed value of 
Q 2 is governed by the so-called BFKL equation [2611262] , Originally developed to 
study Regge processes in high energy scatterings and the Pomeranchuk singularity 
(the QCD Pomeron), it can be combined with the DGLAP equation (for evolution 
in Q 2 ) to provide a quantitative description of the DIS structure functions over an 
enlarged (x,Q 2 ) domain. Several proposals in carrying out the small- a; resummation 
have been considered in the literature, which are comparatively discussed in a re- 
cent working group report [263) . In addition, the evolution in ln(l/x) leads to soft 
gluon enhancements, generating a dense gluonic system over a limited range of the 
nucleon wave function (hot spots). As the gluon occupation number becomes of order 
l/a s , non-linear effects present in the QCD Lagrangian become important, leading 
eventually to the saturation of the gluon density in the nucleons in high energy col- 
lisions [264] . This picture of high energy nucleonic wave functions (a high density, 
nonperturbative gluonic system with a weak coupling constant) is called the Color 
Glass Condensate [265 , and is of great interest in understanding the QCD aspects of 
heavy ion collisions, such as at RHIC and the LHC [266] . 



7.2 77 collisions 

Another area in which jet production has been studied experimentally and theoreti- 
cally is photon-photon collisions in the LEP2 energy range. The two incoming photons 
are produced in e + e~ collisions in the anti-tagged mode, i.e. when both the scattered 
electron and the positron escape detection. This is kinematically analogous to the 
photoproduction process in high energy ep collisions at HERA. In 77 — > hadrons, 
four classes of events have to be distinguished (see, Fig. [40]). The variables used in 
the classification of these events x^ and x~ , which are analogues of the variable x n 
in 7p collisions, are defined as follows: 
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Fig. 39. Measured cross section for da/dx°^ s compared with the NLO QCD predictions 
using the AFG04 [254] (solid line), CJK [255] (dashed line), AFG [257] (dotted line), 
GRV [258I259] (dashed and double-dotted line) and SAL [550] (dashed and single-dotted 
line) photon PDFs. The ratios to the prediction using the AFG04 [254] photon PDF are 
shown at the bottom of the figure. (From [251] V 



where 'part' corresponds to all detected particles and E- ]ct andp z j G t are the two hard- 
jets energy and the component of jet momentum along the z-axis, respectively. The 
four classes are: 

1. Hadron production via vector meson interactions (Vector- Meson Dominance Model 
VDM) [2ri7l2ri8] (Fig.ESt). 

2. The direct domain, where both x^ and x~ are close to 1. This domain is mostly 
populated by the quark-parton model like events 77 — > q + q (Fig. l4T)b). 

3. The single resolved domain, with the presence of a remnant jet, where only one of 
the x+ and x~ is close to 1 and the other is shifted to some lower value (Fig.l40h). 

4. The double-resolved domain, where both x^ are shifted to values below 1 (Fig.HQH). 

Due to the appearance of the double resolved region, jet production in 77 collisions has 
increased sensitivity to the gluon content of the resolved photon. This has enormous 
significance for future high energy 7 — 7 collisions, being entertained in the context 
of a high energy linear e + e~ collider. 

In the past, dijet production in 77 collisions has been studied experimentally at 
^fsZ from 189 to 209 GeV by the OPAL 269 and DELPHI [UD] collaborations 
at LEP. To that end a number of observables (differential distributions) have been 
measured by introducing a^-cuts at 0.75 (OPAL) and x± = 0.85 (DELPHI). These 
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Fig. 40. Main diagrams corresponding to the hadron production in 77 collisions via vector 
meson interactions (VDM-like, a), point-like interactions (QPM-like, b) and with one (c) or 
both (d) photons resolved into partons. (From |270] V 

distributions include, among others, da di]ct / dE^\ with £^ et = 1/2{E%\ + E^ 2 ) and 
dadijet/dXy. The data from both collaborations have been compared with the NLO 
QCD calculations based on the work of [246 and are found to be in good agreement. 
This is shown in Fig. [4]] for the differential distribution fda/dE 3 ^, where the factor 
/ is used to visibly separate the three measurements. Inclusive jet production in 7 — 7 
collisions has also been measured by the L3 |271j and OPAL |272j collaborations at 
LEP. 



7.3 Proton colliders 

7.3.1 Fundamental QCD scattering processes 

In parallel to electron and photon processes in QED, a large number of 2 — > 2 scatter- 
ing processes involving quarks and gluons are predicted in QCD, see e.g. |55] . They 
give rise to jets at hadron colliders. Most interesting are the fundamental abelian 
processes in QED transcribed to the non-abelian extensions in QCD, like 

Rutherford quark scattering : qq' — > qq , 
Rutherford gluon scattering : gg — ¥ gg , 
Super — Compton process : gq^? gq. 

Representative scattering diagrams are depicted in FigH5] The associated cross sec- 
tions scale in the energy squared s for massless initial and final-state quarks, while 
the angular distributions are given by 
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Fig. 41. The dijet cross-section in 77 collisions at LEP as a function of the mean transverse 
energy E 3 ^ of the dijet system, for the three regions in xt f — -space given in the figure. 
The factor / is used to separate the three measurements in the figure more clearly. The 
prediction of the LO program PYTHIA is compared to the data. The NLO calculation is 
from [243. (From OPAL [269] ). 
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and the variables s, u and u have been defined earlier. One should notice the three- 
gluon coupling already in LO. These amplitudes generate the expected Rutherford 
singularities ~ d0* 2 /8* 4 for forward scattering t — > 0, and analogously for backward 
scattering ii 0. 

Calculating the experimentally observed cross sections at hadron colliders requires 
three essential steps, which we have already outlined in the context of calculating the 
DIS cross sections, namely (i) the hard 2 — ^ 2 scattering processes, including NLO 
QCD corrections, (ii) flux of the incoming partons, determined in terms of the PDFs 
of the protons (and antiprotons), discussed earlier in the context of DIS scattering 
at HERA, and (hi) hadronic (non-perturbative) corrections. Here also QCD plays 
an important role in terms of the scale dependence of the PDFs and FFs. Thus, for 
example, the cross section of the hadron-hadron scattering with the four-momenta of 
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Fig. 42. Representative Feynman diagrams for fundamental QCD processes in hadronic 
collisions. 
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Fig. 43. Left-hand frame: Transverse energy distribution of jets at the Tevatron |273| . Right- 
hand frame: Comparison between the measured pr spectra by the CMS collaboration at the 
LHC and theory predictions for calorimeter jets [274]. For better visibility the spectra in 
both frames are multiplied by arbitrary factors indicated in the legend. 

the two colliding hadrons Pi and P2 can be written as [S3] 

(80) 

where the hard interaction between the partons i and j is given by &ij and p± — x\P\ 
and P2 = X2P2 are the momenta of the two partons. 

7.3.2 Jets in hadron colliders and tests of QCD 

An example for inclusive jet production at the Tevatron measured by the DO collab- 
oration is shown in Fig. [43] (left-hand frame) . Similar measurements have been done 
by the CDF collaboration at the Tevatron. The dominant contribution at small px 
can be traced back to Rutherford gluon scattering gg — » gg. This result is naturally 
expected since, on average, the gluon colour charges are significantly larger than the 
quark colour charges and, as discussed earlier, the gluon flux for low values of x by 
far exceeds the quark flux of high-energy protons. These jet cross sections can be 
exploited to determine the gluon distribution of the proton and to measure the QCD 
coupling 275J. By combination with other measurements the two observables are dis- 
entangled in the Tevatron measurements. The gluon flux extracted this way is large, as 
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anticipated, and the QCD coupling is compatible with the world average. This [273] . 
and related measurements |276I277| impact on the proton PDFs and have been used 
in updating this information [231 232 . In particular, they provide constraints on the 
gluon (and quark) distributions in the domain 0.01 < x < 0.5. A detailed discussion 
of jets and comparison of data and theory at the Tevatron can be seen in |278j . 

Very soon, similar but more sensitive analyses will also be undertaken at the LHC 
and a beginning has already been made. In Fig. [43] (right-hand frame), we show a 
comparison between the measured pr spectra by the CMS collaboration |274] at the 
LHC with y/s — 7 TeV and an integrated luminosity of 60 n& _1 for the calorimeter 
jets and theory predictions at the next-to-leading (NLO) order accuracy, using an 
anti-/cT jet algorithm with R = 0.5. Data are divided in several intervals of rapid- 
ity y bins. Theory predictions are based on NLOJET++ [375] with CTEQ-6.6 [2lH] 
sets of parton distribution functions (PDF). The non-perturbative (NP) corrections 
are estimated using two different hadronisation models, PYTHIA [215] and HER- 
WIG++ 280 , with the mean of the two predictions taken as the correction. Despite 
currently modest LHC luminosity, jets having transverse momenta up to 800 GeV are 
measured and the agreement with QCD is excellent. An in-depth review discussing 
the physics basis and use of the general purpose Monte Carlo event generators for 
hadronic collisions at the LHC is available |281] , to which we refer for a comprehensive 
discussion. 

Experiments at the LHC have opened a window to sub-energies in the TeV range 
for studying jet phenomena in QCD, enabling searches for physics beyond-the-SM in 
a number of such extensions. Both ATLAS and CMS have searched for new heavy 
particles, which manifest themselves as narrow resonances in their data collected at 
the LHC at yfs = 7 TeV. Such new states may include an excited composite quark 
q* , expected in theories with quark substructure 282 283 284 ; an axigluon predicted 
by chiral colour- models [2851286] ; a flavour- universal colour-octet coloron 287 288 ; 
or a colour-octet techni-p meson predicted by models of extended technicolor and 
topcolor-assisted tecnicolor [289I290I291T292] . 

The dijet invariant mass (rrijj) is an observable which is particularly sensitive 
to such new objects. This was studied already at the Tevatron in pp collisions with 
negative results, exemplified by the CDF limit on the mass of the excited quarks q* in 
which a mass range 260 < m r < 870 GeV was excluded at 95% C.L. [293 . ATLAS 
has extended this exclusion range to higher q* masses, with the range 0.40 < m q * < 
1.26 TeV now excluded using pp collisions [294 . Fig. l44l shows the predicted signal 
for q* masses of 500, 800, and 1200 GeV satisfying all event selection cuts. No signal 
of q* is found and the data are in excellent agreement with the background estimates 
based on the SM. 

Similar measurements of the dijet invariant mass spectrum and search for new 
particles decaying to dijcts have been performed by the CMS collaboration 295 . 
The highest observed dijet mass by CMS at yfs = 7 TeV is 2.13 TeV. No deviations 
are found from QCD up to this dijet mass. In particular, string resonances with a 
mass less than 1.67 TeV have been excluded by the current CMS measurements at 
95% C.L. The sensitivities to the narrow resonances in the dijet mass will increase 
substantially with the increase in the LHC luminosity and energy. For example, for 
the anticipated luminosity of 1 fb _1 at yfs = 7 TeV, the expected limits are all in the 
range of 2.5 to 3.5 TeV. 



7.3.3 Physics of the top quark and (W ± ,Z) bosons using jets 

Inclusive jet production in pp and pp collisions in association with a Z/j* /W boson 
provides a stringent test of QCD. As these final states are also of great importance in 
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Fig. 44. The data (D) dijet mass distribution (filled points) fitted using a binned back- 
ground (B) distribution described in the text (histogram). The predicted excited quark q* 
signals for excited quark masses of 500, 800 and 1200 GeV are overlaid, and the significance 
of the data-background difference is shown (from ATLAS collaboration [294]). 



the search of the SM Higgs boson arising from the process pp(p) —> W/Z + H(—> bb), 
and in the search of supersymmetry in the missing Et + jets channel, the processes 
PP(p) W/Z/ 7* + jets have received a lot of theoretical and experimental attention. 
In particular, theoretical predictions for vector boson production recoiling against a 
hadron jet at next-to-leading order were presented in [296 297 298 299]. The processes 
p + p — > W/Z/j* + 2 jets to the same level of theoretical accuracy were calculated for 
the Tevatron in [3001 and the corresponding processes p + p — > W/Z/-f* + 2 jets for 
the LHC in [301] , Vector boson production in association with n-jets for n < 4 was 
calculated in 302 303 . A parton- level event generator, called MCFM [304] . which 
gives theoretical predictions for a large number of processes containing W, Z and H 
bosons and jets (including heavy quark jets) is available for the Tevatron and the LHC 
colliders. Similar theoretical tools have been developed which give predictions for the 
transverse momentum distributions of the Z/j* /W produced in hadron collisions, 
based either on fixed order perturbation theory, such as [305 and 306 , or based on 
soft gluon resummations valid at low px 307 , such as RESBOS 308 . They have been 
used in conjunction with the PDFs (23 lj in the analysis of the Tevatron data [3091310] . 
and we show below representative measurements from the CDF Collaboration in 
Figs. 05] The NLO pQCD MCFM framework describes the data rather well over a 
large range of Pj, , as well as the jet- multiplicity. 

The production of heavy gauge boson pairs (WW, WZ, ZZ) in pp and pp collisions 
provides tests of the self-interactions of the gauge bosons and hence deviations from 
the SM-based predictions for the production rate could indicate new physics [311] . 
Since, topologically diboson production is similar to the associated Higgs boson pro- 
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Fig. 45. Left-hand frame: (top) Inclusive jet differential cross section measured by the 
CDF collaboration as a function of p!j? in Z/j* + > 1 jet events (black dots) compared to 
NLO pQCD predictions (open circles), (bottom) Data/Theory versus p^, . Right-hand frame: 
(top) Measured total cross section for inclusive jet production in — > /-i + fi~ events as a 
function of iVjet compared to LO and NLO pQCD predictions, (bottom) Ratio of data and 
LO pQCD predictions versus Nj c t- (From [310] ). 



duction pp(p) — > VH + X (V — W,Z), the experimental techniques developed in 
pp(p) — > VV are important for the Higgs boson searches as well. The process pp — > VV 
with both the vector mesons decaying into lepton pairs (W — > i^Uf,Z — > has 
been observed at the Fermilab Tevatron experiments by CDF [312 313] and DO [314] . 
Diboson production has not been conclusively observed in decay channels involving 
only hadrons. However, evidence for diboson decays into a mixed Iv^qq final state 
[i — e,[i,T;q = u,d,s,c,b) has been presented by DO [315 and CDF [316] . The 
experimental analyses involve large transverse momentum imbalance (due to the es- 
caped neutrino) and two jets whose invariant mass can be reconstructed. Because of 
the limited resolution in the dijet invariant mass, decays of — > 2 jets and Z —> 2 
jets are not distinguished separately. The most significant backgrounds to the diboson 
signals are W(£D)+ jets, Z(uD)+ jets and QCD multijet production. 

In Fig. 1461 we show the dijet mass distribution from the evqq and [ivqq channels 
for the DO data [315] and MC predictions. A clear diboson signal in the dijet invariant 
mass is seen in the lower frame. The resulting cross section <j(WV) = 20.2 ± 4.5 pb 
is consistent with the SM prediction a(WV) = 16.1 ± 0.9 pb at = 1.96 TeV [3T7] . 
Fig. [47] shows the corresponding measurements by CDF [316] . This yields a combined 
WW + WZ + ZZ cross section in pp collisions at \fs — 1.96 TeV: a(pp — > VV) 
In.!) ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) pb, consistent with the SM prediction. 

At the Fermilab Tevatron, top quarks are produced mostly in pairs pp — > ti+X. In 
the SM, top quarks decay into a W boson and a b quark almost 100% of the time. The 
topology of the final states resulting from the ti production depends on whether the 
W boson decays leptonically W — > Ivt, or hadronically W — > qq' leading to two jets. 
Following this, it events have been measured in dilepton l + (r + X, single lepton ^ ± +4 
jets and also in the non-leptonic mode with no energetic leptons. The non-leptonic ti 
final state has the advantage of a large branching ratio (~ 4/9). The major challenge 
of this channel is the large background from QCD multijet production. To increase 
the purity of the candidate sample, methods based on artificial neural networks are 
applied to the data. Further improvement is then obtained from the requirement of 
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Fig. 46. (a) The dijet mass distribution from the combined evqq and \ivqq channels for 
data from the DO collaboration at the Tevatron and MC predictions, (b) A comparison of 
the extracted signal (filled histogram) to background-subtracted data (points), along with 
the ±lcr systematic uncertainty on the background. The residual distance between the data 
points and the extracted signal, divided by the total uncertainty, is shown at the bottom 
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at least one jet identified as originating from a b quark using a secondary vertex b- 
tagging algorithm. These techniques have allowed one to measure the top quark mass 
and the it cross section in spite of the overwhelming QCD multijet production. 

To these ends, a reconstructed top quark mass, m\ cc , is determined by fitting the 
kinematics of the six leading jets from the process pp — > tt+X — > 6 jets. There exists a 
strong correlation between m\ ec and the jet energy scale JES. However, the JES can be 
calibrated using a selected sample of ti candidate events, where a second variable m^- c 
is reconstructed from the jets assigned to the W boson. The variable m^ c is related 
to the boson mass, which is known accurately. Relating m\ ec and m^ c to match 
the experimental data {in situ calibration) reduces significantly the systematic errors. 
Further improvement comes by using a multivariate approach taking advantage of the 
distinctive features of the signal and background events through a neural network. 
Fig. 25] shows the histogram of m\ cc as obtained in the data and compared to the 
distributions in the so-called 1-tag and > 2-tag events from signal and background 
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Fig. 47. Top: Comparison between data in pp ->• VV + X) (V = W ± , Z°) from the CDF 
collaboration at the Fermilab Tevatron at i/s = 1.96 TeV and the fitted background only 
in the Dijet invariant mass. Bottom: Comparison of the diboson signal (solid line) with 
the background subtracted data (points). The dashed lines represent the ±ler statistical 
variations on the signal [CDF |316| ]. 

corresponding to M top = 175 GeV. The best estimates of the top quark mass from 
this analysis is [318] 

M top = 174.8 ± 2.4(stat + JES) GeV . (81) 

The procedure used to measure the top quark mass also returns the average number 
of signal events expected, given the selected data samples. These results can be turned 
into a measurement of the ti cross section, and yield 

o t i = 7.2 ± 0.5(stat) ± l.O(syst) ± 0.4(lum) pb . (82) 



7.3.4 Searches for the Higgs particles 

We have discussed numerous electroweak processes at the Tevatron in which jets play 
an essential role in the analysis. In particular, and Z gauge bosons and top quarks 
have been measured using jets. The last and the most-prized on this list is the Higgs 
boson. This is being searched for at the Tevatron feverishly. For m# < 135 GeV, the 
dominant decay mode is H — > bb [319] ; analyses of this decay mode open a powerful 
new Higgs discovery channel |320) . The dominant production modes are gg — > H and 
qq — > H. The bb signal in this channel is overwhelmed by the QCD bb production. A 
promising production and search strategy is the production of a Higgs boson decaying 
to a pair of bottom quarks in association with a vector boson V (W or Z) decaying to 
quarks or leptons, resulting into a four-jet or a charged lepton + two-jet final states. 
In either case, two of the jets are required to have secondary vertices consistent with 
B-hadron decays. So far Tevatron Run II searches have used signatures where the 
V decays to leptons (see, for example Refs. [3211322] ). Recently, also searches in the 
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four-jet channels have been reported |323j . Using an integrated luminosity of 2 fb _1 , 
Higgs boson searches in this channel provide a weak upper bound. For example, for 
run = 120 GeV, CDF is able to exclude a Higgs production cross section larger 
than 38 times the SM prediction! Hence, establishing the Higgs signal in this channel 
requires much more statistics, but also some fundamental progress in jet algorithms 
to be more efficient in Higgs (and other similar particle) searches. 
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New opportunities are offered by observing b jets in Higgs decays at the LHC. 
The key technique is the 2-jet splitting of a fat bb jet generated in events in which the 
Higgs boson is boosted to large transverse momenta in the processes pp — > W^H and 
ZH . If the "fat" jet is characterised by a jet radius R = ~ Mh/pt, the clustering 
is partially undone by shrinking the radius R until the fat jet R b i decomposes into two 
slim Rb subjets with significantly lower mass, each containing a 6 quark. Additional 
criteria will reduce the contamination by standard QCD processes. Though the boost 
will strongly reduce the event rate, the significance will nevertheless be raised to such 
a level that light Higgs events can clearly be isolated above background. Extending 
the method to the channel pp — > tiH — > tibb, the crucial ttH coupling, apparently 
not accessible at LHC otherwise, can be measured in the light Higgs sector [324] . 

The concept is useful also for the analysis of other processes, for example, the 
search for supersymmetric particles decaying to jets from the hadronic decays of 
the electroweak and Higgs bosons [3 2 5) , or for detecting strongly interacting W , Z 
bosons |326| . or the search for heavy resonances in decays to top-quark jets [327 . 



8 Summary 

Quantum Chromodynamics has been established experimentally in the past four 
decades as the microscopic theory of the strong interactions, formulated as a non- 
abelian gauge theory for coloured quarks and gluons within the Standard Model of 
particle physics. Jet physics has been a crucial instrument for achieving this funda- 
mental result. The beginning was made at SPEAR with the observation of quark jets 
in e + e~ annihilation by the SLAC-LBL collaboration |21) . Subsequent studies un- 
dertaken, in particular at DORIS, PEP and PETRA, involving higher center-of-mass 
energies largely consolidated the phenomenological profile of the quark jets (see [54] 
for a review). In fact, jets provide an irrefutable case for the existence of quarks as dy- 
namical entities directly observable in particle detectors, despite colour confinement, 
convincing even the most die-hard skeptics about the reality of quarks. Moreover, 
making use of the larger masses of the b- and c-quarks, relatively long half-lifes of 
the corresponding hadrons and their characteristic decay patterns, one can efficiently 
flavour-tag the heavy quark jets. In the meanwhile, these techniques have been de- 
veloped to the level of a diagnostic tool to search for new phenomena in which heavy 
quarks play a role. The decays t — > bW and H — » bb are two good cases in point. 

Theoretically, the existence proof of quark jets in fixed order perturbative QCD 
was provided by Sterman and Weinberg |45j using a jet-cone definition which coin- 
cided with the actual process of detection of hadrons in finite segments of hadron 
calorimeters. Subsequently, Sterman |46) provided an all-orders argument for the 
infra-red safety of jet cross sections. Phenomenologically, quark jets follow from the 
observation that the transverse momenta of the hadrons produced in the fragmen- 
tation q* —¥ hadrons is limited, whereas the longitudinal components of the hadron 
momenta scale with energy. A very intuitive and largely accurate quark jet fragmenta- 
tion model was developed along these lines by Field and Feynman 37 , which played 
an important role in the quantitative analysis of jet data. 

Analysis of the decays T(9.46) — > hadrons measured by the PLUTO collabora- 
tion |24| working at DORIS were undertaken in terms of the underlying perturbative 
process T(9.46) — > ggg. In particular, the experimental profile of the most ener- 
getic parton ((£i) ~ 4.1 GeV) was close to the phenomenological expectations of 
a hadron jet. However, a clear three-jet topology using en vogue jet definitions was 
not established in T(9.46) decays for lack of energy [21]. This three-jet topology 
was established later in the e + e~ experiments operating at higher energies (typically 
^/s = 30 GeV), resulting from the energetic quark, anti-quark and gluon from the 



70 



process e + e~ — > qqg. The jet profile of the three jets as well as the inclusive hadronic 
measurements undertaken at PETRA in 1979 followed detailed theoretical expecta- 
tions 30 31 321. Thus, it is fair to conclude that the study of the decay T(9.46) —> 3g, 
initiated by PLUTO at DORIS [23], was an important step in the confirmation of 
QCD which served as a prelude to the unambiguous discovery of the three-jet topology 
by the experiments at PETRA. More detailed and quantitative tests of perturbative 
QCD in the decays of T(9.46) were also presented subsequently by PLUTO [111] . 
Theoretical proofs of the existence of three-jet topologies, in the Sterman- Weinberg 
sense, were provided in 1981, and somewhat later in terms of the next-to- leading 
order calculations of the three-jet cross sections. This was done for inclusive jet dis- 
tributions (such as the Fox- Wolfram shape variable [1421143] . thrust [1441145] and 
energy-energy correlations |153|154|157ll58"] ) and in terms of topological jet cross sec- 
tions [163 164 165 166 167 169 . Confirmation of the non-abelian character of QCD in 
jets [218121912201221] in the four-parton processes e + e~ — > qqgg came from experi- 
ments at LEP 222 223 . In the meanwhile, multijet physics has developed enormously, 
with the NLO calculation of e + e~ — > 7, Z — > 4 jets [210121 1I2T2"] completed around 
1996, and the NLO calculation to five-jet production at LEP reported recently [214] . 

The properties of gluon jets have largely been determined by experiments at PE- 
TRA and subsequently at LEP. The fragmentation of gluon jets was initially conceived 
by treating them as independent partons [3T], or implemented by the perturbative 
process g* — > qq as the first step followed by incoherent quark fragmentation [32] (U 
models) . The resulting picture could largely account for the essential properties of the 
two- and three-jet events seen at PETRA and PEP, and they helped in the discovery 
of three-jet events (and hence, of gluons) at PETRA. However, analysis of the PETRA 
jets saw the emergence of an alternative fragmentation scheme for the e + e~ —¥ qqg 
events - the LUND string model [38] - in which hadronisation was implemented in 
terms of two strings stretched along the colour-anticolour axes which then fragmented 
very much like the quark-antiquark string in e + e~ — » 2-jets. This model provided a 
better phenomenological description of data, in particular the particle flow between 
the quark, antiquark and gluon jets [179118011811182] . The LUND-string effect was 
subsequently understood in perturbation theory in terms of the antenna radiation 
pattern of QCD [48] , reflecting the colour coherence effect of the non-abelian charac- 
ter of this theory. Detailed fragmentation models were built along the angle-ordered 
perturbation theory, which preserve the colour coherence in QCD, and in which par- 
ton showers were included in the form of cascades, which then finally fragmented as 
hadron clusters according to phase space (cluster hadronisation models) [3S]. These 
Monte Carlo models developed for the PETRA jet analysis have played an important 
role in the analysis of all high energy data involving jets. The modern incarnation 
of these fragmentation models are PYTHIA @£], HERWIG [EE] and SHERPA [52] , 
which differ in details on how the parton showers are matched on to the fixed order 
perturbative QCD matrix elements and in the hadronisation schemes. A central role is 
also played by the jet algorithms, which starting from the JADE scheme [H] have now 
evolved as trustworthy tools in the definition of jets, with the modern versions called 
the kr [78] (mostly in e + e~ annihilation processes) and anti-fcr jet algorithms [82] , 

Another large application of QCD is in studying DIS, photoproduction and 77 
collisions. In these cases, initial states are not so well known as in e + e~ annihilation. 
Jets and QCD have played a central role in mapping the PDFs of the proton and 
the photon. We have summarised some of the highlights in this article. In particu- 
lar, DIS measurements at HERA [235I236I237I238I239I24U] have firmly established 
the rise of the structure function F2(x, Q 2 ), which is due to the rapid growth of the 
gluon density g{x,Q 2 ) for low values of x as Q 2 increases. Likewise, high energy pp 
collisions at the Tevatron, in particular the Tevatron Run II data on inclusive jet pro- 
duction [273 276 272, have led to greatly firming up the PDFs of the proton. On the 
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theoretical side, the complete next-to-next-to leading order (3-loop) parton splitting 
functions have been derived by Moch et at [328 3291 . They have been used in working 
out the proton PDFs by the CTEQ [221] and the MSTW collaborations. Thus, 
the HERA and the Tevatron measurements and the progress in the QCD splitting 
functions will prove to be an asset in understanding the forthcoming LHC data. 

In the meanwhile, a fundamental change of paradigm has taken place concerning 
"Jets and QCD". The theory (QCD) is so well controlled (in particular, a s (Mz) 
is known to an accuracy of better than 1% and the crucial property of asymptotic 
freedom is now fully established) that jet-physics can serve as a tool to chart out 
new territories in high energy physics. We have reviewed here some applications of 
jet techniques in quantifying the properties of the top quark and the electroweak 
gauge bosons ( W, Z) . They have already played a significant role in determining the 
properties of the SM particles in experiments at the Tevatron and they will play an 
even more important role in the analysis of data from the experiments at the LHC. For 
example, jets are now increasingly used in developing search strategies for the Higgs 
boson, and even particles suggested in theories beyond the Standard Model. New jet 
techniques, such as the 2-jet splitting of a fat bb jet will be required to disentangle 
the decay H — > bb from an overwhelming QCD background in hadron colliders. 

One problem in jet physics however remains unsolved up to now. While, due to 
asymptotic freedom, the dynamics of quarks and gluons can theoretically be described 
with high accuracy at short distances, matched by numerical lattice calculations for 
static properties of hadrons, the transition from small to large distances in the evolu- 
tion of jets is theoretically not understood. However, a bridge is built, at the present 
time, by intuitively formulated models, which are constrained experimentally so strin- 
gently that hadron jets can be exploited to draw a valid picture of quarks and gluons 
and their interactions at short distances. New theoretical methods may help solve 
this multi-scale problem rigorously in the future. 

Acknowledgements We thank Peter Zerwas for the collaboration in the early 
stages of this work, for numerous helpful discussions that we had with him all along, 
and for his valuable input in good parts of this manuscript. Helpful discussions with 
Hans-Jiirgen Meyer, Hinrich Meyer and Bruno Stella on the PLUTO analysis of the 
T(9.46) data are thankfully acknowledged. We also thank a large number of our 
colleagues and collaborators whose dedicated and painstaking work over decades has 
contributed decisively to the development of QCD and jet physics. This article is 
dedicated collectively to all of them. 



References 

1. H. Fritzsch and M. Gell-Mann, Current Algebra: Quarks and What Else?, Proceedings, 
16th Int. Conference on High Energy Physics, Batavia IL, eConf C720906V2: 135-165 
(1972) 

2. H. Fritzsch, M. Gell-Mann and H. Leutwyler, Advantages of the Color Octet Gluon 
Picture, Phys. Lett. B 47,365 (1973) 

3. D. J. Gross and F. Wilczek, Ultraviolet Behavior of Non-Abelien Gauge Theories, Phys. 
Rev. Lett. 30, 1343 (1973) 

4. H. D. Politzer, Reliable Perturbative Results for Strong Interactions?, Phys. Rev. Lett. 
30, 1346 (1973) 

5. M. Gell-Mann, A Schematic Model of Baryons and Mesons, Phys. Lett. 8, 214 (1964) 

6. G. Zweig, An SU(3) Model for Strong Interaction Symmetry and its Breaking, CERN- 
TH-401 (1964) 

7. K. Nakamura et ai, [Particle Data Group], Review of Particle Physics, J. Phys. G 37, 
075021 (2010) 



72 



8. O. W. Greenberg, Spin and Unitary Spin Independence in a Paraquark Model of Baryons 
and Mesons, Phys. Rev. Lett. 13, 598 (1964) 

9. M. Y. Han and Y. Nambu, Three-triplet Model with Double SU(3) Symmetry, Phys. Rev. 
139, B1006 (1965) 

10. M. Gell-Mann, Quarks, Acta Phys. Austr. Suppl. IX, 733 (1972) 

11. Y. Nambu, A Systematics of Hadrons in Subnuclear Physics, in Preludes in theoretical 
physics, North-Holland Publishing Company, Amsterdam (1966) 

12. C. N. Yang and R. L. Mills, Conservation of Isotopic Spin and Isotopic Gauge Invariance, 
Phys. Rev. 96, 191 (1954) 

13. K. G. Wilson, Confinement of Quarks, Phys. Rev. D 10, 2445 (1974) 

14. C. H. Llewellyn Smith, Inelastic Lepton Scattering in Gluon Models, Phys. Rev. D 4, 
2392 (1971) 

15. D. J. Gross, How to Test Scaling in Asymptotically Free Theories, Phys. Rev. Lett. 32, 
1071 (1974) 

16. R. P. Feynman, The Behavior of Hadron Collisions at Extreme Energies, in Proceedings 
of the 3rd Topical Conference on High Energy Collision of Hadrons, Stony Brook, N. Y., 
237 (1969) 

17. R. M. Barnett et al. [Particle Data Group], Review of particle physics. Particle Data 
Group, Phys. Rev. D 54, 1 (1996) 

18. N. Cabibbo, G. Parisi and M. Testa, Hadron Production in e + e~ Collisions, Lett. Nuovo 
Cim. 4S1, 35 (1970) [Lett. Nuovo Cim. 4, 35 (1970)] 

19. S. D. Drell, D. J. Levy and T. M. Yan, Theory of Deep-Inelastic Lepton-Nucleon Scat- 
tering and Lepton-Pair Annihilation Processes, iii: Deep-Inelastic Electron-Positron Anni- 
hilation, Phys. Rev. D 1, 1617 (1970) 

20. S. M. Berman, J. D. Bjorken and J. B. Kogut, Inclusive Processes at High Transverse 
Momentum, Phys. Rev. D 4, 3388 (1971) 

21. G. Hanson et al. [SLAC-LBL Collaboration], Evidence for Jet Structure in Hadron 
Production by e+e~ Annihilation, Phys. Rev. Lett. 35, 1609 (1975) 

22. G. Hanson et al. [SLAC-LBL Collaboration], Hadron Production by e + e~ Annihilation 
at Center-of-Mass Energies between 2.6 GeV and 7.8 GeV, Part 2: Jet Structure and 
Related Inclusive Distributions, Phys. Rev. D 26, 991 (1982) 

23. J. D. Bjorken and S. J. Brodsky, Statistical Model for Electron-Positron Annihilation 
into Hadrons, Phys. Rev. D 1, 1416 (1970) 

24. Ch. Berger et al. [PLUTO Collaboration], Jet Analysis of the Upsilon(9.46) Decay into 
Charged Hadrons, Phys. Lett. B 82, 449 (1979) 

25. R. Brandelik et al. [TASSO Collaboration], Evidence for Planar Events in e + e~ Anni- 
hilation at High Energies, Phys. Lett. B 86, 243 (1979) 

26. D. P. Barber et al. [Mark- J Collaboration], Discovery of Three Jet Events and a Test of 
Quantum Chromodynamics at Petra Energies, Phys. Rev. Lett. 43, 830 (1979) 

27. C. Berger et al. [PLUTO Collaboration], Evidence for Gluon Bremsstrahlung in e + e~ 
Annihilations at High-Energies, Phys. Lett. B 86, 418 (1979) 

28. W. Bartel et al. [JADE Collaboration], Observation of Planar Three Jet Events in e + e _ 
Annihilation and Evidence for Gluon Bremsstrahlung, Phys. Lett. B 91, 142 (1980) 

29. B. H. Wiik, First Results From Petra, in *Bergen 1979, Proceedings, Neutrino '79, 
Vol.1*, p. 113 

30. J. R. Ellis, M. K. Gaillard and G. G. Ross, Search for Gluons in e + e~ Annihilation, 
Nucl. Phys. B 111,253 (1976) [Erratum-ibid. B 130, 516 (1977)] 

31. P. Hoyer, P. Osland, H. G. Sander, T. F. Walsh and P. M. Zerwas, Quantum Chromo- 
dynamics and Jets in e+e~ , Nucl. Phys. B 161, 349 (1979) 

32. A. Ali, E. Pietarinen, G. Kramer and J. Willrodt, A QCD Analysis of the High-Energy 
e+e" Data from Petra, Phys. Lett. B 93, 155 (1980) 

33. M. L. Mangano and T. J. Stelzer, Tools for the simulation of hard hadronic collisions, 
Ann. Rev. Nucl. Part. Sci. 55, 555 (2005) 

34. S. Bethke, QCD studies at LEP, Phys. Rept. 403-404, 203 (2004) 

35. A. Breakstone et al. [Ames-Bologna-CERN-Dortmund-Heidelberg- Warsaw Collabora- 
tion], High Pt Hadrons as Leading Particles in Jets Produced at the ISR. 1. Momentum 
Distribution of Secondaries in the Trigger Jet, Z. Phys. C 23, 9 (1984) 



73 



36. W. G. Scott [UA1 Collaboration], Jets in the UA1 Experiment, In *Berne 1984, Pro- 
ceedings, Proton Antiproton Collider Physics*, p. 6 

37. R. D. Field and R. P. Feynman, A Parametrization of the Properties of Quark Jets, 
Nucl. Phys. B 136, 1 (1978) 

38. B. Andersson, G. Gustafson and T. Sjostrand, A Three-Dimensional Model for Quark 
and Gluon Jets, Z. Phys. C 6, 235 (1980) 

39. G. Marchesini and B. R. Webber, Simulation of QCD Jets including Soft Gluon Inter- 
ference, Nucl. Phys. B 238, 1 (1984) 

40. S. L. Wu and G. Zobernig, A Method of Three Jet Analysis in e + e~ Annihilation, Z. 
Phys. C 2, 107 (1979) 

41. K. Lanius, Cluster Methods of Jet Analysis, DESY Report 80/36 (1980) 

42. K. Lanius, H. E. Roloff and H. Schiller, Selection of Jets in Multi - Hadron Final States 
Produced by e+e" Annihilation, Z. Phys. C 8, 251 (1981) 

43. H. J. Daum, H. Meyer and J. Burger, A Cluster Algorithm for Jet Studies, Z. Phys. C 
8, 167 (1981) 

44. S. Bethke et al. [JADE Collaboration], Experimental Investigation of the Energy De- 
pendence of the Strong Coupling Strength, Phys. Lett. B 213, 235 (1988) 

45. G. Sterman and S. Weinberg, Jets from Quantum Chromodynamics, Phys. Rev. Lett. 
39, 1436 (1977) 

46. G. F. Sterman, Phys. Rev. D19, 3135 (1979) 

47. B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Parton Fragmentation and 
String Dynamics, Phys. Rept. 97, 31 (1983) 

48. Y. I. Azimov, Y. L. Dokshitzer, V. A. Khoze and S. I. Troian, The String Effect and 
QCD Coherence, Phys. Lett. B 165, 147 (1985) 

49. T. Sjostrand, PYTHIA 8 Status Report. larXiv:0809.0303l [hep-ph] 

50. T. Sjostrand, S. Mrenna and P. Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput. 
Phys. Commun. 178, 852 (2008) 

51. M. B&hr et al., Herwig++ Status Report, Eur. Phys. J., C58 (2008) 639. larXiv:0 809. 2002 
[hep-ph] (2008) 

52. T. Gleisberg, S. Hoche, F. Krauss, M. Schonherr, S. Schumann, F. Siegert and J. Winter, 
Event Generation with SHERPA 1.1, JHEP 0902, 007 (2009) 

53. G. Kramer, Theory of Jets in Electron-Positron Annihilation, Springer Tracts in Modern 
Physics, Nr. 102 (1984) 

54. A. Ali and P. Soding, High Energy Electron-Positron Physics, World Scientific (1988) 

55. R. K. Ellis, W. J. Stirling and B. R. Webber, QCD and Collider Physics, Cambridge 
University Press (1996) 

56. G. S. Bali, H. Neff, T. Duessel, T. Lippert and K. Schilling [SESAM Collaboration], 
Observation of String Breaking in QCD, Phys. Rev. D 71, 114513 (2005) 

57. G. Abbiendi et al. [OPAL Collaboration], Leading particle production in light flavor 
jets, Eur. Phys. J. C 16, 407 (2000) [arXiv:hep-ex/0001054| 

58. S. Albino, B. A. Kniehl and G. Kramer, Fragmentation functions for light charged 
hadrons with complete quark flavor separation, Nucl. Phys. B 725, 181 (2005) 
[arXiv:hep-ph/0502188] . 

59. D. H. Saxon, in Ref. [54] 

60. S. Albino, The Hadronization of Partons, Rev. Mod. Phys. 82, 2489 (2010) 

61. Ya. I. Azimov, L. L. Frankfurt and V. A. Khoze, On Reaction e + e~ — > Hadrons, New 
Particles and e + e~ Annihilation (in Russian), Leningrad 76-222 (1976) 

62. J. D. Bjorken, Properties of Hadron Distributions in Reactions Containing Very Heavy 
Quarks, Phys. Rev. D 17, 171 (1978) 

63. C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, Scaling Violations in Inclusive 
e+e" Annihilation Spectra, Phys. Rev. D 27, 105 (1983) 

64. B. Mele and P. Nason, The Fragmentation function for heavy quarks in QCD, Nucl. 
Phys. B 361, 626 (1991). 

65. J. P. Ma, Perturbative prediction for parton fragmentation into heavy hadron, Nucl. 
Phys. B 506, 329 (1997) [arXiv:hep-ph/9705446] 



74 



66. P. Nason and C. Oleari, A Phenomenological Study of Heavy-Quark Fragmentation 
Functions in e+e~ Annihilation, Nucl. Phys. B 565, 245 (2000) 

67. M. Cacciari, P. Nason and C. Oleari, A Study of Heavy Flavoured Meson Fragmentation 
Functions in e+e" Annihilation, JHEP 0604, 006 (2006) 

68. T. Kneesch, B. A. Kniehl, G. Kramer and I. Schienbein, Charmed-meson fragmentation 
functions with finite-mass corrections, Nucl. Phys. B 799, 34 (2008) |arXiv:071 2.0481 [hep- 
Ph]] 

69. B. A. Kniehl, G. Kramer, I. Schienbein and H. Spiesberger, Finite-mass effects on inclu- 
sive B meson hadroproduction, Phys. Rev. D 77, 014011 (2008) |arXiv:0705.4392l [hep-ph]] 

70. D. Amati and G. Veneziano, Preconfinement as a Property of Perturbative QCD, Phys. 
Lett. B 83, 87 (1979) 

71. V. V. Sudakov, Vertex Parts at Very High-Energies in Quantum Electrodynamics, Sov. 
Phys. JETP 3, 65 (1956) [Zh. Eksp. Teor. Fiz. 30, 87 (1956)] 

72. G. Alexander Proceedings of the XlXth International Conference On High Energy 
Physics (Tokyo, 1978) p. 255 

73. S. Brandt, C. Peyrou, R. Sosnowski and A. Wroblewski, The Principal Axis of Jets: An 
Attempt to Analyze High-energy Collisions as Two-body Processes, Phys. Lett. 12, 57 
(1964) 

74. E. Fahri, A QCD Test for Jets, Phys. Rev. Lett. 39, 1587 (1977) 

75. S. Brandt and H. J. Dahmen, Axes and Scalar Measures of Two- Jet and Three- Jet 
Events, Z. Phys. CI, 61 (1979) 

76. O. Nachtmann and A. Reiter, A New Quantity for Finding Four- Jet Events in Electron- 
Positron Annihilation into Hadrons, Z. Phys. C 14, 47 (1982) 

77. G. P. Salam, Towards Jetography, Eur. Phys. J. C 67, 637 (2010) 

78. S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock and B. R. Webber, New Clustering 
Algorithm for Multi - Jet Cross-sections in e + e~ Annihilation, Phys. Lett. B 269, 432 
(1991) 

79. S.D. Ellis and D.E. Soper, Successive Combination Jet Algorithm for Hadron Collisions, 
Phys. Rev. D 48, 3160 (1993) 

80. S. Catani et al, Longitudinally Invariant kt Clustering Algorithms for Hadron Hadron 
Collisions, Nucl. Phys. B 406, 187 (1993) 

81. M. Wobisch and T. Wengler, Hadronization Corrections to Jet Cross Sections in Deep- 
Inelastic Scattering, |arXiv:hep-ph/9907280] (1999) 

82. M. Cacciari, G. P. Salam and G. Soyez, The Anti-fc t Jet Clustering Algorithm, JHEP 
0804, 063 (2008) 

83. G. C. Blazey et al., Run II Jet Physics, |arXiv:hep-ex/0005012j (2000) 

84. J. E. Huth et al, Toward a Standardization of Jet Definitions, FERMILAB-CONF-90- 
249-E (1990), Published in Snowmass Summer Study 1990:0134-136 (QCD161:D15:1990) 

85. M. H. Seymour, Jet shapes in hadron collisions: Higher orders, resummation and 
hadronization, Nucl. Phys. B 513, 269 (1998) |arXiv:hep-ph/9707338] 

86. R. P. Feynman, Photon- Hadron- Interactions (Benjamin, Reading) (1972) 

87. R. Gatto, Proceedings of the International Symposium, on Electron and Photon Interac- 
tions at High Energies (Deutsche Physikalische Gesellschaft e. V., Hamburg 1965) edited 
by G. Hohler, G. Kramer and U. Meyer-Berkhout, Vol. 1,106 (1965) 

88. Ch. Berger et al. [PLUTO Collaboration], A Study of Jets in Electron Positron Anni- 
hilations into Hadrons in the Energy Range 3.1 GeV to 9.5 GeV, Phys. Lett. B 78, 176 
(1978) 

89. E. Elsen, Multi-Hadron Production in e + e~~ Annihilation at Petra Energies and Com- 
parison with Quantum Chromodynamics (in German), Interner Bericht DESY F22-81/02 
(1981) 

90. C. G. Callan, D. J. Gross, High-Energy Electroproduction and the Constitution of the 
Electric Current, Phys. Rev. Lett. 22, 156 (1969) 

91. T. D. Lee and M. Nauenberg, Degenerate Systems and Mass Singularities, Phys. Rev. 
133, B1549 (1964) 

92. P. M. Stevenson, Comments on Sterman- Weinberg Jet Formula, Phys. Lett. B 78, 451 
(1978) 



75 



93. G. Kramer and B. Lampe, Two- Jet Cross-Section in e + e Annihilation, Z. Phys. C 34, 
497 (1987) [Erratum-ibid. C 42, 504 (1989)] 

94. W. Bartel et al. [JADE Collaboration], Experimental Studies on Multi-Jet Production 
in e+e~ Annihilation at PETRA Energies, Z. Phys. C 33, 23 (1986) 

95. S. W. Herb et al, Observation of a Dimuon Resonance at 9.5 GeV in 400 GeV Proton- 
Nucleus Collisions, Phys. Rev. Lett. 39 (1979) 252 

96. W. R. Innes et al, Observation of Structure in the Upsilon Region, Phys. Rev. Lett. 39, 
1240 (1979) 

97. Ch. Berger et al. [PLUTO Collaboration], Observation of a Narrow Resonance in e + e~ 
Annihilation at 9.46 GeV, Phys. Lett. B 76, 243 (1978) 

98. Ch. Berger et al. [PLUTO Collaboration], Determination of the Electronic Branching 
Ratio of the Upsilon(9.46) and an Upper Limit for its Total Width, Phys. Lett. B 93, 497 
(1980) 

99. C. W. Darden et al, Observation of a Narrow Resonance at 9.46 GeV in Electron- 
Positron Annihilation, Phys. Lett. B 76, 246 (1978) 

100. J. K. Bienlein et al, Observation of a Narrow Resonance at 10.02 GeV in e + e~ Anni- 
hilations, Phys. Lett. B 78, 360 (1978) 

101. D. Andrews et al. [CLEO Collaboration], Observation of Three Upsilon States, 
Phys. Rev. Lett. 44, 1108 (1980) 

102. T. Bohringer et al, Observation of Upsilon, Upsilon' and Upsilon" at the Cornell 
Storage Ring, Phys. Rev. Lett. 44, 1111 (1980) 

103. T. Appelquist and H. D. Politzer, Orthocharmonium and e + e~ Annihilation, 
Phys. Rev. Lett. 34, 43 (1975) 

104. T. Appelquist and H. D. Politzer, Heavy Quarks and Long-lived Hadrons, Phys. Rev. D 
12, 1404 (1975) 

105. K. Roller and T. F. Walsh, Three Gluon Jets as a Test of QCD, Phys. Lett. B 72, 227 
(1977) 

106. K. Roller and T. F. Walsh, Gluons in Quarkonium Decay, Nucl. Phys. B 140, 449 
(1978) 

107. R. Roller, H. Rrasemann and T. F. Walsh, The 3 Gluon Decay of Quarkonium, 
Z. Phys. C 1, 71 (1979) 

108. T. A. DeGrand, Y. G. Ng. and S. H. H. Tye, Jet Structure in e + e~ Annihilation as a 
Test of QCD and the Quark-Confining String, Phys. Rev. D 16, 3251 (1977) 

109. A. de Rujula, J. Ellis, E. G. Floratos and M. R. Gaillard, QCD Predictions for Hadronic 
Final States in e + e" Annihilation, Nucl. Phys. B 138, 387 (1978) 

110. A. Ore and J. L. Powell, Three Photon Annihilation of an Electron- Positron Pair, Phys. 
Rev. 75, 1696 (1949) 

111. Ch. Berger et al. [PLUTO Collaboration], Topology of the Upsilon Decay, Z. Phys. C8, 
101 (1981) 

112. B. R. Stella and H. J. Meyer, T(9.46 GeV) and the Gluon Discovery (A critical recol- 
lection of PLUTO results), Eur. Phys. J. H36, 203 (2011) 

113. R. Roller and H. Rrasemann, Excluding Scalar Gluons, Phys. Lett. B 88, 119 (1979) 

114. T. F. Walsh and P. M. Zerwas, Upsilon(9.46) Decays Do Test QCD, Phys. Lett. B 93, 
53 (1980) 

115. B. Niczyporuk et al. [LENA Collaboration], Charged Hadron Production in e + e~~ An- 
nihilation in the Upsilon and Upsilon' Region, Z. Phys. C 9, 1 (1981) 

116. R. Berkelman, Upsilon Spectroscopy at CESR, Phys. Rep. 98, 145 (1983) 

117. H. Albrecht et al. [ARGUS Collaboration], An Upper Limit for Two- Jet Production in 
Direct Upsilon(lS) Decays, Z. Phys. C 31, 181 (1986) 

118. G. Gidal, B. Armstrong and A. Rittenberg, "Major Detectors in Elementary Particle 
Physics," LBL 91-Supplement (March 1983; updated May 1985) 

119. H.L. Lynch, "Detectors for High Energy e + e~ Physics", in Ref. [54] 

120. H. Georgi and D. V. Nanopoulos, f-Quark Mass in a Superunified Theory, Phys. Lett. 
B 82, 392 (1979) 

121. H. Fritzsch, Quark Masses and Flavor Mixing, Nucl. Phys. B 155, 189 (1979) 



76 



122. R. Brandelik et al. [TASSO Collaboration], e+e' Annihilation at High- Energies and 
Search for the i-Quark Continuum Contribution, Z. Phys. C 4, 87 (1980) 

123. C. Berger et al. [PLUTO Collaboration], Search for a 'Top' Threshold in Hadronic e + e~ 
Annihilation at Energies between 22 GeV and 31.6 GeV, Phys. Lett. B 86, 413 (1979) 

124. N. Cabibbo, Unitary Symmetry and Leptonic Decays, Phys. Rev. Lett. 10, 531 (1963) 

125. M. Kobayashi and T. Maskawa, CP Violation in the Renormalizable Theory of Weak 
Interaction, Prog. Theor. Phys. 49, 652 (1973) 

126. [Tevatron Electroweak Working Group and CDF and DO Collaboration], Combina- 
tion of CDF and DO Results on the Mass of the Top Quark, FERMILAB-TM-2427- 
E, TEVEWWG-TOP-2009-03, CDF-NOTE-9717, DO-NOTE-5899, [arXiv:0903.2"503l [her>- 
ex]] (2009) 

127. A. Ali, J. G. Korner, J. Willrodt and G. Kramer, Indications of bb Production and 
Decay at Petra, Phys. Lett. B 83, 375 (1979) 

128. C. Berger et al. [PLUTO Collaboration], Cross-Sections and Event Topologies in e + e~ 
Annihilation at 13 GeV and 17 GeV Observed with the Pluto Detector, Phys. Lett. B 81, 
410 (1979) 

129. R. Brandelik et al. [TASSO Collaboration], Properties of Hadron final States in e + e~ 
Annihilation at 13 GeV and 17 GeV Center-of-Mass Energies, Phys. Lett. B 83, 261 (1979) 

130. A. Ali, J. G. Korner, G. Kramer and J. Willrodt, Jet-Like Distributions from the Weak 
Decay of Heavy Quarks, Z. Phys. C 1, 203 (1979) 

131. A. Ali, J. G. Korner, G. Kramer and J. Willrodt, Nonleptonic Weak Decays of Bottom 
Mesons, Z. Phys. C 1, 269 (1979) 

132. H. Schopper, First Results from PETRA Provide New Aspects of the Basic Structure 
of Matter, Naturwiss. 67, 161 (1980) 

133. S. L. Wu, Hadron Jets and the Discovery of the Gluon, in the Proceedings of 3rd 
International Symposium on the History of Particle Physics: The Rise of the Standard 
Model, Stanford, California, 2^-27 Jun 1992, p. 600 

134. J. G. Branson, Gluon Jets, Talk given at International Conference on the History of 
Original Ideas and Basic Discoveries in Particle Physics, Erice, Italy, 29 Jul - 4 Aug (1994) 

135. P. Soding, B. Wiik, G. Wolf and S. L. Wu, The First Evidence for Three- Jet Events in 
e+e' Collisions at PETRA: First Direct Observation of the Gluon, Talk given at Award 
Ceremony of the 1995 EPS High Energy and Particle Physics Prize. Brussels, Belgium, 27 
Jul - 2 Aug 1995, Report DESY-96-193 (1996) Published in Brussels EPS HEP 1995:3-14 
(QCD161:I48:1995) 

136. J. Ellis, Those were the Days: Discovering the Gluon, CERN Cour. 49N6, 15 (2009) 

137. P. Soding, On the Discovery of the Gluon, Eur. Phys. J. H 35, 3 (2010) 

138. G. C. Fox and S. Wolfram, Observables for the Analysis of Event Shapes in e+e' 
Annihilation and Other Processes, Phys. Rev. Lett. 41, 1581 (1978) 

139. J. R. Ellis and I. Karliner, Measuring the Spin of the Gluon in e+e' Annihilation, 
Nucl. Phys. B 148, 141 (1979) 

140. A. Ali, J. G. Korner, Z. Kunszt, J. Willrodt, G. Kramer, G. Schierholz and E. Pietari- 
nen, Four-Jet Production in e+e' Annihilation, Phys. Lett. B 82, 285 (1979) 

141. QCD Predictions for Four- Jet Final States in e+e' Annihilation, Nucl. Phys. B 167, 
454 (1980) 

142. R. K. Ellis, D. A. Ross and A. E. Terrano, Calculation of Event Shape Parameters in 
e+e' Annihilation, Phys. Rev. Lett. 45, 1226 (1980) 

143. R. K. Ellis, D. A. Ross and A. E. Terrano, The Perturbative Calculation of Jet Structure 
in e+e' Annihilation, Nucl. Phys. B 178, 421 (1981) 

144. J. A. M. Vermaseren, K. J. F. Gaemers and S. J. Oldham, Perturbative QCD Calcu- 
lation of Jet Cross-Sections in e+e~ Annihilation, Nucl. Phys. B 187, 302 (1981) 

145. R. K. Ellis and D. A. Ross, On the Thrust Distribution in e+e~ Annihilation, Phys. 
Lett. B 106, 88 (1981) 

146. Z. Kunszt, Magnitude of the 0{o? s ) Corrections to Jet Production in e+e' Annihilation, 
Phys. Lett. B 107, 123 (1981) 

147. L. Clavelli and D. Wyler, Kinematical Bounds on Jet Variables and the Heavy Jet 
Mass Distribution, Phys. Lett. B 103, 383 (1981) 



77 



148. A. Ali, The QCD Effective Coupling Constant in e+e Annihilation, Phys. Lett. B 
110, 67 (1982) 

149. R. Brandelik et al. [TASSO Collaboration], Comparison of e+e" Annihilation with 
QCD and Determination of the Strong Coupling Constant, Phys. Lett. B 94, 437 (1980) 

150. H. B. Newman, Results on Jets, QCD and Lepton Production from the Mark-J. (Talk), 
AIP Conf. Proc. 68, 627 (1981) 

151. C. L. Basham, L. S. Brown, S. D. Ellis and S. T. Love, Energy Correlations in Electron 
- Positron Annihilation: Testing QCD, Phys. Rev. Lett. 41, 1585 (1978); 

152. C. L. Basham, L. S. Brown, S. D. Ellis and S. T. Love, Energy Correlations in Electron- 
Positron Annihilation in Quantum Chromodynamics: Asymptotically Free Perturbation 
Theory, Phys. Rev. D 19, 2018 (1979) 

153. A. Ali and F. Barreiro, An 0(a 2 ) Calculation of Energy-Energy Correlation in e+e" 
Annihilation and Comparison with Experimental Data, Phys. Lett. B 118, 155 (1982 

154. A. Ali and F. Barreiro, Energy-Energy Correlations in e+e" Annihilation, Nucl. Phys. 
B 236, 269 (1984) 

155. F. Csikor, Quark Mass Effects for Energy-Energy Correlations in High-Energy e+e~ 
Annihilation, Phys. Rev. D 30, 28 (1984) 

156. K. H. Cho, S. K. Han and J. K. Kim, Energy-Energy Correlations in Electron Positron 
Annihilation: Z Boson and Heavy Quark Effects, Nucl. Phys. B 233, 161 (1984) 

157. D. G. Richards, W. J. Stirling and S. D. Ellis, Second Order Corrections to the Energy- 
Energy Correlation Function in Quantum Chromodynamics, Phys. Lett. B 119, 193 (1982) 

158. D. G. Richards, W. J. Stirling and S. D. Ellis, Energy-Energy Correlations to Second 
Order in Quantum Chromodynamics, Nucl. Phys. B 229, 317 (1983) 

159. W. Bartel et al. [JADE Collaboration], Measurements of Energy Correlations in 
e+e" -> Hadrons, Z. Phys. C 25, 231 (1984) 

160. B. Adcva et al. [MARK-J Collaboration], A Model Independent Second Order Deter- 
mination of the Strong Coupling Constant a s , Phys. Rev. Lett. 50, 2051 (1983) 

161. W. Braunschweig et al. [TASSO Collaboration], A Study Of Energy-Energy Correla- 
tions between 12 GeV and 46.8 GeV CM Energies, Z. Phys. C 36, 349 (1987) 

162. C. Berger et al. [PLUTO Collaboration], A Study of Energy- Energy Correlations in 
e+e" Annihilations at ^ = 34.6 GeV, Z. Phys. C 28, 365 (1985) 

163. K. Fabricius, I. Schmitt, G. Schierholz and G. Kramer, Order a 2 Correction to Jet 
Cross-Sections in e+e" Annihilation, Phys. Lett. B 97, 431 (1980) 

164. K. Fabricius, I. Schmitt, G. Schierholz and G. Kramer, Higher Order Perturbative QCD 
Calculation of Jet Cross-Sections in e+e" Annihilation, Z. Phys. C 11, 315 (1981) 

165. F. Gutbrod, G. Kramer and G. Schierholz, Higher Order QCD Corrections to the 
Three- Jet Cross-Sections: Bare Versus Dressed Jets, Z. Phys. C 21, 235 (1984) 

166. G. Kramer and B. Lampe, Jet Cross-Sections in e+e" Annihilation, Fortsch. Phys. 37, 
161 (1989) 

167. F. Gutbrod, G. Kramer, G. Rudolph and G. Schierholz, Recombination Dependence of 
the 0(a 2 s ) Three- Jet Cross-Section in e+e" Annihilation, Z. Phys. C 35, 543 (1987) 

168. S. Bethke et al. [JADE Collaboration], Experimental Investigation of the Energy De- 
pendence of the Strong Coupling Strength, Phys. Lett. B 213, 235 (1988) 

169. T. D. Gottschalk and M. P. Shatz, A Reassessment of the 0(a 2 ) Three- Jet Cross- 
Section for e+e" Annihilation, Phys. Lett. B 150, 451 (1985) 

170. S. Komamiya et al, Determination of a s from a Differential Jet Multiplicity Distribu- 
tion at SLC and PEP, Phys. Rev. Lett. 64, 987 (1990) 

171. W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Deep Inelastic Scattering 
Beyond the Leading Order in Asymptotically Free Gauge Theories, Phys. Rev. D 18, 3998 
(1978) 

172. R. Brock et al. [CTEQ Collaboration], Handbook of perturbative QCD: Version 1.0, 
Rev. Mod. Phys. 67, 157 (1995) 

173. S. Bethke et al, Studies of Jet Production Rates in e+e" Annihilation at E(cm) = 29 
GeV, Z. Phys. C 43, 325 (1989) 

174. S. Bethke, An Experimental Approach to Optimize and Test Perturbative QCD to 
0(a 2 ), Z. Phys. C 43, 331 (1989) 



78 



175. I. H. Park et al. [AMY Collaboration], Experimental Evidence for the Non-Abelian 
Nature of QCD, Phys. Rev. Lett. 62, 1713 (1989) 

176. W. Braunschweig et al. [TASSO Collaboration], Analysis of Multi-Jet Final States in 
e+e" Annihilation, Phys. Lett. B 214, 286 (1988) 

177. B. Lohr et al., Variation of the Strong Coupling Constant from a Measurement of the 
Jet Energy Spread in e+e" Annihilation, Phys. Lett. B 122, 90 (1983) 

178. H. J. Behrend et al. [CELLO Collaboration], Model Independent Limits on Aqcd from 
e+e" Annihilation in the Energy Range from 14 GeV to 46 GeV, Z. Phys. C 44, 63 (1989) 

179. W. Bartel et al. [JADE Collaboration], Experimental Study of Jets in Electron - 
Positron Annihilation, Phys. Lett. B 101, 129 (1981) 

180. W. Bartel et al. [JADE Collaboration], Test of Fragmentation Models by Comparison 
with Three- Jet Events Produced in e+e" -> Hadrons, Phys. Lett. B 134, 275 (1984) 

181. H. Aihara et al. [TPC/Two Gamma Collaboration], Tests of Models for Quark and 
Gluon Fragmentation in e+e" Annihilation at yfs = 29 GeV, Z. Phys. C 28, 31 (1985) 

182. M. Althoff et al. [TASSO Collaboration], A Study of Three- Jet Events in e+e" Anni- 
hilation into Hadrons at 34.6 GeV Center-of-Mass Energy, Z. Phys. C 29, 29 (1985) 

183. P. D. Sheldon et al, A Comparison of the Particle Flow in Three- Jet and Radiative 
Two- Jet Events from e+e" Annihilation at E(cm) = 29 GeV, Phys. Rev. Lett. 57, 1398 
(1986) 

184. G. Altarelli and G. Parisi, Asymptotic Freedom in Parton Language, Nucl. Phys. B 
126, 298 (1977) 

185. L. N. Lipatov, The Parton Model and Perturbation Theory, Sov. J. Nucl. Phys. 20, 94 
(1975) [Yad. Fiz. 20, 181 (1974)] 

186. V. N. Gribov and L. N. Lipatov, Deep Inelastic ep Scattering in Perturbation Theory, 
Sov. J. Nucl. Phys. 15, 438 (1972) [Yad. Fiz. 15, 781 (1972)] 

187. Y. L. Dokshitzer, Calculation of the Structure Functions for Deep Inelastic Scattering 
and e+e" Annihilation by Perturbation Theory in Quantum Chromodynamics, Sov. Phys. 
JETP 46, 641 (1977) [Zh. Eksp. Teor. Fiz. 73, 1216 (1977)] 

188. Y. I. Azimov, Y. L. Dokshitzer, V. A. Khoze and S. I. Troyan, Humpbacked QCD 
Plateau in Hadron Spectra, Z. Phys. C 31, 213 (1986) 

189. P. Abreu et al. [DELPHI Collaboration], Measurement of the Gluon Fragmentation 
Function and a Comparison of the Scaling Violation in Gluon and Quark Jets, Eur. Phys. 
J. C 13, 573 (2000) 

190. C. P. Fong and B. R. Webber, Higher Order QCD Corrections to Hadron Energy 
Distributions in Jets, Phys. Lett. B 229, 289 (1989) 

191. M. Z. Akrawy et al. [OPAL Collaboration], A Study of Coherence of Soft Gluons in 
Hadron Jets, Phys. Lett. B 247, 617 (1990) 

192. P. Abreu et al. [DELPHI Collaboration], The Scale-dependence of the Hadron Multi- 
plicity in Quark and Gluon Jets and a Precise Determination of C(A)/C(F), Phys. Lett. 
B 449, 383 (1999) [arXivihep ex/9903073| . 

1 93. I. M. Dremin and J . W. GaryT^Hadron multiplicities, Phys. Rept. 349, 301 (2001) 
[arXiv:hep-ph/0004215] 

194. Y. L. Dokshitzer, V. A. Khoze and S. I. Troian, Specific Features of Heavy Quark 
Production. 1. Leading Quarks, Leningrad Report LU-TP-92-10 (1992) 

195. C. Pahl, S. Bethke, S. Kluth, J. Schieck and the Jade collaboration, Study of Moments 
of Event Shapes and a Determination of a s using e+e" Annihilation Data from Jade, Eur. 
Phys. J. C 60, 181 (2009) [Erratum-ibid. C 62, 451 (2009)] [arXiv:0810.2933l [hep-ex]] 

196. P. Zerwas, W & Z Physics at LEP, Eur. Phys. J. C 34, 41 (2004) and in Prestigious 
discoveries at CERN, ed. R. Cashmore, G. Maiani and J. -P. Revol, Springer Verlag, (2004) 

197. P. Abreu et al. [DELPHI Collaboration], Energy-dependence of Event Shapes and of 
a s at LEP2, Phys. Lett. B 456, 322 (1999) 

198. P. Pfeifenschneider et al. [JADE collaboration and OPAL Collaboration], QCD Anal- 
yses and Determinations of a s in e+e" Annihilation at Energies between 35 GeV and 189 
GeV, Eur. Phys. J. C 17, 19 (2000) 

199. P. Achard et al. [L3 Collaboration], Determination of a s from Hadronic Event Shapes 
in e+e" Annihilation at 192 GeV < </s < 208 GeV, Phys. Lett. B 536, 217 (2002) 



79 



200. A. Heister et al. [ALEPH Collaboration], Studies of QCD at e+e~ Centre-of-Mass 
Energies between 91 GeV and 209 GeV, Eur. Phys. J. C 35, 457 (2004) 

201. S. Catani, G. Turnock, B. R. Webber and L. Trentadue, Thrust Distribution in e + e~ 
Annihilation," Phys. Lett. B 263, 491 (1991) 

202. S. Catani, G. Turnock and B. R. Webber, Heavy Jet Mass Distribution in e + e~ anni- 
hilation," Phys. Lett. B 272, 368 (1991) 

203. S. Catani, L. Trentadue, G. Turnock and B. R. Webber, Resummation of Large Loga- 
rithms in e+e" Event Shape Distributions, Nucl. Phys. B 407, 3 (1993) 

204. L. Lonnblad, Ariadne Version 4: A Program for Simulation of QCD Cascades Imple- 
menting the Color Dipole Model, Comput. Phys. Commun. 71, 15 (1992) 

205. Y. L. Dokshitzer and B. R. Webber, Calculation of Power Corrections to Hadronic 
Event Shapes, Phys. Lett. B 352, 451 (1995) [arXiv:hep- ph/9504219] . 

206. Y. L. Dokshitzer, A. Lucenti, G. Marchesini and G. P. Salam, Universality of 1/Q 
Corrections to Jet-shape Observables Rescued, Nucl. Phys. B 511, 396 (1998) [Erratum- 
ibid. B 593, 729 (2001)] 

207. Y. L. Dokshitzer, A. Lucenti, G. Marchesini and G. P. Salam, On the Universality of 
the Milan Factor for 1/Q Power Corrections to Jet Shapes, JHEP 9805, 003 (1998) 

208. M. Dasgupta and G. P. Salam, Event shapes in e+ e- annihilation and deep inelastic 
scattering, J. Phys. G 30, R143 (2004) [arXiv:hep-ph/0312283| 

209. G. Dissertori, A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, G. Heinrich 
and H. Stenzel, Precise Determination of the Strong Coupling Constant at NNLO in QCD 
from the Three- Jet Rate in Electron-Positron Annihilation at LEP," Phys. Rev. Lett. 104, 
072002 (2010) 

210. A. Signer and L. J. Dixon, Electron Positron Annihilation into Four Jets at Next-to- 
Leading Order in a s , Phys. Rev. Lett. 78, 811 (1997) 

211. L. J. Dixon and A. Signer, Complete O(al) Results for e + e~ — > ("f,Z) — > Four Jets, 
Phys. Rev. D 56, 4031 (1997) 

212. Z. Nagy and Z. Trocsanyi, Next-to-Leading Order Calculation of Four- Jet Observables 
in Electron Positron Annihilation, Phys. Rev. D 59, 014020 (1999) [Erratum-ibid. D 62, 
099902 (2000)] 

213. P. D. Acton et al. [OPAL Collaboration], A Global Determination of a s (Mz) at LEP, 
Z. Phys. C 55, 1 (1992) 

214. R. Frederix, S. Frixione, K. Melnikov and G. Zanderighi, NLO QCD Corrections to 
Five- Jet Production at LEP and the Extraction of a s (M z ), JHEP 11:050 (2010) 

215. T. Sjostrand, S. Mrenna and P. Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 
0605, 026 (2006) 

216. G. Corcella et al, HERWIG 6.5: An Event Generator for Hadron Emission Reactions 
with Interfering Gluons (including Supersymmetric Processes), JHEP 0101, 010 (2001) 

217. S. Catani, F. Krauss, R. Kuhn and B. R. Webber, QCD Matrix Elements + Parton 
Showers, JHEP 0111, 063 (2001) 

218. M. Bengtsson and P. M. Zerwas, Four- Jet Events in e + e~ Annihilation: Testing the 
Three-Gluon Vertex, Phys. Lett. B 208, 306 (1988) 

219. O. Nachtmann and A. Reiter, A Test for the Gluon Self Coupling in the Reactions 
e + e" -> Four Jets and Z -> Four Jets, Z. Phys. C 16, 45 (1982) 

220. J. G. Korner, G. Schierholz and J. Willrodt, QCD Predictions for Four- Jet Final States 
in e + e~ Annihilation: 2. Angular Correlations as a Test of the Triple Gluon Coupling, Nucl. 
Phys. B 185, 365 (1981) 

221. Z. Nagy and Z. Trocsanyi, Four- Jet Angular Distributions and Color Charge Mea- 
surements: Leading Order versus Next-to-Leading Order, Phys. Rev. D 57, 5793 (1998) 
[arXiv:hep-ph/9712385l 

222. B. Adeva et al. ]L3 Collaboration], A Test of QCD Based on Four- Jet Events from Z 
Decays, Phys. Lett. B 248, 227 (1990) 

223. G. Abbiendi et al. [OPAL Collaboration], A Simultaneous Measurement of the QCD 
Colour Factors and the Strong Coupling, Eur. Phys. J. C 20, 601 (2001) 

224. M. Czakon, The Four-Loop QCD Beta-function and Anomalous Dimensions, Nucl. 
Phys. B 710, 485 (2005) 



80 



225. T. van Ritbergen, J. A. M. Vermaseren and S. A. Larin, The Four-Loop Beta Function 
in Quantum Chromodynamics, Phys. Lett. B 400, 379 (1997) 

226. T. Gehrmann, M. Jaquier and G. Luisoni, Hadronization Effects in Event Shape Mo- 
ments, Eur. Phys. J. C 67,57 (2010) 

227. G. Dissertori, A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, G. Heinrich, 
G. Luisoni and H. Stenzel, Determination of the Strong Coupling Constant Using Matched 
NNLO+NLLA Predictions for Hadronic Event Shapes in e + e" Annihilations, JHEP 0908, 
036 (2009) 

228. S. Bethke, The 2009 Wolrd Average of a a (M z ), Eur. Phys. J. C 64, 689 (2009) 

229. J. C. Coll ins, D. E. S oper, G. F. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1-91 
(1988) hep-ph/0409313; 

230. We follow here the review by G. Dissertori and G.P. Salam on "Quantum Chromody- 
namics" in Ref. [7] 

231. P.M. Nadolsky et al., Implications of CTEQ Global Analysis for Collider Observables, 
Phys. Rev. D 78, 013004 (2008), and earlier references given therein 

232. A.D. Martin et al, Parton Distributions for the LHC, Eur. Phys. J. C 63, 189 (2009), 
and earlier papers referenced therein 

233. S. Catani and M.H. Seymour, A General Algorithm for Calculating Jet Cross-Sections 
in NLO QCD, Nucl. Phys. B 485, 291 (1997); [Erratum: Nucl. Phys. B 510, 503 (1998)] 

234. Z. Nagy and Z. Troscanyi, Multijet Cross-Sections in Deep Inelastic Scattering at Next- 
to-Leading Order, Phys. Rev. Lett. 87, 082001 (2001) 

235. C. Amsler et al. (Particle Data Group), Review of Particle Physics, Phys. Lett. B 667, 
1 (2008) 

236. F.D. Aaron et al. [HI Collaborate] , Jet Production in ep Collisions at High Q 2 and 
Determination of a s , Eur. Phys. J. C 65, 363 (2010) 

237. S. Chekanov et al. [ZEUS Collaboration], Inclusive Jet Cross Sections in the Breit frame 
in Neutral Current Deep Inelastic Scattering at HERA and Determination of a s , Phys. 
Lett. B 547, 164 (2002) 

238. S. Chekanov et al. [ZEUS Collaboration], Inclusive- Jet and Dijet Cross Sections in 
Deep Inelastic Scattering at HERA, Nucl. Phys. B 765, 1 (2007) 

239. A. Aktas et al. [HI Collaboration], Measurement of Inclusive Jet Production in Deep- 
Inelastic Scattering at High Q 2 and Determination of the Strong Coupling, Phys. Lett. B 
653, 134 (2007) 

240. S. Chekanov et al. [ZEUS Collaboration], Jet- Radius Dependence of Inclusive- Jet 
Cross-Sections in Deep Inelastic Scattering at HERA, Phys. Lett. B 649, 12 (2007) 

241. C. F. von Weizsacker, Radiation Emitted in Collisions of Very Fast Electrons (in Ger- 
man), Z. Phys. 88, 612 (1934) 

242. E.J. Williams, Nature of the High- Energy Particles of Penetrating Radiation and Status 
of Ionization and Radiation Formulae, Phys. Rev. 45, 729 (1934) 

243. S. Frixione et al., Improving the Weizsacker- Williams Approximation in Electron - 
Proton Collisions, Phys. Lett. B 319, 339 (1993) 

244. C. H. Llewellyn Smith, QCD Predictions for Processes Involving Real Photons, Phys. 
Lett. B 79, 83 (1978) 

245. S. J. Brodsky, T. A. DeGrand, J. F. Gunion and J. H. Weis, Production of Large 
Transverse Momentum Jets in Photon - Photon Collisions, Phys. Rev. Lett. 41, 672 (1978) 

246. M. Klasen and G. Kramer, Inclusive Two- Jet Production at HERA: Direct and Re- 
solved Cross-Sections in Next-to-Leading Order QCD, Z. Phys. C 76, 67 (1997 

247. M. Klasen, T. Kleinwort and G. Kramer, Inclusive Jet Production in 7p and 77 Pro- 
cesses: Direct and Resolved Photon Cross-Sections in Next-to-Leading Order QCD, Eur. 
Phys. J. direct C 1, 1 (1998) 

248. B.W. Harris and J.F. Owen, Jet Photoproduction and the Structure of the Photon, 
Phys. Rev. D 57, 5555 (1998) 

249. S. Frixione and G. Ridolfi, Jet Photoproduction at HERA, Nucl. Phys. B507, 315 
(1997) 

250. A. Aktas et al. [HI Collaboration], Photoproduction of Dijets with High Transverse 
Momenta at HERA, Phys. Lett. B 639, 21 (2006) 



81 



251. S. Chekanov et al. [ZEUS Collaboration], Eigh-E T Dijet Photoproduction at HERA, 
Phys. Rev. D 76, 072011 (2007) 

252. S. Chekanov et al. [ZEUS Collaboration], High Mass Dijet Cross-Sections in Photopro- 
duction at HERA, Phys. Lett. B 531, 9 (2002) 

253. S. Chekanov et al. [ZEUS Collaboration], Dijet Photoproduction at HERA and the 
Structure of the Photon, Eur. Phys. J. C 23, 615 (2002) 

254. P. Aurenche, M. Fontannaz and J. P. Guillet, New NLO Parametrizations of the Parton 
Distributions in Real Photons, Eur. Phys. J. C 44, 395 (2005) 

255. F. Cornet, P. Jankowski and M. Krawczyk, A New 5 Flavour NLO Analysis and 
Parametrizations of Parton Distributions of the Real Photon, Phys. Rev. D 70, 093004 
(2004) 

256. S. Chekanov et al. [ZEUS Collaboration], Scaling Violations and Determination of a 3 
from Jet Production in jp Interactions at HERA, Phys. Lett. B 560, 7 (2003) 

257. P. Aurenche, J. P. Guillet and M. Fontannaz, Parton Distributions in the Photon, Z. 
Phys. C 64, 621 (1994) 

258. M. Chick, E. Reya and A. Vogt, Parton Structure of the Photon Beyond the Leading 
Order, Phys. Rev. D 45, 3896 (1992) 

259. M. Cluck, E. Reya and A. Vogt, Photonic Parton Distributions, Phys. Rev. D 46, 1973 
(1992) 

260. W. Slominski, H. Abramowicz and A. Levy, NLO Photon Parton Parametrization Using 
e+e" and ep Data, Eur. Phys. J. C 45, 633 (2006) 

261. V. S. Fadin, E. A. Kuraev and L. N. Lipatov, On the Pomeranchuk Singularity in 
Asymptotically Free Theories, Phys. Lett. B 60, 50 (1975) 

262. I. I. Balitsky and L. N. Lipatov, The Pomeranchuk Singularity in Quantum Chromo- 
dynamics, Sov. J. Nucl. Phys. 28, 822 (1978) [Yad. Fiz. 28, 1597 (1978)] 

263. M. Dittmar et al., Summary Report for the HERA - LHC Workshop Proceedings: 
Parton Distributions, Report FUM-936-FT, |arXiv:0901.2504l [hep-ph] (2009) 

264. L. McLerran, Gluon Evolution and Saturation Proceedings, presented at the Gmbov 
Memorial Workshop on Quantum Chromodynamics and Beyond, |arXiv:1011. 3202 [hep- 
ph]] (2010) 

265. For the concept and definitions of the terms Color Glass Condensate, see, L. McLerran, 
The CGC and the Glasma: Two Lectures at the Yukawa Insitute, arXiv:1011.3204 [hep- 
ph]] (2010) 

266. See, for a review, E. Iancu and R. Venugopalan, The Color Glass Condensate and High 
Energy Scattering in QCD, in Quark gluon plasma, edited by R. Hwa and X.N. Wang, 
World Scientific, 2003, [arXivihep ph /0303204| (2003) 

267. S. J. Brodsky, F. E. Close and J. F. Gunion, Phenomenology of Photon Processes, 
Vector Dominance and Crucial Tests for Parton Models, Phys. Rev. D6, 177 (1972) 

268. J. Kwiecinski, B. M. Badelek, Analysis of the Electroproduction Structure Functions 
in the Low Q**2 Region combining the Vector Meson Dominance and the Parton Model 
with possible Scaling Violation," Z. Phys. C43, 251 (1989) 

269. G. Abbiendi et al. [OPAL Collaboration], Di-jet Production in Photon Photon Colli- 
sions at y/a ee = from 189 GeV to 209 GeV, Eur. Phys. J. C 31, 307 (2003) 

270. J. Abdallah et al. [DELPHI Collaboration], Di-jet Production in Gamma-Gamma Col- 
lisions at LEP2, Eur. Phys. J. C 58, 531 (2008) 

271. P. Achard et al. [L3 Collaboration], Inclusive Jet Production in Two-Photon Collisions 
at LEP, Phys. Lett. B 602, 157 (2004) 

272. G. Abbiendi et al. [OPAL Collaboration], Inclusive Jet Production in Photon-Photon 
Collisions at ^s ee = from 189 to 209 GeV, Phys. Lett. B 658, 185 (2008) 

273. V. M. Abazov et al. [DO Collaboration], Measurement of the Inclusive Jet Cross-Section 
in pp Collisions at y/s = 1.96 TeV, Phys. Rev. Lett. 101, 062001 (2008) 

274. Measurement of the Inclusive Jet Cross Section in pp Collisions at 7 TeV, CMS Physics 
Analysis Summary, CMS PAS QCD-10-011 (2010) 

275. A. Bhatti and D. Lincoln, Jet Physics at the Tevatron, FERMILAB-PUB- 10-028- E- 
PPD [arXiv:1002.1708l [hep-ex]] (2010) 



82 



276. T. Aaltonen et al. [CDF Collaboration], Measurement of the Inclusive Jet Cross Section 
at the Fermilab Tevatron pp Collider Using a Cone-Based Jet Algorithm, Phys. Rev. D 
78, 052006 (2008) [Erratum- ibid. D 79, 119902 (2009] 

277. A. Abulencia et al. [CDF - Run II Collaboration], Measurement of the Inclusive Jet 
Cross Section Using the k,T Algorithm in pp Collisions at y/s — 1.96 TeV with the CDF 
II Detector, Phys. Rev. D 75, 092006 (2007) [Erratum-ibid. D 75, 119901 (2007)] 

278. S. D. Ellis, J. Huston, K. Hatakeyama, P. Loch and M. Tonnesmann, Jets in hadron- 
hadron collisions, Prog. Part. Nucl. Phys. 60, 484 (2008) |arXiv: 0712. 24471 [hep-ph]] 

279. Z. Nagy, Three-jet Cross Sections in Hadron Hadron Collisions at Next-to-leading Or- 
der, Phys. Rev. Lett. 88, 122003 (2002) 

280. M. Bahr et al, Herwig++ Physics and Manual, Eur. Phys. J. C 58, 639 (2008) 

281. A. Buckley et aL larXivTlOl. 25991 [hep-ph], 

282. E. Eichten, K. D. Lane and M. E. Peskin, New Tests for Quark and Lepton Substruc- 
ture, Phys. Rev. Lett. 50, 811 (1983) 

283. U. Baur, I. Hinchliffe and D. Zeppenfeld, Excited Quark Production at Hadron Collid- 
ers, Int. J. Mod. Phys. A 2, 1285 (1987) 

284. U. Baur, M. Spira and P. M. Zerwas, Excited Quark and Lepton Production at Hadron 
Colliders, Phys. Rev. D 42, 815 (1990) 

285. P. H. Frampton and S. L. Glashow, Chiral Color: An Alternative to the Standard 
Model, Phys. Lett. B 190, 157 (1987) 

286. J. Bagger, C. Schmidt and S. King, Axigluon Production in Hadronic Collisions, Phys. 
Rev. D 37, 1188 (1988) 

287. R. S. Chivukula, A. G. Cohen and E. H. Simmons, New Strong Interactions at the 
Tevatron?, Phys. Lett. B 380, 92 (1996) 

288. E. H. Simmons, Coloron phenomenology, Phys. Rev. D 55, 1678 (1997) 

289. K. D. Lane and M. V. Ramana, Walking Technicolor Signatures at Hadron Colliders, 
Phys. Rev. D 44, 2678 (1991) 

290. K. Lane and S. Mrenna, The Collider Phenomenology of Technihadrons in the Techni- 
color Straw Man Model, Phys. Rev. D 67, 115011 (2003) 

291. R. Foadi, M. T. Frandsen, T. A. Ryttov and F. Sannino, Minimal Walking Technicolor: 
Set Up for Collider Physics, Phys. Rev. D 76, 055005 (2007) 

292. A. Belyaev et al, Technicolor Walks at the LHC, Phys. Rev., D 79, 035006 (2009) 

293. T. Aaltonen et al. [CDF Collaboration], Search for New Particles Decaying into Dijets 
in Proton-Antiproton Collisions at yfs = 1.96 TeV, Phys. Rev. D 79, 112002 (2009) 

294. G. Aad et al. [ATLAS Collaboration], Search for New Particles in Two- Jet Final States 
in 7 TeV Proton-Proton Collisions with the ATLAS Detector at the LHC, Phys. Rev. Lett. 
105, 161801 (2010) 

295. V. Khachatryan et al. [CMS Collaboration], Search for Dijet Resonances in 7 TeV pp 
Collisions at CMS, Phys. Rev. Lett. 105, 211801 (2010) 

296. R. K. Ellis, G. Martinelli and R. Petronzio, Lepton Pair Production at Large Transverse 
Momentum in Second Order QCD, Nucl. Phys. B 211, 106 (1983) 

297. P. B. Arnold and M. H. Reno, The Complete Computation of High p t W and Z 
Production in 2nd Order QCD, Nucl. Phys. B 319, 37 (1989) [Erratum-ibid. B 330, 284 
(1990)] 

298. P. B. Arnold, R. K. Ellis and M. H. Reno, High p t W and Z Production at the Tevatron, 
Phys. Rev. D 40, 912 (1989) 

299. W. T. Giele, E. W. N. Glover and D. A. Kosower, Higher Order Corrections to Jet 
Cross-Sections in Hadron Colliders, Nucl. Phys. B 403, 633 (1993) 

300. J. M. Campbell and R. K. Ellis, Next-to-Leading Order Corrections to W+ 2Jet and 
Z+ 2Jet Production at Hadron Colliders, Phys. Rev. D 65, 113007 (2002) 

301. J. M. Campbell, R. K. Ellis and D. L. Rainwater, Next-to-Leading Order QCD Predic- 
tions for W+ 2Jet and Z+ 2Jet Production at the CERN LHC, Phys. Rev. D 68, 094021 
(2003) 

302. F. A. Berends, W. T. Giele, H. Kuijf, R. Kleiss and W. J. Stirling, Multi - Jet Produc- 
tion in W, Z Events at pp Colliders, Phys. Lett. B 224, 237 (1989) 



83 



303. F. A. Berends, H. Kuijf, B. Tausk and W. T. Giele, On the Production of a W and 
Jets at Hadron Colliders, Nucl. Phys. B 357, 32 (1991) 

304. J. M. Campbell and R. K. Ellis, MCFM for the Tevatron and the LH C, Nucl. Phys. 
Proc. Suppl. 205-206, 10 (2010); MCFM homepage, [http://mcfm.fnal.gov| 

305. K. Melnikov and F. Petriello, Electroweak Gauge Boson Production at Hadron Colliders 
through 0(o? a ), Phys. Rev. D 74, 114017 (2006) 

306. S. Catani, L. Cieri, G. Ferrera, D. de Florian and M. Grazzini, Vector Boson Production 
at Hadron Colliders: A Fully Exclusive QCD Calculation at NNLO, Phys. Rev. Lett. 103, 
082001 (2009) 

307. G. A. Ladinsky and C. P. Yuan, The Nonperturbative Regime in QCD Resummation 
for Gauge Boson Production at Hadron Colliders, Phys. Rev. D 50, 4239 (1994) 

308. C. Balazs and C. P. Yuan, Soft Gluon Effects on Lepton Pairs at Hadron Colliders, 
Phys. Rev. D 56, 5558 (1997) 

309. V. M. Abazov et al. [DO Collaboration], Measurement of the Normalized Z/'y* — > 
Transverse Momentum Distribution in pp Collisions at y/s — 1.96 TeV, Phys. Lett. B 693, 
522 (2010) 

310. Measurement of Z/^f* — > fi + + Jets Production Cross Section, CDF Note 10216, 
July 10, 2010 (CDF Collaboration); URL: |http: / / www-cdf-fnaJ.gov 

311. For a recent discussion of the standard electroweak model, see the review by J. Erler 
and P. Langacker: "Electroweak model and constraints on new physics" in Ref. \7\. 

312. D. E. Acosta et al. [CDF Collaboration], Measurement of the W + W~ Production Cross 
Section in pp Collisions at %/s = 1.96 TeV using Dilepton Events, Phys. Rev. Lett. 94, 
211801 (2005) 

313. T. Aaltonen et al. [CDF Collaboration], First Measurement of ZZ Production in pp 
Collisions at yfs = 1.96 TeV, ibid 100, 201801 (2008) 

314. V. M. Abazov et al. [DO Collaboration], Measurement of the WW Production Cross 
Section in pp Collisions at </s = 1.96 TeV, Phys. Rev. Lett. 94, 151801 (2005) [Erratum- 
ibid. 100, 139901 (2008)] 

315. V. M. Abazov et al. [DO Collaboration], Evidence of WW + WZ Production with 
Lepton + Jets Final States in Proton- Antiproton Collisions at ^/s =1.96 TeV, Phys. Rev. 
Lett. 102, 161801 (2009) 

316. T. Aaltonen et al. [CDF Collaboration], First Observation of Vector Boson Pairs in a 
Hadronic Final State at the Tevatron Collider, Phys. Rev. Lett. 103, 091803 (2009) 

317. J. M. Campbell and R. K. Ellis, An Update on Vector Boson Pair Production at Hadron 
Colliders, Phys. Rev. D 60, 113006 (1999) 

318. T. Aaltonen et al. [The CDF Collaboration], Measurement of the Top Quark Mass and 
pp — > ti Cross Section in the All-Hadronic Mode with the CDFII Detector, Phys. Rev. D 
81, 052011 (2010) 

319. A. Djouadi, M. Spira and P. M. Zerwas, QCD Corrections to Hadronic Higgs Decays, 
Z. Phys. C 70, 427 (1996) 

320. J. M. Butterworth, A. R. Davison, M. Rubin and G. P. Salam, Jet Substructure as a 
New Higgs Search Channel at the LHC, Phys. Rev. Lett. 100 (2008) 242001 

321. T. Aaltonen et al. [CDF Collaboration], Search for Standard Model Higgs Boson Pro- 
duction in Association with a W Boson at CDF, Phys. Rev. D 78, 032008 (2008) 

322. T. Aaltonen et al. [The CDF Collaboration], Search for Doubly Charged Higgs Bosons 
with Lepton-Flavor- Violating Decays involving Tau Leptons, Phys. Rev. Lett. 101, 121801 
(2008) 

323. T. Aaltonen et al. [CDF Collaboration], A Search for the Associated Production of the 
Standard-Model Higgs Boson in the All-Hadronic Channel, Phys. Rev. Lett. 103, 221801 
(2009) 

324. T. Plehn, G. P. Salam and M. Spannowsky, Fat Jets for a Light Higgs, Phys. Rev. Lett. 
104, 111801 (2010) 

325. J. M. Butterworth, J. R. Ellis and A. R. Raklev, Reconstructing Sparticle Mass Spectra 
using Hadronic Decays, JHEP 0705 (2007) 033 

326. J. M. Butterworth, B. E. Cox and J. R. Forshaw, WW Scattering at the CERN LHC, 
Phys. Rev. D 65, 096014 (2002) 



84 



327. U. Baur and L. H. Orr, Searching for tt Resonances at the Large Hadron Collider, 
Phys. Rev. D 77, 114001 (2008) 

328. S. Moch, J. A. M. Vermaseren and A. Vogt, The Three-loop Splitting Functions in 
QCD: The Non-Singlet Case, Nucl. Phys. B 688, 101 (2004) 

329. A. Vogt, S. Moch and J. A. M. Vermaseren, The Three-loop Splitting Functions in 
QCD: The Singlet Case," Nucl. Phys. B 691, 129 (2004) 



